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I .  Int roduct Ion 


1,  One  of  the  main  subjects  of  Analysis  is  the  integra¬ 
tion  of  differential  equations  ,  ordinary  and  partial  .  However, 
except  for  a  few  simple  cases  ,  the  integri  vioii  is  very  diffi¬ 
cult  . 

The  theory  of  analytic  functions  in  one  or  more  complex 
variables  ,  initiated  by  Cauchy  ,  Weierstrass  and  Riemann  was 
applied  by  them  to  the  study  of  differential  equations. 

Subsequent  researches  gave  important  results  which  will 
be  described  very  briefly  . 

Consider  a  system  of  ordinary  differential  equations 
(1)  y  =  f(x,y)  , 

where  y  and  f  are  vectors  in  an  n-dimenslonal  space  .  Suppose 
that  the  components  of  f  are  holomorphic  in  a  neighborhood 
of  X  =  x^  )  y  -  ;  then  ,  the  differential  system  (1)  has 

one  and  only  one  solution  y(x)  =  y(x;y^,x^)  such  that  y(x^)=*y^. 
This  function  y(x)  may  be  continued  analytically  and  gives 
rise  to  an  analytic  function  of  x  ,  the  singular  points  of 
which  have  to  be  determined  . 

When  (!)  is  linear  in  y  ,  the  results  are  as  follows  ; 
i.  the  singular  points  of  y(x)  are  among  the  singular  points 
of  the  coefficients  of  the  system  ;  thus  they  are  fixed  and 
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determined  directly  by  the  equations  ; 

ii.  the  set  of  all  solutions  of  (l)  is  a  linear  vector  space 
so  that  the  dependence  of  y(x)  on  the  initial  values  is  known 

Suppose  that  x  =  is  a  regular  singular  point  in  the 
sense  of  L,  Fuchs  .  Then  the  analytical  character  of  the  inte 
gral  in  the  neighborhood  of  x  =  x^  is  completely  determined 
by  the  indicial  equation  and  accordingly  is  given  by  a  set 
of  constants  vmich  are  easily  obtained  from  the  given  system. 

2.  More  complicated  circumstances  arise  when  (1)  is  not 
linear  in  y  ,  In  general  ,  it  is  not  possible  to  recognize 
if  a  given  value  x^  of  x  is  regular  or  singular  for  y(x)  ; 
the  analytical  character  of  y(x)  at  x^  may  depend  on  the 
value  y^  of  y(x^)  .  Moreover  ,  the  singular  points  of  y(x) 
may  be  parametric  (or  movable)  i.e.  may  depend  on  y^  ;  they 
may  also  be  algebraic  or  transcendental  ,  or  even  essential 
points  and  may  be  isolated  or  not  . 

Therefore  ,  the  following  problem  initiated  by  L.  Fuchs 
as  a  consequence  of  his  researches  on  linear  differential 
equations  is  of  importance  :  to  determine  all  the  equations 
(2)  R(  y  ,  y  ,  x)  =  0  , 

where  R  is  a  polynomial  in  y  and  y  with  analytic  coefficients 
in  X  ,  whose  integral  has  no  parametric  critical  points  .  In 
the  neighborhood  of  every  movable  singularity  ,  y(x)  must  be 
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one-valued  .  To  abb.  viate  ,  we  shall  say  that  an  equation 
of  this  type  is  stable  together  with  its  integrals  . 

This  problem  was  considered  by  Abel  and  Jacobi  when 
R=y  -(l-y)(l-ky)  , 

(k  a  constant  )  and  by  Briot  et  Bouquet  when  R  is  Independent 
of  X  . 

The  solution  of  this  problem  was  given  by  L,  Fuchs  , 
Poincare  and  Painleve  and  is  as  follows  :  the  Integral  of 
the  stable  equations  (2)  are  determined  algebraically  ,  or 
by  quadratures  ,  or  depend  on  a  Riccati  equation  .  Consequently, 
they  are  reducible  to  classical  transcendents  ,  i.e,  to 
algebraic  and  elliptic  functions  or  to  functions  defined  by 
linear  differential  equations  . 

3.  To  define  new  transcendental  functions  ,  one  must 

therefore  consider  differential  equations  of  higher  order  , 

However  ,  new  complications  appear  ,  of  which  the  most  serious 

is  the  possible  existence  of  parametric  transcendental  or 

even  essential  singular  points  .  For  instance  ,  in  the  latter 

case  ,  y(x)  comes  arbitrarily  close  to  any  complex  value  in 

every  neighborhood  of  x  so  that  lim  y(x)  is  not  uniquely 

x=x 

o 

determined  , 

Few  and  Inconclusive  results  were  obtained  by  Picard  and 
Mlttag-Leffler  concerning  stable  equations  of  the  second  order. 
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Palnleve  was  the  first  to  attack  successfully  this  problem 
and  to  overcome  the  difficulties  .  To  this  effect  ,  Painleve 
has 

i.  to  obtain  a  set  of  necessary  conditions  for  the  absence  of 
parametric  critical  points  ; 

ii.  to  show  that  the  necessary  conditions  thus  obtained  are 
or  are  not  sufficient  . 

Consider  for  instance  the  equation 
(3)  y  =  R(  y  ,  y  »  X  ) 

where  R  is  a  rational  function  of  y  and  y  ,  with  analytic 
coefficients  in  x  ,  In  order  to  obtain  a  set  of  necessary 
conditions  for  stability  ,  Painleve  introduces  in  (3)  by  a 
suitable  transformation  of  y  and  x  ,  a  parameter  e  in  such 
a  way  that  the  new  equation  has  the  same  fixed  critical  points 
as  (3)  and  is  integrable  when  e  =  0  .  Then  y(x,e)  may  be  de¬ 
veloped  into  a  series  of  ascending  powers  of  e  , 

p 

y{x,e)  =  y^Cx)  +  e  y^(x)  +  e  yg(x)  +  ...  , 

the  coefficients  of  which  are  also  uniform  and  determined  by 
quadratures  ;  the  conditions  that  these  functions  be  stable 
are  necessary  conditions  for  the  stability  of  (3)  . 

This  method  was  used  by  Painleve  and  Gambler  to  determine 
all  the  stable  equations  (2)  .  The  result  is  as  follows  : 

The  stable  equations  (2)  are  Integrable  in  terms  of  classi¬ 
cal  transcendents  or  may  be  reduced  by  a  transformation 


5. 


(4) 


ay  +  b 
u  =  ^ 


cy  +  d 


z  =  qp  (x) 


where  a,b,c,d,  tp  are  analytic  functions  of  x  ,  to  one  of  the 
six  canonical  equations  in  the  following  table  : 

Table  I 


»i  p 

11.  y  =  6y''  +  x  , 

II  o 

12.  y  =  2y-^+xy+a 


I  5. 


*2 


13.  y  = 


14.  y  = 


_  y 


y 

2y 


+  «y  +  P  +  Y  y3  +  , 

XX  y 


+  ^  y^+4xy^+2(x^-  a)y+  P  , 


**  /  1  1 

y=  __  +  ___ 

V  2y  y  -  1 


*2  y  , (y-1) 


2/ 


y--  “ 

X 


ay+_l_ 
y 


> 


y  (y+i) 


y  _ 

X  y-1 


,  II 

I  6.  y 

1 

f  W  b  +  1 

2 

Vy  y-1  y-xy 

X  x-l  x-y  j 

+ 

y(y-i: 

1  (y-x)  [*  +  P  X  y  (x-l) 

6x(x-l)  , 

2/ 

X  (x- 

«\2  L  2  /  i\2 

■1)  y  (y-1) 

,  \2  J 

(y  -  x) 

a>  P  » 

T  »  6 

are  constants  . 

As  was  shown  by  Painleve  ,  the  integrals  of  the  equations 
of  Table  I  are  not  reducible  to  classical  transcendents 
(except  possibly  for  particular  values  of  a  ,  (5  ,  y  ,  ft  ) 
and  thus  define  new  transcendental  functions  . 
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2+.  The  method  xised  by  Painleve  ,  although  theoretically 
simple  ,  implies  intricate  calculations  .  Because  of  these 
difficulties  and  because  of  the  great  number  of  cases  that 
have  to  be  considered  ,  Painleve  found  only  three  equations 
of  Table  I  (  namely  I  1*2,3)  ; the  three  other  equations  were 
discovered  by  Gambler  . 

The  method  of  Painleve  may  also  be  applied  to  equations 
of  order  higher  than  two  ;  however  ,  the  intricacies  increase 
with  the  order  of  the  equation  . 

To  reduce  these  difficulties  to  a  minimum  ,  another  method 
will  be  developped  in  this  paper  ;  it  will  systematically  avoid, 
as  far  as  possible  ,  any  integration  of  systems  of  linear 
differential  equations  and  will  reduce  the  problem  to  elemen¬ 
tary  algebraic  processes  ,  This  method  was  initiated  by  the 
author  in  1939  ;  a  few  examples  were  given  at  that  time  . 

The  ultimate  purpose  of  the  actual  theory  is  :  to  find 
out  if  a  given  equation  of  the  form  (3)  is  stable  and  if  so  , 
to  integrate  this  equation  . 

To  solve  this  problem  ,  one  may  write  all  the  equations 
(3)  obtained  from  the  canonical  equations  by  a  transforma¬ 
tion  (4)  .  However  ,  this  new  Table  of  equations  is  far  too 
extended  to  be  written  in  full  .  According  to  Gambler  ,  certain 
sub-classes  give  rise  to  more  than  one  hundred  equations  . 

Fortunately  ,  our  method  associates  with  each  given 
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equation  ,  a  small  number  of  integers  ,  easily  determined 
and  giving  immediately  the  corresponding  sub-class  ;  to  perform 
the  integration  ,  one  has  only  to  refer  to  the  appropriate 
paragraph  of  the  paper  ,  The  situation  is  thus  similar  to  the 
theorem  of  Fuchs  in  case  of  regular  integrals  of  a  linear 
differential  equation  . 

Finally  ,  v/e  observe  that  our  method  applies  also  to 
equations  of  order  higher  than  two  ,  not  necessarily  linear 
in  the  derivative  of  higher  order  .  This  will  be  the  subject 
of  another  paper  , 
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I .  Statement  of  the  problem  . 

1 ,  Let  X  be  a  complex  independent  variable  and  e  a  complex 
parameter  .  Let  y  =  {yj^,...,y^)  be  a  set  of  functions  satis¬ 
fying  the  system  of  ordinary  differential  equations 

(1.1)  =  A  ( X ;  y ;  e  ) 

dx 

where  A  =  (A^,.,,,A^)  is  a  given  vector-function  of  x  ,  y  ,  e, 

holomorphic  within  a  certain  domain  D  containing  the  point 

0/0  o  ^  _ 

X  =  ,  y  =  y  =  (y^  ,  ,  y^)  ,  e  =  0  . 

When  e  =  0  ,  the  system 

(1.2)  =  A(x,y;0) 

dx 

is  called  the  reduced  (  or  non  perturbed  ,  or  undisturbed) 
system  corresponding  to  the  perturbed  system  (l.l)  . 

It  is  known  that  the  system  of  differential  equations 
(1.1)  admits  a  unique  solution  y(x)  =  y(x;e)  =  y(x-x^;y^;e) 
which  is  holomorphic  within  a  certain  domain  containing  the 
point  X  ,y  ,e  =  0  and  which  reduces  to  y  when  x  =  x  , 

On  using  a  theorem  of  Poincare  ,  basic  in  perturbation 
theory  ,  one  may  write 

oo 

(1.3)  y(x)  =  y(x;e)  =  v(x)  +  2  ^  (x) 

p=l  P 

where  v(x)  =  (v^,..,,v^)  is  a  solution  of  the  reduced  equation 
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(1.4) 


dv 


=  A(x;v;0)  . 


dx 


This  solution  v(x)  depends  on  n  arbitrary  constants 
(for  instance  y°  and  is  holomorphic  for  certain 

values  of  these  constants  Cj^  when  x  varies  along  a  cur^'^e  L 
joining  two  points  of  the  domain  D  . 

On  setting  the  matrix 

(1.5) 


r./  \  8A.  (x;v;0) 

B(x;v)  =  , 


3v, 


one  may  write 


(1.61 


A(x;v+  E  e^u  ;e)  =  A(x;v;0)+e  T  B(x;v)u-+C,  (x;v)'l 
p=l  P  L  1  1  J 

|^B(x;v)  u^  +  Cj^(x;v,u^,  . . .  ,Uj^_^)J  , 


00 

^=2  kl 


where  Cj^(x; v,u^ , . , .  ,u^  is  a  vector  with  components  depen¬ 
ding  only  on  the  components  of  the  vectors  v  ,  u^,  .  . ,  , 

Therefore  ,  Uj^  is  a  solution  of 

(1.7)  '^^k  =  B(x;v)  Uj^  +  Cj^(x;v,u^,  ...,Uj^_^)  , 


( k  1^2f.«*)  • 

The  solution  of  this  non  homogeneous  system  of  linear 
differential  equations  may  be  obtained  by  the  method  of 
variation  of  parameters  when  the  solution  of  the  homogeneous 
system 
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(i.e) 

is  known  . 


du 

dx 


B  ( X ;  V )  u 


Moreover  ,  if  v(x;c) 
solution  of  (1.4)  »  depending  on  the 


c^)  ,  is  the  general 

Chi.  jr  w  k;}  w  (.xix  V  o 


c  ^  I  . . . 

(1.9) 


c  .  one  has 

n  ’ 

d  0v 
dx  0Cj^ 


(k=l, . . . ,n) 


so  that  ________  ,  , , .  ,  ______  are  n  linearly  independent  solu- 

0Ct  0c 

1  n 

tions  of  (1,8)  . 


2,  Classification  of  singularities  ,  The  solution  y(x) 

=  y(x-x^,y^;e)  of  (1.1)  together  with  the  series  obtained 

by  analytical  continuation  defines  a  function  -  call  it  again 

y(x)  -  which  is  a  solution  of  (1,1)  and  in  which  x^  ,  y^ 

appear  as  parameters  .  This  function  y(x)  may  have  singular 

points  depending  or  not  on  the  initial  value  y^  of  y  .  Those 

singular  points  of  y{x)  independent  of  y  are  called  fixed 

o 

or  intrinsic  singular  points  of  y(x)  or  of  the  differential 
equations  (1.1)  ;  those  singular  points  of  y(x)  which  depend 
on  y^  are  known  as  the  movable  or  parametric  singular  points 
of  y(x)  or  of  the  differential  equations  (1,1). 

When  the  integral  y(x)  has  no  parametric  critical  points 
(i.e.  branch  points  or  essential  singularities  )  ,  then  we 
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say  that  y(x)  is  a  stable  integral  of  (1.1)  ;  when  the  diffe¬ 
rential  system  (l.l)  has  only  stable  integrals  ,  we  say  that 
this  system  is  stable  .  In  the  other  cases  ,  the  differential 
system  and  the  corresponding  integrals  are  unstable  . 

The  following  problem  then  arises  :  to  find  necessary 
and  sufficient  conditions  in  order  that  a  given  system  (1.1) 
be  stable  . 

To  solve  this  problem  ,  the  following  theorem  is  essential 

General  theorem  of  stability  ;  If  the  general  solution  of 
the  differential  system  (1.1)  is  uniform  in  x  for  all  values 
of  e  ^  D  except  possibly  e  =  0  ,  then  it  will  also  be  uniform 
for  e  =  0  ,  Moreover  ,  the  coefficients  Up(x) , (p=l,2, . , . ) ,  of 
the  series  (1,3)  are  also  uniform  . 

To  prove  this  theorem  ,  consider  in  the  x-plane  ,  a  closed 
path  L  beginning  and  ending  at  x^  ,  on  which  y(x,e)  is  ana¬ 
lytic  .  Let  u,  (x)  be  the  first  of  the  functions  v(x)  ,u  (x)  , 

K  p 

(p=*l,2,...)  ,  which  takes  on  two  (or  more  )  values  at  x^  ; 
then  ,  one  has 

y(x,e)  =  v(x)+euj+, ,  ,+e^”^Up_^+e^  t  ’^p'*’^'^p+l'^*  *  • "] 

=  Yj(x,e)  +  Yg(x,e)  . 

When  e  is  small  enough  ,  the  values  of  Yg(x;E)  are  very  near 
to  the  values  of  Up(x)  ;  therefore  ,  if  Up(x)  is  multi-valued 
at  x^  ,  then  Yg(x,e)  is  also  multi-valued  at  x^  .  Furthermore, 
the  number  of  values  of  y(x,e)  at  x^  is  not  less  than  the 
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numbers  of  values  of  any  of  the  functions  v(x)  ,Uj^(x)  ,  (k=l,2,  . .  , )  , 
at  x^  .  Consequently  ,  if  y(x,e)  is  uniform  ,  then  the  functions 
v(x)  ,  ,(k=l,2,...)  are  also  uniform  . 

The  method  initiated  by  P.  Painleve  to  solve  the  problem 
breaks  up  into  two  parts  .  First  ,  a  set  of  necessary  condi¬ 
tions  for  the  absence  of  parametric  critical  points  is  obtai¬ 
ned  and  a  set  of  equations  satisfying  these  necessary  conditions 
is  derived  ;  second  ,  it  is  shown  by  direct  integration  or  by 
using  an  appropriate  method  that  these  particular  equations 
are  stable  . 

To  obtain  the  necessary  conditions  ,  a  parameter  e  is 
introduced  into  the  differential  system  under  consideration  , 
in  such  a  way  that  the  new  system  has  only  fixed  critical 
points  and  is  integrable  when  e  =  0  ,  The  functions  v(x),Uj^(x) 
are  solutions  of  systems  of  linear  differential  equations  and 
thus  determined  by  quadratures  .  The  conditions  that  -  their 
parametric  branch  points  are  fixed  give  rise  to  necessary 
conditions  for  the  stability  of  the  given  system  . 

The  necessary  conditions  thus  obtained  enable  one  to 
simplify  the  given  system  ;  the  same  method  is  again  applied 
until  no  further  necessary  conditions  of  stability  are  obtained. 

This  method  ,  theoretically  simple  ,  requires  heavy  and 
tedious  calculations  . 

To  avoid  these  combersome  calculations  ,  another  method 
is  developed  in  what  follows  ;  it  systematically  avoids  ,  as 
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far  as  possible  ,  any  integration  of  systems  of  linear  diffe 
rential  equations  and  reduces  the  problem  to  elementary  alge 
braic  processes  . 

II .  The  theorems  of  stability  , 

3 .  In  the  next  paragraphs  ,  we  shall  apply  the  general 
theorem  of  stability  to  various  important  particular  cases  ; 
the  results  thus  obtained  will  prove  very  useful  in  what  fol 
lows  . 

Throughout  this  paper  ,  we  denote  by  dobs  differentia¬ 
tions  with  respect  to  the  independent  variable  x  . 

Let  k  0  be  an  integer  and  A  ^  0  a  constant  , 

IfelESl®  ^  •  In  order  that  the  equation 

(3.1)  y^^  y  =  A 

be  stable  ,  it  is  necessary  and  sufficient  that  k  =  0  . 
Indeed  ,  the  general  solution  of  (3.1)  is  given  by 

(3.2)  y^^^  =  +  A(x  -  x^) 

where  x^  and  y^  are  arbitrary  constants  ;  it  is  clear  that 
y(x)  defined  by  (3.2)  is  stable  if  and  only  if  k=0  . 

Theorem  II,  In  order  that  the  equation 

(3.3)  y  =  A  y^ 

be  stable  ,  it  is  necessary  and  sufficient  that  k  ^  2, 

Indeed  ,  the  general  solution  of  (3,2)  is  given  by 
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(3.4)  =  (i-k)  [‘a(x  -  x^)  +  , 

where  x^  and  c  are  arbitrary  constants  .  In  order  that  y(x) 

defined  by  (3.4)  be  stable  ,  it  is  necessary  and  sufficient 

that  1  -  k  =  ^  where  n  is  an  integer  ;  then  k  ^  2  . 

n 

Theorem  III,  Let  P(y)  be  a  holomorphic  function  of  y  in  a 
neighborhood  of  y  =  0  ,  In  order  that  the  system  of  differen¬ 
tial  equations 


(3.5; 


y  =  z  » 

y^  z  =  A  [1  +  P(y)^ 

be  stable  ,  it  is  necessary 

i.  that  k  =  0  or  1  ; 

ii ,  and  if  k  =  1  ,  that  A  =  1  -  ^  ,  where  n  ^  0  is  a 

n 

positive  or  negative  integer  or  n  =  oo(i.e.  A  =  1)  . 

Indeed  ,  set 

■  k 

y  =  e  u  =  e  Uj  +  t  u^  +  . , .  , 

k  k  k+l  , 

z=e  v=e  Vj+e  v^f...  ; 

the  differential  system  (3.5)  becomes 

;  =  V  , 

u^  V  =  A  v^  1  +  P(eu)^  . 

It  then  follows  that  u^^  ,Ug, .  satisfy  the  following  equa¬ 
tions 


therefore  ,  if  a  ,  b  ,  c  are  arbitrary  convstants  with  the 
restriction  that  b  ,  c  0  ,  one  has 
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u^(x)  =  b  , 


k 

b  c 


-  Ac(x-a) 


u„(x)  =  I  V  (t)  dt  =  _  Ig 

’  a  Ac 


b*^  -  Ac^x  -  a) 


In  order  that  u^lx)  be  stable  ,  it  is  necessary  that  k  =  0 


or  1 


Now  ,  suppose  k  =  1  and  set  x=a+Et  ,y=eu, 
z  =  z  ,  where  a  is  a  constant  .  Substitution  in  (3.5)  gives 


du 

~ 


=  z 


u 


dz 


=  A  z 


dt 


^  [l  +  P(eu)]  . 


(3.6) 


When  e  =  0  ,  one  finds  the  reduced  system 
du 


dz  A  2 
u  =  A  z 


dt  dt 

The  function  u(t)  defined  by  (3.6)  satisfies  the  equa¬ 


tion 


^  A 

u  =  A 


dt' 


(•^) 


and  implies  : 

i.  if  A  =  1  , 

ii,  if  A  ^  1 


ct  +  Ct 
u  =  e  1 


u  =  (c^  +  ct ) 


1/(1-A) 


c  and  c^  are  arbitrary  constants  .  Therefore  ,  in  order  that 
u(t)  be  stable  ,  it  is  necessary  that  A  =  1  or  ^  n  , 


1  -  A 


where  n  ^  0  is  a  positive  or  negative  integer  note  that 
n  =  CO  gives  A  =  1  J  . 
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5.  In  this  and  In  the  next  paragraphs  ,  let  p  denote  a 
constant  ,  k  0  an  Integer  ,  P{x;z)  a  polyxiomial  in  z  of 
degree^  k-2  ,  and  H^(x;z;u)  ,(i=l,2,...)  ,  a  polynomial  in 
u  ,  Further  ,  we  suppose  that  h(x;z)  ,  H^(x,z;u)  are  holomor- 
phic  functions  of  z  at  z  =  0  and  that  P(x;z)  ,  h(x;z)  ,H^(x;z;u) 
are  analytic  functions  of  x  in  a  given  domain  D  , 

We  consider  the  differential  system 

z  =  1  +  z  P(x;z)  +  z  u  , 

. 

z  u=  pu+h(x;z)  +  z  H^(x;z;u)  , 

concerning  which  we  have  the  following  basic  theorem  : 

Theorem  IV  .  In  order  that  the  differential  system  (5.1) 
be  stable  .  it  is  necessary 

i,  that  p  be  an  Integer  ,  positive  ,  negative  ,  or  zero  , 
li.  and  if  p  =  0  ,  that  h(x;0)  =  0  . 

Proof  .  Set  x=a+et  ,  z=ev  ;  a^Disa  constant  . 
Substitution  in  (5«1)  yields 

j  k 

=  1  +  e  V  P(a+E  ,  t  ;  Ev  )  +  e  v  u  , 


=  p  u  +  h(a+et;ev)  +  ev  H^(a+et,ev;u)  . 


When  E  =  0  ,  one  finds  the  reduced  system 


(5.2) 


=  p  u  +  h(a;0)  . 


Suppose  p  0  .  Then  the  general  solution  of  (5.2)  is 
v(t)  =  t  -  b  ,  u(t)  =  c(t-b)P  -  ; 


17 


b  and  c  are  arbitrary  constants  . 

In  order  that  u{t)  be  stable  ,  it  is  necessary  that  p 
be  an  integer  ,  positive  ,  negative  ,  or  zero  , 

Now  suppose  p  =  0  ;  the  non  perturbed  system  i5»2)  becomes 


dt  dt 


whose  general  solution  is 

v(t)  =  t  -  b  ,  u(t)  =  c  +  h(a;0)  Ig(t-b)  , 

where  b  and  c  are  again  arbitrary  constants  ,  Therefore  ,  in 
order  that  u(t)  be  stable  ,  it  is  necessary  that  h(a;0)  =  0  ; 
because  a  D  is  arbitrary  ,  this  condition  becomes  h{x;0)  =  0. 


6.  In  what  follows  ,  we  shall  apply  our  basic  theorem  IV 
to  the  differential  system 

z  =  1  +  z  u 

(6.1) 


z  u  =  p  u  +  h(x;z)  +  z  H^(x;z;u)  , 
where  p  is  an  integer  and  h(x;z)  ,  H^(x,z;u)  have  the  properties 
indicated  above  |^see  ^5^  • 

When  p  is  a  negative  integer  ,  no  condition  for  stability 
follows  from  theorem  IV  , 

Suppose  p  ^  0  and  set 
(6.2)  u  =  P(x;z)  +  z^  V  , 

where  v  is  an  unknown  function  and  P(x;z)  is  a  polynomial  in 
z  of  degree  p-1  whose  coefficients  will  be  determined  later  on. 
Substituting  in  (17)  ,  one  finds 


le. 


z  =  l4-zP(x;z)+  z^"^^  V 
so  that  ,  from  (6,2)  ,  it  follows  that 

u  =  Q{x;z)  +  z^  V  +  p  z^"  V  +  z^  Hg(x,z;v)  ; 
Q(x;z)  is  a  polynomial  in  z  of  degree^  p»l  ;  H2'(x;z;v)  has 
the  same  properties  as  H^(x,z;u)  , 

Now  ,  one  may  write 

H^(x,z;u)  =  Hj^(x,z  ;P+z^v) 

=  R(x,z)  +  z^  H^lxjzjv)  ; 

R(x,z)  is  a  polynomial  in  z  of  degree^  p-1  and  is  an  ana- i 

lytic  function  of  x  ;  H^(x,z;v)  has  the  same  properties  as 

H^(x,z;u)  .  Then  ,  the  second  equation  (6.1)  becomes 

z^"*"^  V  =  E(x;z)  +  H^(x,z;v)  , 

where  H.  ==  H_  -  H„  and 
4  3 

E(x;z)  =  p  P(x;z)  -  z  Q(x;z)  +  h(x;z)  +  z  R(x;z) 
is  a  polynomial  in  z  of  degree  p  , 

Suppose  that  we  determine  the  p  coefficients  of  the 
polynomial  P(x;z)  in  order  that  E(x;z)  5  ^(x)  z^  ;  then  , 
differential  system  (6.1)  becomes 

\  z  =  1  +  z  P(x;z)  +  z^"^^  V  , 


(6.3) 


V  =  (^(x)  +  Z  H^(x,z  ;v)  . 


Then  the  condition  for  stability  ,  £  (*)  -  0  *  follows 
from  theorem  IV  . 

It  is  clear  that  when  £  (x)  =.  0  ,  the  system  (6.3)  has 
a  unique  holomorphic  solution  z(x)  ,  v(x)  such  that  z(x  )=0 


v(x  )  =v- ,  ,  where  v  is  an  arbitrary  constant  . 

o  p.  o 


the 


» 
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Note  that  the  condition  £(x)  =  0  is  also  sufficient  in 
order  that  the  general  solution  of  the  differential  system 
(6.3)  have  no  branch  point  in  .  Therefore  ,  no  additional 
condition  for  stability  can  be  obtained  by  the  preceding  method; 
its  efficacy  is  exhausted  . 


7.  It  is  particularly  easy  to  dispose  of  the  above  method 
when  p  is  small  .  For  future  use  ,  it  will  be  useful  to 
consider  the  values  p=0,  1,2,3» 

Suppose  that  the  second  equation  of  (6,1)  is  written  more 
explicitly  in  the  form 


(7.1; 


z  u  =  p  u  +  A  +  A,  z  +  A„  z^  +  A„  z-^  + 
o  1  2  3 

+  u(B^  z  +  +  ...  )  +  k  z  u^  , 


where  the  A^’s  and  B^’s  are  analytic  functions  of  x  and  k 


0  . 


is  a  number  Independent  of  x  . 

i,  p  =  0,  The  condition  for  stability  is  A^ 
il.  p  =  1  .  Set  u  =  a  +  z  V  ,  where  a  is  to  be  determined 
later  on  .  Substituting  in  (7»1)  »  one  obtains  the  two  relations 


(7.2) 


“  +  *0  *  ° 


a  •  =  A^  +  a  B^  +  k  « 


On  eliminating  a  ,  it  follows  that  the  condition  for 
stability  is 


(7.3) 


A  +  At  -  A  B-  +  k  A'"  =  0  . 
o  1  o  1  o 


iii,  p  =  2  .  Set  u=o+p  z  +  z^v,  where  o  and  p  are  to 
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be  determined  later  on  ,  Substituting  in  (7.1)  »  one  obtains 
the  three  relations 


{  ^  2  a  +  =  0  , 

(7.4)  1  a  =  P  +  +  a  +  k  a  ^  , 

I  p  =  A^  +  a  +  P  +  (2k-l)  ap 

To  obtain  the  condition  for  stability  ,  one  must  eliminate 
a  and  p  from  this  set  of  equations  . 

2  3 

iv.  p  =  3  ,  Set  u=a+pz+Yz  +zv;  the  same 
method  gives  rise  to  the  four  relations 

3  a  +  A  =0  , 

o  * 

*  2 

a  =  2p+A^  +  otB^  +  kot  , 

(3  =  Y  +  Ag  +  a  Bg  +  p  B^  +(2k-l)  ap  , 

Y  =  A^  +  «  B^+  P  Bg  +  Y  B^+(k-l) ( P  ^+2  CY  ) , 


(7.4) 


from  which  the  condition  for  stability  follows  by  the  elimi¬ 
nation  of  a  ,  p  ,  Y  . 

In  general  ,  the  process  of  elimination  is  complicated  . 
Fortunately  when  applied  to  particular  equations  ,  the  condi¬ 
tion  for  stability  is  ,  in  most  cases  ,  readily  obtained  and 
easy  to  work  with  .  Examples  will  be  given  later  on  , 


III,  Applications  . 


S.  To  find  all  the  st.able  equations  of  the  form 


(d.i) 


*2 

y 
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where  P(y)  p  Q(y)  are  polynomials  in  y  independent  of  x  , 

1.  Q(y)  Is  a  constant  or  has  only  simple  roots  . 

Let  y  =  a  be  a  root  of  Q{y)  of  multiplicity  1  ,  The 

equation  (8,1)  is  equivalent  to  the  differential  system 


(6.2) 


y  =  z 


A  z' 


(y  -  a) 


1  +  gly)"] 


where  A  is  a  constant  and  g(y)  is  a  holomorphic  function  of 

y  in  a  neighborhood  of  y  =  a  and  such  that  g(a)  =  0  , 

According  to  theorem  III  ,  a  necessary  condition  for  this 

system  (26)  to  be  stable  is  k  =  1  and  A  =  1  or  1  -  ^  , where 

n 

n  5^  0  is  a  positive  or  negative  integer  . 

ii.  Let  a^,,,.,ap  be  the  p  simple  roots  of  Q(y)  ,  Denote 
by  A,  the  residue  of  at  y  =  a,  ;  one  has  A,  =  1  or 

"WT  ^ 

Aj^  =  1  -  ____  I  where  n^^  ^  0  is  a  positive  or  negative  integer 
(cf.i).  Then  (8.1)  may  be  rewritten  as 


(8.3)  y  =  y^r  S  \  ■  +  R(y)  1  , 

''  k=l  y  -  aj^ 


where  R(y)  is  a  polynomial  in  y  . 


Now  ,  set  y 

=  u-1 

;  equation 

(8.3) 

becomes 

P 

s 

^k 

+ 

f— 

L 

k=l  1 

u  ' 

1  u  / 

u^  B(u)  . 


According  to  theorem  III 


1 


/  1 


u  =  0  must  be  a  simple  pole 


I 


of  B(u)  so  that 


u 


R 


=  b  +  0(u)  ,  where  b  is  a  constant 
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Therefore  ,  the  polynoinial  R{y) 


^  ^  ^  ^  is  identically 

y  V  y  ; 


zero  . 


Moreover  ,  It  follows  from  theorem  III  that 


A  =  S  A, 


k=l 


is  equal  to  1  or  to  1  -  ^  ,  where  n  ^  0  is  a  positive  or 

n 

negative  integer  , 


9.  It  follows  from  the  preceding  paragraph  ,  that  to  obtain 
all  the  stable  equations  of  the  form  (6,1)  ,  one  has  to  find 
the  set  of  integers  u  ,  ,  (k=l,,,,,p)  different  from  zero 

and  satisfying  the  equation 

(9.1)  S  f  1  -  ^  \  =  1  -  ^  . 

k=l  V  ”-k  ^  ^ 

Note  that  n  and  may  be  infinite  . 

1.  Suppose  A  -  1  or  n  =  oo  ;  then 


therefore  p  4  2  and 

a.  p  *  1  f  ^ 

b.  p  =  2  ,  n^  =  n^  =  2  , 

ii.  Suppose  A  >>  1  ;  then  A  =*  1  -  ^  4.  2  .  Because 


n 
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1  - 


^  ^  ^  f  one  has  p  4;^  4  • 


Hk  //  2 


a.  p  =  4  ;  then  A  =  2  ,  A,.  =  ^  ,  (k*l,2,3»4)  t  i»e» 

^  "2" 

n  =  -1  ,  n^  =  n^  =  n^  =  =  2  , 

b.  p  =  0  ;  then  Ak  =  0  ,  (k=l,2,3»4)  • 

c.  p  =  1  ;  then  A  =  A.  =  1  -  ^  ,  n  is  an  integer  ^0,1, 

n 

d.  p  =  2  ;  then  ^  ^  ^  -  1  0  so  that  n  is  negative 

n  n^  n^ 

On  setting  n=-m,m^0,  one  has 


(9.2) 


1  t  1 


m 


=  1 


n. 


1  “2 

whose  solutions  are 

(3»3»3)  >  (2^i+>4)  »  (2»3»6) 


(2,2,00  ) 


( l,n,-n 


where  n  is  an  integer  , 


e,  p  =  3  ;  then 


1  +  1  + 


n 


n. 


n^ 


^  -2<f^-2/0 


1  2 

and  n  is  negative  ,  On  setting  n  =  -m,ni^0,  one  has 

1 
n 


(9.3) 


1  =2-  1 


m 


1  1 

"nj"  ^  T 


“1  2 

j^note  that  Uk  ^  1  ^  ;  therefore  m  =  1  or  2  . 

If  m  =  1  ,  equation  (9.3)  reduces  to  (9.2)  with  n  =  n^. 
If  m  =  2  ,  one  has  =  n^  =  n^  =  2  .  Note  that  n^  >  2 


is  impossible  because  in  that  case 


n. 


implies 


1 

"nl 


+  ^  ^  1  ,  an  impossibility  since  n^  ,  n^  are  different 

n„  ^ 
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from  1  .  Thus  the  solution  is  unique  . 

Corresponding  to  the  solutions  of  equation  (9.1)  t 
are  nine  possible  types  of  equations  of  the  form  (6.1)  . 
convenience  ,  we  set 

A  =  ^  ,  C  =  _  ,  E  *  ^  , 

ay  +  b  cy  +  d  ey  +  f 

where  a,b,c,d,e,f ,g,h  are  constants  and  eventually  zero 

The  nine  possible  values  of  are  tabulated  below 

“oTTr 


there 

For 


G=  g 
gy  +  b 


Table  II  . 

0  ;  A  ;  A  +  C  ; 

A(1  +  ^  )  +  C(1  -  ^  )  ,  n  >  1  an  integer  ; 
n  n 

_2_  (A  +  C)  +  E  ;  _f_(A  +  C  +  E)  ; 

2  3 

^A+^(C+E)  ;  ^A+^C+^E  ; 

2  4  3  2 

_J__  (  A  +  C  +  E  +  G  )  . 

2 

The  corresponding  differential  equations  are  shown  to  be 
stable  by  direct  integration  .  Indeed  ,  they  are  equivalent  to 

;-o  , 

y  =  K(ay  +  b)  , 

• 

y  =  K(ay  +  b) (cy  +  d)  , 
y’^  =  K(ay  +  b)’^’*’^(cy  +  d)*^”^ 
y^  =  K(ay  +  b)(cy  +  d)(ey  +  f)^ 


f 
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=  K(ay  +  b)^(cy  +  d)^{ey  + 
y^  =  K(ay  +  b)^(cy  +  d)^(ey  +  f)^ 
y^  =  K(ay  +  b)^(cy  +  d)^(ey  +  f)^ 
y^  =  K(ay  4-  b)(cy  +  d)(ey  +  f)(gy  + 
where  K  is  an  arbitrary  constant  . 

The  solutions  of  these  equations  are  known  and 
only  elementary  or  elliptic  functions 
they  are  all  stable  ;  the  given  list  is  exhaustive  . 

These  results  will  be  used  later  on  . 


h)  , 

involve 
,  Therefore, 


10.  For  future  use  ,  we  need  the  following  theorem  . 
Denote  by  a  and  b  two  constants  ,  not  both  zero  ,  and  by 
k  and  n  two  integers  .  Suppose  k  1  ;  let  n  ^  0  be  positive 
or  negative  and  eventually  infinity  , 

ShggESffi  V  .  In  order  that  the  equation 


be  stable  ,  it  is  necessary  that  k  =  1  and  n  ^  1  . 

Proof  ,  Note  that  equation  (10.1)  is  invariant  under  the 

transformation  y=ez  ,x=et  ,  where  e  is  a  parameter  ; 

note  also  that  under  this  transformation  ,  one  has 

y.1  -  k  dx  _  ^1  -  k  dt 
dy  dz 


Accordingly  ,  we  set 


26 


u  =  y 


1  «  k  dx 
dy 


then  equation  (10,1)  is  equivalent  to  the  system 
k-1  dy  _  1 


dx 


u 


y‘‘ 


where 

P(u)  ~  ^  -  k  -  a  u  -  b 

n 

This  system  is  also  equivalent  to  the  system 


(10.2) 


yk-1  dy  ^  1 

dx  ”Tr 


(10.3)  y  =  u  P(u)  . 

dy 

Suppose  k  1  ;  then  ,  the  equation  (10.3)  has  a  solution 

u  =  h  ,  where  h  ^  0  is  such  that  P(h)  =  0  .  Then  ,  from  (10.2), 
k*  1 

one  has  dx  =  h  y  dy  and 

k  k  ,  .  v 

y  =  _  (x  +  c)  , 

h 


where  c  is  an  arbitrary  constant  ;  this  function  y(x)  is  not 
uniform  . 

Suppose  k  =  n  =  1  ;  then  ,  equation  (10,1)  is  equivalent 
to  the  differential  system 

(10,4)  y  =  ^  ,  y  u  +  a  u  +  b  u^  =  0  . 

u 


Set  x=et  ,u=ev  where  e  is  a  parameter  ;  then  , 
the  differential  system  (10,4)  may  be  rewritten  as 


dy  _  1 

dt  "  V 


dv 


+  e  v(a  +  bev)  =  0  . 
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ViTe  may  now  apply  the  general  theorem  of  stability  .  To 
do  this  ,  y(t)  ,  v(t)  are  developed  as  series  of  ascending 
powers  of  the  parameter  e  ,  If  ,  y^  ,  are  arbitrary 
constants  ,  one  finds 

y  =  Yq  +  ^  ~  +  0(e)  , 


V  =  V  -  e 


a  Ig 


+  O(e^) 


if  a  ^  0 


V  =  -  e^b  Ig  ^  y^  +  ^  ~  j  +  O(e^)  if  a  *  0  . 

o 

It  then  follows  readily  that  the  differential  system  (10. if) 
and  the  equation  (10.1)  are  not  stable  .  Our  theorem  is 
proved  , 


IV,  The  equation  y  =  R(x,y,y)  . 

11.  Consider  the  equation 

(11.1)  y  =  R(x,y,y):  , 

where  R(x,y,y)  is  a  rational  and  irreducible  function  of  y 
and  y  ,  with  coefficients  analytic  in  x  ,  Our  purpose  is  to 
find  all  the  stable  equations  of  the  type  (11,1)  . 

The  equation  (11,1)  is  equivalent  to  the  differential 
syst em 

•  • 

(11.2)  y  =  z  ,  z  =  R(x,y,z)  . 

Let  x^  ,  y^  be  particular  (  arbitrary  )  values  of  x  ,  y  ; 


the 


let  z  ==  a(x  ,y  )  be  a  pole  of  order  m  >  0  of  R(x  ,y  ,z)  . 

oo  ^  o  o’ 

Set  z  —  u  +  a(x,y)  ;  then  ,  in  a  neighborhood  of  > 

differential  system  (11,2)  may  be  rewritten  as 

y  =  u  +  a{x,y)  , 

(11.3)  I  . 

u  =  u  P(x,y,u)  , 

where  P(x,y,u)  is  a  holomorphic  function  of  x,y,u  such  that 
P(x^,y^,0)  ^  0  . 

Let  e  be  a  parameter  .  On  setting 


(11.4)  X  =  x^  +  t  , 


y  =  y  +  w  ,  u  *  v  , 


the  differential  system  (11.3)  becomes 

/  £<1“P  V  +  a(x  +c^t,y  +eS/)  , 


(11.5) 


dt 


V 


r  ( m+ 1 )  — p  m  d v  _  <  ,  p.  ,  a  r  > 

^  V  _  =  P(x  +e^t,y  +e^,e  v) 

dt  °  ° 


so  that  ,  assuming  p  —  q  »  p  =  r(m+l)  ,  r  =  1  ,  one  obtains 

a(x  +  eP  t  ,  y^+ePw  )  +  e  v  , 


(11.6) 


dw 

dt 


®  'Iv  _/  ,  p  .  .  P  , 

_  =  P(x  +  t,y  +  w  ,  e  V  )  . 


dt 


When  E  =  0  ,  one  finds  the  reduced  system 


dw  ,  ,  m  dv  _ , 


for  this  system  to  be  stable  ,  it  is  necessary  ,  according  to 

theorem  I  ,  that  m  =  0  . 

♦  • 

Therefore  ,  R(x,y,y)  is  a  polynomial  in  y  of  degree  s  ; 


let  us  write 
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R(x,y,y)  «  R^(x,y)  +  R^(x,y)  y  +  ...  +  Rg(x,y)  y®  . 

To  determine  the  value  of  the  positive  integer  s  ,  observe 
that  the  equation  (11.1)  is  now  equivalent  to  the  differential 
system 


(11.7) 


r 

t 


y  =  z  , 

*  g 

z  =  RQ(x,y)  +  ...  +  Rg(x,y)  z  . 


On  setting 

(11.6)  X  =  x^  +  t  »  y  =  y^  +  u.  ,  z  =  e”^  v 
one  obtains 


(11.9) 


^q-p+r  du  ^  ^ 

“dt" 


g(s-l)r-p  dv  =  R  (x  +e^t,y  +e%.)v®  +  0(e). 
dt 


s  o 


When  q-p+r  =0  ,  (s-l)r-p=0  ,  r*  1  ,  i.e.  ,  p*s-  1, 


q  =  s-2  ,  one  finds  the  reduced  system 


au 

dt 


=  V 


dv 

dt 


-  v' 


for  this  system  to  be  stable  ,  it  is  necessary  ,  according  to 
theorem  II  ,  that  8^2. 

Therefore  ,  in  order  that  equation  (ll.l)  be  stable  ,  ^ 

•  • 

is  necessary  that  R(x,y,y)  ^  ,  with  respect  to  y  ,  a  polynomial 
of  the  second  degree  . 

For  convenience  ,  we  shall  write  equation  (11.1)  as 

(11.10)  y  =  A(x,y)  y^  +  B(x,y)  y  +  C(x,y)  , 

where  A,B,C  are  rational  functions  of  y  with  coefficients  analy- 


30 


tic  in  X  . 


12.  The  next  step  in  our  investigation  is  to  characterize 
the  functions  A,B,C  regarded  as  functions  of  y  .  Precisely  , 
we  shall  show  that 

1.  All  the  poles  of  A{x,y)  regarded  as  a  function  of  y  are 
simple  ; 

ii.  The  poles  of  B(x,y)  and  C(x,y)  regarded  as  functions  of 
y  are  included  among  those  of  A(x.y)  and  are  simple  ; 

iii.  The  degrees  of  the  polynomials  ^  and 

A(x,yy  A(x,y ) 

in  y  are  at  most  4  and  6  ,  respectively  . 

Proof  of  1  .  Set  x  =  x^  +  e  t  ,  where  x^  is  arbitrary  ; 
equation  (11,10)  becomes 

y  -  A(x„  +  e  t  ,  y)  /  ^  +  0(e)  . 


In  order  that  the  reduced  equation 


y 


A(x^»y) 


2 


be  stable  ,  it  is  necessary  that  A{x^,y)  be  identical  to  one  of 
the  nine  types  enumerated  above  J^see  ^  9  »  Table  II  .  Because 
x^  is  arbitrary  ,  the  a  ,  b  ,  c  ,  d  ,  e  ,  f  ,  g  of  these  nine 
types  are  now  to  be  regarded  as  analytic  functions  of  x  ;  these 
functions  may  eventually  be  identically  zero  or  constants  , 


Now  ,  set 
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u 


u  = 


u 


_  c  y  +  d 


af-be  cy  +  d 

u  =  _________  .  _ 

cf-de  ay  +  b 

according  as  A(x,y)  coincides  respectively  with  the  first  type 
of  ^  9  »  Table  II  ,  or  the  second  ,  or  the  third  ,  or  the 
fourth  ,  or  with  one  of  the  other  types  ;  then  ,  one  finds 
that  A(x,y)  may  be  assumed  to  be  one  of  the  following  eight 
distinct  rational  functions 


Table 

III 

i. 

0 

f 

ii. 

1 

9 

y 

ill. 

(1  - 

J_) 

n  y 

9 

n  ^  1  an  integer  ; 

iv. 

1 

+  1 

9 

2  y  y  - 

“iT 

V. 

2 

f  — +  • 

1 

.)  : 

\  ~ 

y  - 

1/ 

Vi. 

3 

1 

= 

"4" 

\  7  7 

-  1 

/ 

vii. 

2 

.  +. 

1 

• 

1 

~ 

7 

2 

y  -  1 

viii , 

1 

-  l—^ 

1 

+  '  'I  . 

“2" 

\  ^ 

y  - 

1  7  -  a  j 

In  viii 

,  H  may 

be  a  constant 

(  0  ,  1  )  or  may  depend 

on  X  ;  in  the  latter  case  ,  on  taking  a  new  variable  t  =  H(x) 


* 
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one  may  suppose 
I^As  an  example 


that  H(x)  coincides  with  x  . 

,  we  verify  the  result  for  u  = 


1 

a  y  +  b 


t 


where  a  and  b  are  analytic  functions  of  x  ;  then  , 


y 


y 


M(x,u) 


+  N(x,u) 


u  +  P{x,u) 


t 


where  M  ,  N  ,  P  are  rational  functions  of  u  with  coefficients 

analytic  in  x  , 

Therefore  ,  the  equation 

*  2  • 
y  =  ^  +  B(x,y)  y  +  C(x,y) 

ay  +  b 

assumes  the  form 

u  =  ________  +  B^(x,u)  u  +  C^(x,u)  , 

where  are  rational  functions  of  u  with  coefficients 

analytic  in  x  . 

Note  also  that  the  third  type  of ^  9  »  Table  II  is  now 
considered  as  a  special  case  of  the  fourth  type  (  for  n  =  oo 

Proof  of  ii  ,  Let  y  =  a(x)  be  a  pole  of  order  j  of  A(x,y) 
of  order  k  of  B(x,y)  ,  and  of  order  of  C(x,y)  ,  One  has 
j  =  0  or  1  ;  we  suppose  that  at  least  one  of  the  integers  k  and 
P  is  greater  than  j  ,  On  replacing  y  +  a(x)  by  y  ,  one  may 
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assume  a(x)  s,  0  , 

In  a  neighborhood  of  the  pole  y  =  0  ,  equation  (11.10) 
may  be  written  as 


(12.1) 


y  =  1  “  JL  +  0(y)}  +  I^B(x.O)  +  0(y)j 

[  C(x,0)  +  0(y)]  , 


y 

1 


where  n  ^  0  is  a  positive  or  negative  integer  which  eventually 
may  be  Infinite  .  Note  that  n  =  1  when  j  =  0  . 


i. 


Now  set  y  =  e  u  and 
k 


X  +  e*”  t 
o 


f  ^  2k  -  1  , 

t  when  6  2k  -  1  . 
Then  ,  equation  (12.1)  becomes 


ii.  X  =  X  +  e 
o 


when 

(1+  P  )/2 


1 

1 

)  - 

• 

+  a  +  0(e)  .  when  P  <^'  Ek-1 

n 

k 

u 

u-  ( 

1 

1 

1  - 

+  ^  +  0(e)  ,  when 

P>  2k-l  , 

v 

n 

u 

a.( 

1 

1 

\ 

• 

+  a  ^ _  +  0(e)  , 

when  P  =  2k' 

n 

i  u 

uP 

here  ,  a 

B(x  , 
o 

0)  ,  b 

=  C(x^  ,0)  . 

The 

reduced 

equations  are  contained  in 

(12.2) 

u 

- 

1  ' 

.p 

\  +  a  “  +  *’ 

V 

”n”; 

'  U  ,k  ,2k-l 

^ if  a  ^  0  ,  k  is  necessarily  a  positive  integer  ;  if  a 
one  has  2k- 1  =  £  ,  where  £  is  a  positive  integer  3  • 


9 


3k 


In  order  that  equation  (12»2)  be  stable  ,  it  is  necessary, 
according  to  theorem  V  ,  that  k  =  1  and  n  1  .  This  proves 
li. 

Proof  of  iii  .  The  result  ii  is  essential  for  it  enables 
one  to  give  upper  bounds  for  the  degrees  of  B  and  C  .  Indeed, 
let  D(x,y)  be  the  least  common  denominator  of  the  partial 
fractions  of  A(x,y)  ,  l.e,  (  see  Table  III)  , 

D(x,y)  =  y  for  ii  ,  iii  , 

=  y(y  -  1)  for  Iv  -  vii  , 

=*  y(y  -  l)(y  -  H)  for  viil  . 

Then  ,  one  may  write 


B(x,y)  - 

D(x,y) 


C(x,y)  - 

D(x,y ) 


where  M(x,y)  and  N(x,y)  are  polynomials  in  y  of  degrees  m  and 
n  ,  respectively  ,  with  coefficients  analytic  in  x  , 

On  making  the  substitution  y  =  u"^  ,  equation  (11,10) 
becomes 

A 


u 


u^ 


u 


[^- 


Now  ,  we  consider  two  cases  according  as  to  whether  A  *  0 
or  A  3^  0  . 

i.  A  =  0  .  Then  ,  u  =  0  is  a  simple  pole  for  the  coefficient 

*2 

of  u  so  that  necessarily  m^l,n^3, 

li,  A  3^  0  .  On  taking  into  account  the  various  types  of 
A(x,y)  given  above  (  Table  III  ,  i-viii)  ,  one  sees  that 
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*2 

u  =  0  is  still  a  simple  pole  of  the  coefficient  of  u  ; 
therefore  ,  if  q  is  the  degree  of  D  in  y  ,  one  has 
+  l  ,  n^q  +  3. 

Thus  ,  q^i4.»ni^24;,n^6  .  More  precisely  ,  if  A(x,y) 
coincides  with  one  of  the  types  i-viii  (Table  III)  ,  the 
degrees  (m,n)  are  (1,3)  for  i  ,  (2,4)  for  ii  and  iii  ,  (3»5) 
for  iv-vii  ,  and  (4,6)  for  viii  . 

13.  The  above  method  enables  one  to  obtain  in  a  simple 
manner  conditions  of  stability  which  restrict  the  form  of 
the  coefficients  A  ,  B  ,  C  of  equations  (11,10).  To  obtain 
other  necessary  conditions  for  stability  ,  one  must  ,  accor¬ 
ding  to  Painleve  ,  introduce  a  parameter  e  into  the  equation 
and  develop  the  solution  y(x,e)  of  the  new  equation  as  a 
series  in  powers  of  e  .  As  we  have  seen  ,  the  coefficients  of 
this  development  are  determined  by  quadratures  ;  the  conditions 
that  their  critical  points  are  fixed  are  necessary  conditions 
for  stability  .  This  method  must  be  applied  to  all  the  poles 

of  A(x,y)  and  to  y  =  00  .  For  y  =  oo  ,one  may  set  y  =  ^ 

u 

or  use  the  transformation  y  =  u/e'3  . 

However  ,  to  avoid  the  numerous  and  arduous  computations 
involved  in  the  method  of  Painlev^  ,  we  ehiftll  use  another  , 
simpler  method  , 
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14.  The  values  of  y  for  which  the  general  existence 
theorem  of  Cauchy  does  not  apply  are  y  =  od  ,  0  ,  1  ,  H{^  0,1), 
or  X  (  according  to  the  type  of  the  equation  )  .  In  the 
following  paragraphs  ,  we  shall  consider  particularly  the 
case  y  =  oo  ;  the  other  cases  may  be  settled  directly  or  by 
appropriate  substitutions  of  the  independent  variable  . 

In  a  neighborhood  of  y  =  oo  ,  one  has 


(14.1) 


/  M{x,y) 

N(x,y) 

"’D(x"iy] 


M 

a  y  t  a.  +  o  , 

^  "d 

by^  +  b,y^  +  b-y  +  b_  +  .  ^o 
1  2  3 


where  are  polynomials  in  y  ,  each  degree  of  which  is 

less  than  the  degree  of  D  ,  and  where  the  a»s  and  b’s  are 
analytic  functions  of  x  , 

Further  , 


(14.2)  A{x,y)  =  +  0(1)  , 

y 

where  the  constant  A  has  the  following  values 


(  0(type  i)  ;  1  (type  ii)  ;  1-  ,n  > 1  an  integer 


n 


(14.3) 


(type  iii)  ; 


o  2i  7 

(types  iv.vi.viii);  ^  (type  v)  ;  (type  vii)  . 

“2“  3  "T" 


These  values  correspond  to  A  =  1  -  ,  n  an  integer  , 

n 


where  n  has  the  following  values 
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( 14»4) 


n  =  l(type  i)  ;  n  =  oo (type  ii)  ;  n  ^ l(type  iii) ; 


n  =  -2(types  iv,vi,viii) ;n  =  -3(typev) ;n=-6(type  vli). 
Then  ,  equation  (11,10)  may  be  rewritten  as 


(14.5)  y  +  0(y)^  ^ _  +|^ay  +  a^^  +  O(y)']  y  +  by^  +bjy^+0(y) 


Let  us  now  determine  necessary  conditions  for  the  absence 
of  parametric  critical  points  for  the  equation  (14.5).  To  do 
this  ,  suppose  that  in  a  neighborhood  of  x  =  ,  one  has 

s(x) 


(14.6) 


y(x) 


(x  -  x^) 


where  r  >  0  and  s(x  )  /  0  ;  s(x)  is  a  holomorphic  function  of 


X  . 


that 


Substitute  y{x)  given  by  (14.6)  into  (14.5)  and  note 


y(x)  =  - 


s(Xo) 


(x-x„) 


r+1 


1  +  0(x  -  * 


y(x)  =  ^^^o^  r  1  +  0(x  -  X  )]  . 

f  \r+2 

(x  -  x^) 

First  ,  suppose  that  a(x^)  or  b(x^)  is  not  zero  ;  then  , 
the  dominant  terms  arise  from  y  ,  y  .»  y  y  »  and  are 

y 

respectively  proportional  to 

I  \-r-2  ,  x-r-2  ,  >-2r-l  ,  \-3J' 

(x-x^)  ,  • 

Therefore  ,  to  obtain  an  identity  at  least  two  of  the  nvunbers 
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r+2  ,  2r+l  ,  3r  must  be  equal  ;  this  gives  r  =  1  . 

Second  ,  suppose  that  a(x^)  and  b(x^)  are  both  zero  ; 

*•  “2*2 

then  the  dominant  terms  arise  from  y  ,  y  »  y  »  y  and  are 

y 

respectively  proportional  to 

(x-x  (x-x  (x-x  (x-x  ) 

o  ''o  '  o  '  o 

so  that  r  =  2  . 


15.  Suppose  r  =  1  and  set 
(15.1)  y  -  »  z  =  1  +  u  z  . 


z 

For  convenience  and  for  future  reference 
following  formulas 


(15.2) 


9 


we  note  the 


(15.5) 

(15.6) 


2r 
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in  addition  , 


(15.7) 

(15.6) 


(15.9) 


A(x,y)  =  A  ^  +  O(z^)  ; 

s 

JL  =  a  +  0(1) 

D  z 


=  a^  +  0(z) 


when  a  ^  0 

when  a  =  0 
when  b  ^  0 


when  b  =  0 


> 


> 


Set 

(15.10)  P  =  b  -  a  s  +  A  +  2  , 

(15.11)  p®3-2A+as  , 

Introduce  y(x)  given  by  (15.1)  into  the  equation  (14.5) 
and  rewrite  the  result  as 


•  p 

(15.12)  z  u  =  -  +  p  u  +  h(x,z)  +  z  H(x,z;u)  , 

z 

where  h(x, z)  and  H(x,z;u)  are  analytic  functions  of  x  and 
holomorphic  functions  of  z  in  a  neighborhood  of  z  =  0  ; 
H(x,z;u)  is  a  polynomial  in  u  . 

Now  ,  we  determine  s  by  setting  P  =  0  ;  it  then  follows 
from  theorem  IV  that  a  necessary  condition  for  the  equation 
(11.10)  to  be  stable  is  that  p  be  an  integer  . 

These  conditions  will  determine  the  possible  values  of 


a  and  b 
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16,  Consider  the  relations 

(16.1)  bs^-as-1-  =  0  , 

n 

(16.2)  p=l+_f_+as  , 

n 

where  p  must  be  an  integer  . 

i.  Suppose  b  =  0  and  a  ^  0  ;  then  a  s  +  1  +  ^  =  0  and 

n 

p  =  ^  .  In  order  that  p  be  an  integer  ,  one  must  have  n  *=  1 

n 

or  n  =  00  ,  so  that 


(16.3) 

n  =  1 

If 

a 

1  , 

as  +  2  =  0  ; 

(16.4) 

n  =  00 

»  P  = 

0  , 

a  s  +  1  »  0  . 

ii . 

Suppose  a 

=  0  and 

b  ^  0 

;  in  order  that  p  ® 

be  an  integer  ,  one  must  have  n=2,  l,-2oroo.  Therefore 
the  following  cases  arise 


(16.5) 

n  =  1 

f 

P  -  3 

f 

b 

2 

s  - 

2 

=  0  ; 

(16.6) 

n  =  2 

$ 

p  =  2. 

f 

2 

b 

- 

3  =  0  ; 

(16.7) 

n  =  -  2 

9 

p  =  0 

9 

2 

b 

- 

1  =  0  ; 

(16.6) 

n  =  oo 

f 

p  =  1 

f 

b 

2 

s  = 

1 

• 

iii. 

Suppose  a 

b  ^ 

0  .  On 

eliminating 

s 

between  (l6.1 

and  (16.2)  ,  one  obtains 


(16.9) 

b(  P  -  1-^) 

2  2  i 

“^  =  1 

■ 

^  n  J 

Taking 

into  account  the 

value  of 

b  given  by  (16,9)  ,  one 

rewrites  (16.1)  as 


2+1 


(p  -  ^)  s"  -  (  p-1-^  )  -  (  l+^)(  P-l-_Lj  '-0. 

This  equation  determines  the  product  as  ;  its  roots  are 

P  -  1  -  -L  ^nd  -  (n+l)(np«n-2) 

n  np  “  1 

so  that  cf.  (16.2)J 

.  (  '  ’ 

1  +  +  as  * 

^  p  +  2  +  n 

_ =  q  • 

\  pn  -  1 

Therefore  ,  the  integers  p  and  q  satisfy  the  Diophantine 
equation 

(16,10)  p+q+n+2=pqn. 


The  integral  solutions  of  (l6,10)  are  given  in  Appendix  I 
and  determine  all  the  distinct  possibilities  that  we  shall 
have  to  consider  later  on  . 

These  possibilities  are  given  below  together  with  the 
corresponding  relation  (16,9)  and  the  related  values  of  p 
and  s  , 


(16,11) 

(16,12) 

(16.13) 

(16.14) 


a  =  0  t 
&  0  , 
9  b  +  a^ 


n  =  1 
b  0 


P  =  3  , 

p  =  1  , 

0  ,  as 

-  3  »  as 

2  ,  as 

5  »  as 


b  s^  -  2  -  0  ; 

a  s  +  2  *  0  ; 

-  3  » 

-  6  ; 

-  1  » 

2  . 


42 


(16.15) 

(16.16) 


a  =  0 


a  ^  0 


n  =  oo 

b  ^  0  ,  P  “  1 
b  =  0  ,  p  =  0 


bs*^  -  1  =  0  ; 

a  s  +  1  =  0  , 


(16.17) 


n  >y  Z  , 

b(n+2)^  +  n  a^  =  0  ,  (  p*0  ,  nas=-(n+2)  , 

1  p»-(n+2) ,nas=-(n+l) (n+2) 


(16.16)  2b  «  a* 


(16.19) 


a  =  0 


n  =  2 
(P-i 
i,  P  -  5 
b  ^  0 


,  as  +  1  * 
,  a  s  -  3 

p  =  2  , 


Sbs*^  -3  =  0. 


(16.20)  3  =  2  b 


n  =  3  . 

p=l  ,  3as  =  -2  , 

p  =  3  ,  3as  =  4. 


(16.21)  5  a*^  =  b 


IL^  . 

^p=l  5as*-2  , 

p  =  2  5as  =  3. 


n  =  -  2  . 


(16.22) 


a  =  0  ,  b?^0  ,  p  =  0  ,  2bs^=l 
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(16.23)  b  =  3 


3  a  s  *  2  , 

3  a  s  +  1  =  0  . 


n  =  ~  6  . 

(16.24)  8  b  =  3  a^  ,  C  p  *  0 

1  P  =  4 


3  a  s  “  -  2 
3  a  s  =  10  . 


17.  Now  suppos-^  r  -  Z  i  then  a(x)  2  0 


b(x)  £  0  .  Set 


(17.1) 


z 


1  +  u  z 


so  that 


(17.2) 

(17.3) 


(17.4) 


Introduce  y(x)  given  by  (17.1)  into  the  equation  (14.5) 
and  rewrite  the  result  as  (15.12)  with 

(17.5)  P  =  b  t  -  2  -  Jl  , 

n 

(17.6)  P  =  1  +  JL  • 

n 


Now  ,  we  determine  t(x)  by  setting  P  =  0  .  In  order  that  the 
equation  (11,10)  be  stable  ,  it  is  necessary  that  p  be  an 
integer  ;  we  suppose  p  ^  0  ,  The  only  five  distinct  possibi¬ 
lities  are 


(17.7) 

n  = 

00 

9 

p  =  1 

9 

t  = 

2 

• 

f 

(17.0) 

n  = 

1 

9 

P  =  5 

9 

t  = 

6 

» 

(17.9) 

n  = 

2 

9 

fO 

II 

a. 

9 

^1 

t  = 

4 

f 

(17.10) 

n  = 

4 

9 

p  =  2 

9 

^1 

t  = 

3 

» 

(17.11) 

n  = 

-4 

9 

o 

II 

p. 

9 

t  = 

1 

When 

n  is 

an 

arbitrary  integer 

9 

one 

has 

also  to  consider 

the  case  asbsb^sO  , 


10,  We  have  also  to  take  into  account  the  values  y  =  0  ,  1 


or  X  , 


For  y  =  0  [_^resp,  y  =  1^  ,  one  may  set  y  =  u*^ 

.  u-I  ] 


resp , 


and  use  the  results  of  the  preceding  paragraphs  . 

One  may  also  proceed  in  a  more  direct  manner  ,  For  instance 
for  y  =  0  ,  one  may  set 
(10.1)  y=sz  ,  z=l+uz 

so  that 


(10.2) 


y  =  s  +  sz  (u+  ^  )  , 


(  y  =  2s  -f  s 


u  +  sz  ^  u 


+  2  u  +  + 

s  s 


and  use  the  above  method  , 

In  what  follows  ,  we  shall  have  to  consider  movable  poles 


$ 


^  or  zeros  or  unities"^  of  order  one  or  two  ;  to  abbreviate 
we  use  the  obvious  notations  |^or  or 

and  for  example  ,  speak  of  an  equation  of  class  (P^  ,  Z^^  ,  U^) 
to  signify  an  equation  having  movable  poles  ,  zeros  and 
unities  ,  all  of  the  first  order  . 

To  continue  our  investigations  ,  we  must  consider  sepa¬ 
rately  each  of  the  eight  types  of  equations  corresponding  to 
the  eight  possible  values  of  A(x,y)  see  Table  III  ,  i-viii"^  . 

Observe  also  that  to  obtain  canonical  forms  for  the  stable 
equations  ,  it  is  often  most  convenient  to  use  a  transforma¬ 
tion  call  it  T(  ^  ,  f  )'\ 

(Id. 3)  y(x)  =  X  (x)  u  +  n  (x)  ,  t  =  <p(x) 

which  does  not  alter  the  main  features  of  the  equations  con¬ 
sidered  L  ^(x)  »  (x)  »  T  (x)  are  analytic  functions  of  x  ; 

in  some  cases  ,  one  has  X  (x)  =  1  ,  r(x)  =  0  or  '?(x)  =  x  , 
For  convenience  ,  and  future  use  ,  we  note  the  follo¬ 


wing  formulas  where  primes  denote  differentiations  with 

du 

"d^ 


respect  to  t  ,  i,e.  u* 


u 


It  = 


d^u 


dt* 


(Id. 4) 


y=Xu+Xcpu'+  |i 


y  =  Xu+{2  X  cp+X  qp)u*  +  X(p^  u”+  ii 


When  (i  =  0 

*2 

,  one  finds 

o  2  V 

*2  V 

••  .  y 

y  -  _1_ 

=  X  f  u"  - 

u’  \ 

J—  + 

(  ^  " 

^  \  u  +  X  v  u»  . 

y 

\ 

^  ) 

V 

^  / 

Throughout 

the  remainder 

of  this 

paper 

,  we  denote  by  a,b,c,d, 
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e,f,g,h,k,  ^  ,  with  or  without  subscripts  ,  analytic  func¬ 
tions  of  X  ,  (not  always  the  same  )  and  by  H  and  K  ,  with  or 
without  subscripts  ,  (arbitrary  )  constants  ,  (not  always 


the  same  ) 

.  Moreover  ,  we  set  A 

II 

h' 

*  =  ^  so  that 

X 

• 

(16.5) 

A  =  ^  . 

qp 

X 
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Agggg^iX  I  • 

The  Blophantlne  equation  p+q+n+2=pqn, 

19.  The  problem  of  finding  all  the  possible  types  of 
stable  equations  (11.1)  depends  on  the  problem  of  finding  all 
the  integl'al  solutions  of  the  Diophantine  equation 

(19.1)  p+q+n+2=pqn. 

We  suppose  n^lorn==-2,-3or-6  and  p  0  , 
i,  n  =  1  .  Equation  (19.1)  is 

(19.2)  p+q+3=pq  . 

When  p  =  0  »  one  has  q  *  -  3  ;  when  p  =  1  ,  equation  (19.2) 
has  no  solution  . 

Suppose  p  ^  2  ,  Equation  (19.1)  may  be  rewritten  as 

(19.3)  -L  (  1  +  _L'|  ^  1  -  — 

q  p  /  p 

so  that  q  >  0  ;  we  may  even  suppose  q  2  because  when  q  =  1  , 
equation  (19.3)  has  no  solution  . 

Since 


one  finds  p  ^  5  1  there  are  four  cases  to  consider  ,i.e. 

P  =  2  ,  3  >  4  ,  5  . 

It  is  readily  seen  that  for  p  =  4  »  equation  (19.2)  has 
no  integral  solution  .  There  remain  the  three  distinct  solutions 


P«2,q  =  5  ;  p  =  3  »q  =  3  ;p  =  5fq==2  . 

Therefore  ,  when  n  =  1  ,  the  only  integral  solutions 
(p,q)  of  (19.2)  are 

(19.4)  (0,-3)  ,  (2,5)  ,  (3,3)  ,  (5,2)  . 

ii.  n  =«  oo  .  Equation  (19.1)  reduces  to  p  q  -  1  and  has 
the  only  integral  solution  p=l,q=l. 

iii.  n  ^  2  .  For  p  =  0  ,  we  have  q  =  -(n+2)  . 

Suppose  p  5  0  ;  equation  (19.1)  rewritten  as 


shows  that  q  1  , 

From  (19,1)  ,  it  also  follows  that 


since  n;^2  ,  q!^l  ,  qn;^2  ,  one  has  p  ^  5  •  Thus  ,  we  have 
to  consider  five  possibilities  ,  namely  p=l  ,2, 3, 4, 5. 

1  1  O 

a,  p  =  1  ;  the  equation  (19.1)  is  now  ^  ~  1 

n  q  n  q 

and  has  been  considered  in  i  ;  its  only  distinct  integral  solu¬ 
tions  are 

(n,q)  =  (00,1)  ,  (5,2)  ,  (3,3)  ,  (2,5)  . 

b,  p  =  2  ;  equation  (19.1)  becomes  ^  ^  ^  =  2 

n  q  n  q 

and  has  only  the  two  integral  solutions 
(n,q)  =  (5,1)  ,  (2,2)  . 

c,  p  =  3  ;  equation  (19.1)  becomes  ^  ^  ^  =  3 

n  q  n  q 
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and  has  only  the  integral  solution  (n,q)  -(3»1)  . 

d.  p  =  4  ;  equation  (19.1)  has  no  integral  solution  . 

e.  p  =  5  ;  equation  (19,1)  is  ^  ^  =  5 

n  q  q  n 

and  has  only  the  integral  solution  {n,q)  *  (2,1)  . 

iv,  n  =  -2  ;  equation  (19»l)  isp+q+2pq=0  and  its 
only  integral  solution  isp=0,q=0. 

V.  n  =  -  3  ;  equation  (19.1)  isp+q+3pq=l  and  has 
only  the  two  integral  solutions  p=l  ,  q=0  ;p=0  ,q=l. 

vi.  n  =  -  4  ;  equation  (19.1)  isp+q+4P9=2  and  has 
only  the  two  integral  solutions  p=0  ,  q=2  ;  p=2  ,  q=0, 

vii.  n  =  -  6  ;  equation  (19,1)  isp+q+6pq=4  and 
has  only  the  two  integral  solutions  p=0,q=4;p“4» 
q  =  0  . 


We  summarize  the  distinct  integral  solutions  of  (I9.I)  ; 


we  give  the  values  of  n 

n 

and  the  corresponding  values  of  (p>q).  i 

1 

!  (p>q) 

00 

1,1 

1 

2  ,  5  ;  3  ,  3 

2 

1  ,  5  ;  2  ,  2 

3 

1  ,  3 

5 

1,2  ; 

-2  1 

j 

0,0 

0,1  1 

-4  i 

i  0,2  1 

-6 

0,4  1 

n  2 

i  j 

1  0  ,  “ ( n+2 ) 

I 
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V,  Equations  of  the  type  A{x,y)  *=  0  . 

20.  Set 

(20.1)  F(x,y)  =ey+ey^+fy+g  . 

The  stable  equations  of  this  type  are  of  the  form 


(20.2) 

y 

=  F(x  ,  y)  , 

(20.3) 

y 

• 

==^ayy  +  F(x,y)  , 

(20.4) 

y 

•  2 

=  a  y  y  -  _  y^  +  F(x  ,  y)  , 

9 

(20.5) 

y 

=  b  y^  +  F(x  ,  y)  , 

(20.6) 

y 

•  S3 

=ayy+a  y^+F(x,y)  . 

Equation  (20,2)  has  parametric  poles  of  the  second  order 
(P2)  ;  the  other  equations  have  movable  poles  of  the  first 
order  (PI)  .  Therefore  ,  the  stable  equations  of  the  type 
A(x,y)  5  0  may  be  restricted  to  five  distinct  equations  which 
may  be  represented  symbolically  by 

P2  ;  (PI  ,  p  =  1)  ;  (PI  ;  p  =  0  ,  p  =  -3)  ; 

(PI  ;  p  =  3  ,  p  =  3)  ;  (PI  ;  P  =  2  »  p  =  5) 

Equations  (20,2-6)  belong  to  the  type 
(20.7)  y*ayy+by^+  F(x  ,  y)  . 

The  general  transformation  £  see  (16.3)"] 

(20,6)  y  =  A.  (x)  u  +  H  (x)  ,  t  =  (p(x) 

does  not  alter  the  form  of  equation  (20,7)  which  becomes 


I 
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u»*  = 
where  A  ,  B 

(20.9) 

i 

I 

21.  Equation  (20.2).  We  explain  the  particulars  of  our 
method  on  this  equation  which  is  of  class  P2. 

According  to  section  17  p  set  * 

(21.1)  y  =  p  z  =, 1  +  u  z  . 

2 

z 

For  equation  (20.2)  p  one  has  n  =  1  and  et  *  6  ,  p  =  5  i 

[  see  (17.0)]  ,  ! 

It  follows  from  (20.9)  that  a  transformation  y  =  X  u  may  | 

be  chosen  so  that  s=  1  ,  i.e.  e=6  ;  equation  (20.2)  gives  ! 

I 

(21.2)  zu  =  5u.  +  c-  ^  z  -  ^  ^  +cuz+2zu^  . 

2  2 

To  apply  our  main  theorem  of  stability  (  theorem  IV)  ,  we 
must  set  see  section  6  J 

(21.3)  u  =  P(x  ;  z)  +  z^  V  ,  j 


Auu»+Bu^  +  Cu»+Eu^  +  Fu  +  G  , 

,C,E,F,G  are  given  by 

(pA  =  aX  ,  ep^B*bX^  , 

<pC=(c  +  aji  -2A  -  2$)  , 

E  =  [  a  X  +  3  b  \  IX  +  e  X]  , 

<p^F  =  a(  A|i  +  |i  )+2en  +  cA  +f+3b  jx  ^  ^ 

X 

Xqt^G  »  an  n  +b  n^+c  H  +e  +f|i  +g-  Ji 


where  P  is  a  polynomial  in  z  of  degree  4  .  To  determine  its 
coefficients  ,  we  derive  from  (21.2)  p  (21.3)  p  five  relations; 
a  sixth  relation  giving  the  condition  for  stability  is  then 
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deduced  from  theorem  IV  , 

To  simplify  the  problem  ,  It  is  most  convenient  to 
transform  equation  (20.2)  by  a  general  transformation  T(  H,  <p) 
into  an  equation  of  the  same  form  for  which  c  =  f  =  0  *  This 
is  always  possible  ,  for  one  has  to  determine  k  ,  n ,  <p  so  as 
to  satisfy  ^  see  (20.9)"^  E  =  6,C  =  0,F  =  0;  then  k, 

(p  are  given  by 

« 

+  A-2$  -  0  ,  e-2A-®=0 

e 

and  are  determined  by  quadratures  ;  n  is  determined  by 
2eji  +  cA  +  f-  k  =  0  , 

With  these  simplifications  ,  our  problem  is  readily  solved  , 

Indeed  ,  equation  (20,2)  becomes 

y  =  6  y^  +  g(x) 

and  (21,2)  reduces  to 

z  u  =  5  u  -  ^  +  2  z  u^  . 

2 

3  h  ‘S 

Now  ,  set  u  =  a  z  +  P  z  +  z  v  and  determine  «  and  P  by 
•  • 
ifa*g»4P  —St  the  condition  for  stability  is  then  p  =  0 

(theorem  IV)  so  that  g(x)  =  0  or  g(x)  =  K  x  +  H  . 

The  stable  equations  (20,2)  are  thus  reducible  to 

(21.4)  y=6y^+Kx+H, 

By  trivial  changes  of  variables  ,  this  equation  may  be 

brought  into  one  of  the  following  canonical  forms 
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i.  when  K  =  H  =  0  , 

(21.5)  y  =  6  ; 

ii.  when  K  »  0  ,  H  ^  0  , 

(21.6)  y  *  6  ^ 


ill.  when  K  ^  0  , 

(21.7)  y  =  6  y^  +  X  . 

Equations  (21.5)  and  (21.6)  may  be  integrated  by  means 
of  elliptic  functions  so  that  their  general  integral  is  a 
one- valued  function  of  x  .  Equation  (21.7)  is  not  integrable 
in  terms  of  classical  transcendents  ;  its  solution  is  one  of 
the  new  transcendents  discovered  by  P.  Painlev4  .  (see  table  I 
equation  I  ) , 


22.  We  consider  now  equations  (20.3-6)  all  of  class  PI  . 

six) 

Equation  20.3.  Set  y  =  _______  where  s(x)  is  given  by 

z 

as  +  2  =  0  ;  note  that  p  =  1  .  A  transformation  T(  X  ,  ,  <f  ~x) 

may  be  chosen  so  as  to  secure  a=-2(i.e.  s=  1)  ,  and  esc. 
Then  ,  one  has 

*  ■  2  2 
ZU  =  U-  fz-gZ  +  CUZ  +  ZU  . 

On  setting  u  =  z  v  ,  one  obtains  the  condition  for  stability 
f  =  0  .  The  stable  equations  (20.3)  are  thus  reducible  to 
(22.1)  y=-2yy+c(y+y^)+g 

and  are  equivalent  to  the  differential  system 

*  2  ’ 

z  =  l  +  zv  ,  v=cv-g 
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The  general  solution  of  (22.1)  is  given  by 

*  2 

w  =  y  +  y  , 

where 

w  =  c  w  +  g  . 

Therefore  ,  w  is  determined  by  quadratures  and  y  is 
given  by  a  Riccati  equation  ;  equation  (22.1)  is  thus  stable  . 

23.  Equation  (20,4) .For  this  equation  ,  p  =  0  ,  as  =  -3 ; 
p  =  -3  ,  as  =  -6  ,  We  have  only  to  consider  p  =»  0  since  a 
negative  p  does  not  result  in  a  condition  of  stability  .  Through 
a  transformation  T(  X  )  ,  one  may  assxjme  a  =  -3  and  s  =  1  , 

The  associated  equation  in  u  is 

•  2  2 

zu  =  C“e-fz-gz  +  CZU  +  ZU 

which  gives  the  condition  for  stability  c  =  e  , 

The  stable  equations  (20.4)  are  of  the  form 

(23.1)  y=-3yy-y^+c(y+y^)+fy+g. 

On  setting  y  =  ^  ,  this  equation  reduces  to  the  linear 

V 

equation  of  the  third  order 

* 

(23.2)  v-cv-fv-gv=0 

which  was  investigated  by  E,  Vessiot  (Ann.  Fac.  Toulouse  1695»F2). 

24.  Equation  (20,5).  We  have  p=3,bs=2.A  trans¬ 
formation  T(  X  ,  ji  ,  )  may  be  chosen  so  that  b  =  2  (i.e, 

s  =  +  1  )  and  e  =  c  =  0  , 

Indeed  ,  one  has  to  determine  X  ,  ,  T  by  (  set  again 


55 


^  *  A  ,  ®  )  , 

X  <p 

b\^  =  2  q)^  ,  c-2A-2®=0  ,  3b|x  +  e  >=  0 

X  and  ^  are  given  by  quadratures  through 

^  +2A  -2$=0  ,  c-2A~2'5*0. 

b 

With  the  values  ass\imed  for  b  ,  c  ,  e  ,  equation  (20.5) 
becomes 

(24.1)  y=2y^+fy+g  ; 

Its  associated  equation  In  u  is 

(24.2)  zu  =  3u.-fz-gsz^  +  zu^  . 

Now  ,  set 

u=Pz+YZ  +z-^v 
and  determine  p  and  y  by 

2p  =f  ,  p=Y-gs  ; 

the  condition  for  stability  is  y  ~  0  ,  Therefore  ,  f  +  2  g  s  =  0; 

* 

because  s  =  +  1  ,  this  relation  splits  up  into  f  =  0  ,  g  =  0  . 
Hence  ,f=Kx+H;gis  also  an  arbitrary  constant  (  say  K^) . 
The  stable  equations  (20.5)  are  thus  reducible  to 

(24.3)  y=2y^+(Kx+H)y+K^  . 

When  K  =  0  ,  this  equation  becomes 

(24.4)  y=2y^+Hy+K^ 

and  is  at  once  integrable  in  terms  of  elliptic  functions  . 

When  K  0  ,  a  trivial  change  of  variable  brings  (24.3) 
into  the  canonical  form 

(24.5)  y=2y^+xy+K2  . 
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This  equation  is  not  integrable  in  terms  of  classical 
transcendents  ;  its  solution  is  one  of  the  new  transcendents 
discovered  by  P,  Painlevd  (see  table  I  ;  equation  II)  , 

25 •  Equation  (20.6).  One  has  p  =  2  ,  as  ~  -1  and  p  =  5  » 
as  =  2  .  A  transfoxTiiation  T(  X  ,  [i  ,  ip  )  may  be  chosen  so  that 
a*-l  ,  e  =  c  =  0  .  Indeed  ,  one  has  to  determine  X  ,  n  ,  by 
aX=-  (p  ,  c+ati-2A  -f=0  , 

aA  +  3bii-fe  =  0 

or  by 

• 

^  +  A  -  $  =  0  ,  c+a  -2  A  -  ®  *0  ,  e+aA  +3bn  *  0  ; 

a 

these  equations  have  one  and  only  one  solution  . 

Equation  (20.6)  assumes  the  form 

(25.1)  y=-yy+y^+fy+g  ; 
the  associated  equation  in  u  is 

(25.2)  zu=(3-s)u-fz-  ^  +  z  u^  . 

s 

a.  First  ,  consider  p  =  2  ,  s  =  1  ;  the  equation  (25*2) 
becomes 

(25.3)  zu=2u-fz-gz^+zu^  . 

On  setting  u  =  |3z  +  z  v,  one  obtains  the  relations 

P  =  f  ,  P  =  -  g 

and  hence  the  condition  for  stability  f  +  g  =  0  . 

To  simplify  the  notations  ,  we  set 
f  =  -12V  ,  g=12V. 
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b.  It  follows  from  a  ,  that  the  stable  equations  (25.1) 
are  of  the  form 

(25.5)  y  *  -  y  y  +  y3  -  12  V  y  +  12  V  . 

Now  ,  consider  p  =  5>s==“2;  the  condition  for  the  stabili 
ty  of  (25.5)  gives  rise  to  a  condition  on  V  . 

The  equation  in  u  associated  with  (25.5)  ie 
2u=5u+12Vz+6Vz^+zu^  . 


P  z  +  Y  z^  +  e  z^  +  z^  V 


On  setting 
u 

one  obtains  the  relations 

4P  =  -12V  ,  p’  =  6W3y 

Y*2  ,  5+pY*e 

e+Y^  +  2p?.  =0  ; 

these  relations  easily  give 

P  =  -3V  ,  Y=-3V 

3 


5  =  -  3  y 

2 


e  =  - 


V  +  9  V  V 


and  hence  the  condition  for  stability 

( It )  *  2  ** 

-  12  -  12  V  V  =  0 


or 


(  V  -  6  V^)  =  0 


dx 


and  finally 

(25-.6)  V  =  6V^  +  Kx  +  H  . 

Thus  ,  the  stable  equations  (25.1)  are  of  the  form  (25.5) 
where  V  is  a  solution  of  (25.6)  ,  i.e.  ,  an  elliptic  function 


58 


or  a  solution  of  the  irreducible  Painleve’s  equation  I, 

c.  That  the  equation  (25.5)  where  V  is  given  by  (25.6)  t  i® 
really  a  stable  equation  remains  to  be  shown  .  To  do  this  , 
we  employ  a  method  which  shall  often  be  used  in  what  follows. 

Observe  that  the  stability  condition  corresponding  to 
p  *=  2  (see  a)  is  readily  obtained  and  brings  equation  (25.1) 
into  the  form  (25.5)  *  Moreover  ,  one  has 

y  =  ^  ,  z=l+zu  ,  u=-12Vz+z^v 

z 


so  that  ^  +  u  ;  it  then  follows  that  ^  +  y  =  u 

y  z  y 

is  regular  at  the  parametric  poles  corresponding  to  p  ==  2  and 
has  simple  poles  at  the  parametric  poles  corresponding  to  p  =  5. 
Accordingly  ,  set 


(25.7)  y  +  y‘=^  =  V  ; 

V  has  double  poles  at  the  parametric  poles  corresponding  to 
p  =  5  ;  it  is  thus  expected  that  the  stability  condition  for 

V  is  connected  with  equation  (20.1)  . 

Indeed  ,  the  equation  (25.5)  is  equivalent  to  the  differen¬ 
tial  system 


(25.8) 


y  +  y 


|^v  =  vy-12Vy  +  12V  . 

Hence  ,  v  is  a  rational  function  of  y  ,  y  and  ,  vice  versa  , 
y  is  a  rational  function  of  v  ,  v  ;  therefore  ,  if  y  is  stable  , 
V  is  also  stable  and  conversely  , 

On  eliminating  y  between  the  relations  (25.8)  ,  one  finds 
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(25.9)  V  =  -  12  V  V  +  12  V 

r  write  ^  +  y  =  ^  ,  and  use  the  logarithmic  derivatives 

L  y 


of  y 


y 

a  • 

V  -  12  V 


1  • 


V  -  12  V 

The  transformation  v  -  6  V  =  6  w  brings  (25.9)  into  the 


canonical  form  see  !  21-] 


w  =  6  +  V  -  6 

in  order  that  this  equation  be  stable  ,  one  must  have  [|see  (21,4 

V=6V^  +  Kx  +  H 

and 


(25.10) 


w=6w‘^  +  Kx  +  H  . 


Therefore  ,  the  general  solution  of 

y  =  -  y  y  +  y^  -  12  V  y  +  12  V  , 

where  V  is  a  particular  solution  of  (25.10)  ,  is 

•  • 

w  -  V 


(25.11) 


y  = 


w  -  V 

w  is  a  solution  of  (25.10)  distinct  from  V  . 


26.  To  sum  up  ,  the  following  set  of  equations  may  be 
considered  as  a  set  of  canonical  equations  of  the  type  A(x,y)=0. 

i.  y  =  0  ; 

ii.  (P2  ;  p  =  5) 

a.  y=6y^+K 
integrable  by  elliptic  functions  ; 

b,  y  =  6  y*^  +  X 
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not  integrable  in  terms  of  classical  transcendents  ; 

iii.  (PI  ;  p  =  1) 

y=.-2yy  +  c{y+y^)+g 
integrable  by  quadratures  ; 

iv.  (PI  ;  p  =  0) 

y=-3yy-y^+c(y+y^)  +fy+g 
may  be  reduced  to  a  linear  equation  of  the  third  order  ; 
•V.  (PI  ;  p  =  3  ,  3) 

a.  y=2y3+Ky+H 
integrable  by  elliptic  functions  ; 

b.  y  »  2  y^  +  X  y  +  K 

not  integrable  in  terms  of  classical  transcendents  ; 

Vi.  (PI  ;  p  =  2  ,  5  ) 

y  =  -  y  y  +  y3  -  12  y  +  12  , 

where  is  a  solution  of  ii  ,  a  or  ii  ,  b  , 

Integration  : 


where  V  and  are  distinct  solutions  of  ii  ,  a  or  ii  ,  b  . 
In  iii  and  iv  ,  c  ,  f  ,  g  are  arbitrary  analytic  functions 
of  X  , 

To  obtain  the  most  general  stable  equations  of  the  type 


A(x,y)  =  0  ,  one  has  to  use  a  general  transformation  T(  X  ,  n  ,  (p  ) 
where  \  ,  |i  ,  q>  are  arbitrary  analytic  functions  of  x  . 


f 
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VI,  Equations  of  the  type  A(x,y)  «  ^ 

y 


27,  The  stable  equations  of  this  type  are  of  the  form 

*2 

(27,1)  y  -  JL_  =  a  y  y  +  a 


^  +  b  y^  +  +  F(x,y)  , 


y  y  y 

where 

(27,2) 

the  a’s  and  b’s  satisfy  conditions  to  be  given  later  on  , 

The  substitution  y  v  *  1  transforms  the  equation  (27,1) 


F(x,y)  =  a^  y  +  bj  y*^  +  bg  y  +  b^ 


into  an  equation  of  the  eeime  form  ,  namely 

(27,3)  V  -  ^  ^  ^  "  — — 

V  V  ^  V 


^  +  a,v  -  b^v^-b_v-b. 


therefore  ,  the  stability  conditions  corresponding  to  the 
parametric  zeros  of  y(x)  may  be  deduced  from  the  stability 
conditions  produced  by  the  parametric  poles  ,  according  to  the 
table  of  equivalence 


(y)  : 

(v)  ; 

The  transfonnation  y  -  \  (x)  u  ,  t  »=  9  (x)  does  not  alter 


^2  ^1 


U 

-'’3  -'’S  -‘’1 


"4 

-b 


the  form  of  equation  (27.1)  which  becomes 


u«  - 


u 


,2 


=  Au»u+A, 


u 


B 


u 


+  Bu-^  +  +  A,u»  +  BiU*"  +  B«  u  +  B_  , 

,,  ,,  X  1  2  J 


U 


U 


where  on  setting  A 


X 


,  one  has 
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•  *  m  * 

(p  A  =  a  X.  ,  ^  ^2  *  ^"2  ' 

« 

<P 

=  bX^  ,  <p^Bj=b^\  +a  X  ,  Ji^B^^bg+a^  A  -  A  , 

To  simplify  the  problem  ,  one  may  choose  ^ ,  T  so  that 

a^  =  bg  =  0  .  One  has  only  to  determine  X  ,  (p  in  order  that 

.•  •  •  “1 

(p*-  a^  cp  ,  A  -  a^A  -  =  oj  , 


26 «  The  general  solution  of  (27.1)  may  have  simple  or 

double  parametric  poles  and  simple  or  double  paraonetrlc  zeros  ; 

we  consider  these  cases  separately  , 

1.  Double  parametric  poles  ,  One  has  a  =  b  *  0  •  The  transfor- 
s  (x)  * 

matlon  y  =  _______  ,  z  =  1  +  u  z  ,  where  s  is  given  by  b.s  *  2  , 

2  ^ 
z 

changes  (27.1)  into 


zu 


u  + 


Ig  s  +  ^2  z^ 

s 


b,  5  a^  3 

4  z^  +  2  z^  u 


2s 

Because  p  =  1  ,  one  sets  u  =  z  v  and  finds  the  condition  for 


stability 


Ig  s  =  0  or 


Ig  bj^  =  0  (note  that  b^^s  =  2) 


dx 


dx^ 


■and  finally 

(26.1) 


e 


K^x 
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li.  Double  parametric  zeros  .  One  has  =•  0  ;  according 

to  the  table  of  equivalence  ,  the  condition  for  stability  is 


(28.2) 


’03  -  Hg  e 


iii.  Simple  parametric  poles:.  ,  Set  y  =  ,z=»l+uz, 

z 

P 

where  s  is  given  by  bs  -  as  -  1  =  0  and  p  =  1  +  as  ;  equation 
(27»1)  becomes 


(28  .3)  zu=pu-b^s-as  +^2“  ^  ^  ^ 


-  z^+  ^2  z^u 


with 

(28.4) 


A  -  _  Ig  s  +  °2 

2  s 


dx 


It  follows  from  (l6,15)  and  (I6.I6)  that  we  have  to  consider 
two  cases  according  as  to  whether  a=90  ,  b^f^Oora^O  ,  b  =  0. 

a)  a  ^  0  ,  b  =  0  ;  then  p=0,as+l*0.  The  condition  for 
stability  is  bj^s  +  as  =  0  or 

(28.5)  b^  =  a  . 

p 

b) a=0,b^0;  then  p  =  1  ,  bs  =  1  .  Set  u  =  «  +  z  v  in 
(28,3)  ;  bhe  condition  for  stability  is  given  by 

a=b^s  ,  a=A  , 

or  by 

(28.6)  _  (b  s)  =  A  . 

dx  ^ 

From  bs^  =  1  ,  it  follows  s  *  eb”^^^  ,  e  +  1  and  from 


(28.4)  and  (28.5) 


f 
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b,  b 


dx 

^  dx' 

or 

(26.7) 

d 

dx 

so  that 

(28.6) 

«  Kx 

b  =  H  e  . 

-1/2  .  .  J_  ig  fc.  +  e 

2 


Ig  b  »  0 


dx 


iv.  Simple  parametric  zeros  .  From  111  a  ,  b  and  the  table 
of  equivalence  ,  we  easily  obtain  the  following  conditions 
for  stability  . 

a) ag^0,b^*0;p  =  0,  ^  condition  for 

stability  is 

(2e.9)  b^  4-  =  0  . 

b)  a^  =  0  ,  i  P  =  1  »  ^4®^  +1  =  0.  The  condition  for 

stability  Is  given  by  [^see  { 26.6)1 


(26.10) 


d 

dx 


(b^s)  + 


Ig  s  +  ^ 
s 


0  . 


29,  It  follows  from  the  preceding  paragraphs  that  the 
stable  equations  of  the  type  A(x,y)  =  ^  ,  may  be  restricted 

y 

to  six  distinct  equations  which  may  be  represented  symbolically 
by 

( {P2;Z2)  ;  (F2;ZI,p  =  0  or  p  =  1)  ; 

(29.1)  \ 

j  (Pl,p  =*  0  ;  Z1  ,  p  =  0  or  p  =  1)  ;  (Pl,p  *  l;Zl,p  =  1). 
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Indeed  ,  the  other  equations  are  deducible  from  one  or 
the  other  of  the  equations  (29* 1)  hy  setting  y  =  u”^  . 

We  consider  separately  each  of  the  equations  (29.1)  • 


30.  P2:Z2.  One  has  a=b=a_=b.  =0.  The  conditions 
■  2  4 

for  stability  are  given  by  (26.1)  for  P2  and  by  (26,2)  for  Z2, 
The  equation  is 


(30.1) 


2  K  X 

y  -  e  ^  e 

y 


and  is  a  particular  case  of  the  Irreducible  equation  III  (see 
table  I  ) , 


31.  P2;Zl.p  =»  0  ,  One  has  a  =  b  =  b^»0>  0  . 

For  P2  ,  the  condition  for  stability  is  given  by  (26,1)  . 

For  Z1  ,  p  =  0  ,  this  condition  is  given  by  (26.9)  .  Set  a^  *  q, 
b^  =  -  q  and  write  the  equation 

(31.1)  . 

y  y 


Kx 

Now  ,  def  ine  y  =  X  u  ,  where  X  is  given  by  X  H  e  =*  1  and  set 

•  • 

r  =  ;  then  because  r  =  ^  ^  ,  (31.1)  becomes 

X  ^  ^  X  x^ 

(31.2)  ii  -  _  =  r  +  u^  -  r  . 

u  u 

*  • 

•  u.  4"  R 

For  convenience  ,  set  r  =  R  and  v  *  i  equation  (31,2) 

u 

is  equivalent  to  the  differential  system 


V  =  u 
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On  eliminating  u  between  these  relations  ,  one  finds  v  +  R  =  vv 

*  18 

so  that  V  is  determined  by  the  Riccati  equation  v  +  R  =®  _____  v 

2 

Therefore  ,  u  is  stable  , 


32.  P2;Zl.p  =  1  .  One  has  a  =  b  =  =  0  ,  b^  ^  0  . 

For  P2  ,  the  condition  for  stability  is  given  by  (26.1)  . 

For  Z1  ,  p  =  1  ,  one  has  s  =  b~^^^  e  ,  e  =  +  i  so  that  (26.10) 
becomes 


dx 


b:'/"  b,  -  .i 


Ig  b^  =  0 


dx‘ 


and 


b  ,  =  H,  e 
4  4 


2K3X 


b3  =  H3  e 


K^x 


2K2X 


Equation  (3I.I)  is  written 

K,x  _  K_x  jj 

(32.1)  y  -  7 _  =  H.  e  ^  y"^  +  e  + 

—  1  3  y 

and  is  a  particular  case  of  the  irreducible  equation  III 

{  see  table  I  )  . 


33.  PI  .p=0  ;Z1  .p=0  .  One  has  b  =  b^  =  0  ,  aa^  ^  0. 

•  r 

The  conditions  for  stability  are  b^  =  a  ,  b3  +  a^  =  0  see 
(26,5)  and  (26.9)"^  so  that  the  equation  takes  the  form 

.0  • 

••  *  V  *  2  * 

(33.1)  y  -  J, _  =  a  y  y  +  a  +  a  y  -  a 

y  y 

•  2 

Integration  .  It  is  easily  seen  that  y  -  a  y  =  v  is  regular 
for  the  poles  PI  ,  p  =  0  ;  the  equation  (33.1)  gives 
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V  +  (  V  +  ag) 

y 

or  V  +  a^  =  K  y  ,  Therefore  ,  equation  (33.1)  is  equivalent 

to  the  differential  system 

*  2 

y-ay  =  v  ,  v  +  a2  =  Ky 

so  that  y  is  determined  by  the  Riccati  equation 

*  2 

y  =  ay  fKy-a^ 
and  is  stable  . 

34.  PI  .  P=0  ;  Z1  .  p=l  .  For  PI  ,  one  has  a  ^  0  , 
b  =  0, as  +  1=0  and  the  condition  for  stability  is  b^  =  a  , 

For  Z1  ,  one  has  b,  0  ,  b,  s  +  1  =  0  :  a  transformation 
y  —»  X  y  enables  one  to  assxime  =  -  lors  =  +  l  ;  the  condi- 
tion  for  stability  is  then  b^  +  a  =  0  . 

• 

For  convenience  ,  set  b^  =  q  so  that  a  =  -q,bj^  =  -q; 

the  equation  then  assumes  the  form 

•  2  ,  , 

(34.1)  y  *  _ =  -  q  y  y  -  -  q  y^  +  q  • 

y  y 

The  transformation  y  u  =  1  brings  this  equation  to 

(34.2)  u  -  ^  =  -  q  ^  +  u^  -  q  u^  +  q 

u  u 

which  belongs  to  the  class  (Zl  ,p=0;Pl,p=l)  . 

To  integrate  the  equation  (34.1)  ,  observe  that 

y  *  1 

V  =  ^  +  q  y  -  _  +  q 

y  y 

is  regular  for  the  parametric  poles  and  for  the  parametric 
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zeros  corresponding  to  s  =  1  , 

Equation  (34«1)  is  equivalent  to  the  differential  system 
y  +  qy^-l  +  qy  =  vy  ,  yv  =  v  . 

On  eliminating  y  between  these  relations  ,  one  finds 

•  •  • 

v  +  vv-qv-qv  =  0 

and  V  is  a  solution  of  the  Riccati  equation 
*  12 

V  +  v"-  -  q  V  =  H  . 

2 

w 

Set  V  =  2  ____  ;  then  v  is  a  solution  of 
w 

••  *  H 

w  -  q  w  «  _____  w  =  0  . 

2 

Therefore  y(x)  is  determined  by 


w  w 


and  is  stable  , 


35.  PI  .  P  =  1  ;  Z1  .  p  »  1  .  One  has  a*ag=0,bb^^0. 
According  to  (28,7)  and  (28.8)  and  the  table  of  equivalence  , 


the  conditions  for  stability  are 

2KiX 


for  FI 
for  Z1 


b  =  e 


2K3X 


bi  =  e 


K^x 


e 


K^x 


The  equation  is 
.*2 


(35.1) 


2K.X  _  K.x  p  K—X  „  2K^x 

y  -  _  =  H,  e  ^  Y  tHp  e  ^  y  +H  e  e 

y  ^  V 


A  transformation  y  \  (x)  y  ,  where  \  is  defined  by 
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(K^-K  )x 

X'^e  =1  and  m  by  +  Kg  =  2  m  ,  brings  equation  (35»l) 

into  the  fonn 


(35.2)  y 


+  e”“  (  +  Hj)  . 


Now  ,  we  consider  two  cases  according  as  to  whether  m  =  0  or 
m  0  , 

i.  m  =  0  .  The  equation  (35 .2)  is 

*2 

(35.3)  y  -  ^  =  H  +  ^4  +  H  y^  +  H 

y  y 


and  may  be  integrated  by  a  process  which  may  be  considered  as 
a  method  of  "  variation  of  parameters  ”,  As  this  method  will  be 
used  often  later  on  ,  we  shall  explain  its  particulars  on 
equation  (35.3)  . 

y^ 

The  general  solution  of  y  -  ^  =  0  is  y  =  u  y  where 

y 

u  is  an  arbitrary  constant  . 

Suppose  now  that  u  is  a  function  of  y  and  substitute  y  *  u  y 
into  (35.3)  ;  one  obtains 


u  y  =  H,  y^  +  y^  +  H  + 

^  y 


or 


^  H, 

u  =  H.  y  +  Hj,  +  _ ^  ^ 

dy  ? 


Therefore  , 

u^  =  y^  +  2  H  y  -  ^^3  -  \  +  H- 

y  „2  ^ 


Finally  ,  y(x)  is  given  by 
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+  2  y^  +  y^  -  2  y  - 


(35.4) 

and  Is  an  elliptic  function  of  x  ;  thus  ,  y{x)  is  stable  , 


ii,  m  ^  0  .  A  linear  transformation  on  x  transforms  equa¬ 
tion  (35.2)  into 

*2 


(35.5) 


y  - 


2x  ( 
=  e  I 


.)  -*  ( 


y-  +  H3 


)  ' 


this  equation  is  not  Integrable  in  terms  of  classical  transcen¬ 
dents  and-  is  equation  III  of  Table  I  , 


A  transformation  t  =  «f>  (x)  = 


mx 


brings  equation 


m 


?  _ 


(35.2)  into  the  other  canonical  form 

(35.6)  y  - 

[_  Note  that  0  ^  . 


£_  +  y^  + 

X  y 


mx 


(Hg  y^  +  H3  ) 


36,  To  sum  up  ,  the  following  set  of  equations  may  be 
considered  as  a  set  of  canonical  equations  of  the  type  A(x,y)*=__^ 


i. 


*2 

y  -  y  =  0 


ii.  (P2;Zl,p  =  0) 

*2 

y 

y  -  = 

y 


=  r  +  y^  -  r 


reducible  to  a  Riccati  equation  ; 


iii.  (Pl,p  =  0;Zl,p  *  0) 

•z 

y  -  ^  =  a  y  y  +  a 


y  *  2  * 

2  _  +  a  y  -  a^ 


reducible  to  a  Riccati  equation  ; 
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iv.  (Pl,p  =  0;Zl,p  *  1) 


y  -  JL _ =  -  q  y  y  -  _ - 

y  y 

reducible  to  two  Riccatl  equations  ; 

••  y^  3  y  2 

y  -  ^  =  ay-^  +  ^  +  py 

y  y 

integrable  by  elliptic  functions  ; 

*2 

y  -  ^  (  a  y^  +  0  I  + 


(“  *  -^) 


or 


^  +  a  y^  +  6  + 

X  y 


•*  2 

q  y  +  q 


+  Y 

(  P  +  r  ) 

. -  {  p  y^  +  Y  ^ 

m  X 


where  the  constants  a  ,  p  ,  Y  >  ^  eventually  be  zero  ; 

m  0  is  also  a  constant  . 
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VII. 


Equations  of  the  type  A 


1 


37.  The 
(37.1)  y  -  ^ 


stable  equations  of  this  type  are 

1  -  ^  \  ^  =  ayy  +  a>  ^  +  by^ 

n  /  y  y 


of  the  form 
+  +  F(x,y)  , 

y 


where 

(37.2)  F(x,y)  =  a^  y  +  b^  y^  +  b^y  +  b^  ; 


the  a’s  and  b’s  satisfy  conditions  to  be  given  later  on  . 

The  substitution  yv  =  1  transforms  equation  (37.1)  into 
an  equation  of  the  same  form  ,  namely 


( 1+1-1 

*p 

i  ''  ■( 

a  v+a  +_^\ 

’  3  1.  2  ,  ,  b 

v-b.v^-b_v  -b_v-b,-  ; 

n  ' 

V  ' 

k  3  Z  1  — 

therefore  ,  the  conditions  for  stability  corresponding  to  the 
parametric  zeros  of  y(x)  may  be  deduced  from  the  conditions 
for  stability  produced  by  the  parametric  poles  ,  according  to 
the  table  of  equivalence 

(y)  :  n  a  b  b^  b^  b^  b^^ 

(v)  :  -n  a^  a^  a  -b^  -b^  -b^  -b^  -b 


This  table  shows  that  without  loss  of  generality  ,  we  may  res¬ 
trict  n  to  be  positive  ,  and  accordingly  do  so  • 

The  transformation  y  =  X  (x)  u  ,  t  *  <p  (x)  does  not  alter 
the  form  of  equation  (37.1)  which  becomes 
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u" 


-^1  -  ^  j  »^Au-fA^+  ^2  j  u»+Bu^+B^U^+B2U+B2+_^ 


+Jk. 


where 


<p  A  =  a  X 


(p^B  =  b  X  ^  , 


n 


(jAj^=a^-  _  A  -  ® 

<p^  =  b^  X  +  a  X 


X  cp  A-  =  a 


*2  *  1  2 
<p  Bg  —  bg  +  a^A  —  A  —  ____  A  > 


Xqj^B^  =  b^  +  a^  A 


,  2  *2n  _  , 

X  -  b^  . 


For  future  use  >  we  also  note  the  following  system  of 

equations  equivalent  to  equation  (37.1)  ;  on  setting  y  =  ^  , 

z 


z  *=  1  +  uz  ,  one  has 


z  \  n 


+  as  -  bs' 


1  +  ___  +  as  1  u 
n  ■  / 


•  *  2 

(37,4)  +a^-b^s-as-  ^  ^  +Az-  ^  b^+a^  ^ 

ri  s  s  •  S’ 

2s\  .3.o2,z< 

+  /  a  -  - _ 1  z  u  +  2  z  u  +  u 

V  ns/  s  -  "n 


b,  z' 


with 


(37.5)  A  =  ^  1  -  J_j 


-  a. 


2  -  b 


2  * 


Two  equations  corresponding  to  one  another  by  the  transformation 
yu  =  1  are  not  considered  as  distinct  ;  therefore  ,  we  have 
only  to  consider  four  classes  of  equations  represented  symboli¬ 
cally  by  (P2;Z2)  ,  (P1;Z2)  ,  (P2;Zl)  ,  (P1;Z1)  . 

We  consider  each  class  separately  and  summarize  the  various 
stability  conditions  already  obtained  above  (  see  ^  17)  with 
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tables  where  the  corresponding  values  of  p  are  also  Indicated, 


P2 ! 22 .  For  P2  ,  one  has  a  =  b  =  0  and  for  Z2  ,  a^  =  b^  =  0  . 
The  various  cases  to  be  considered  are  given  in  the  follov/ing 
table  : 


n 

P2 

t 

P  : 

Z2 

P 

n 

bi  =  0 

_  i 

1 

X 

bj  -  0- 

- 

2 

b^  7^  0 

3 

b^  =  0 

4 

b^  7^  0 

2 

b3  7^  0 

0 

4 

bl  =  0 

1 

! 

b-  7^  0 

0 

Thus 

,  the  stable  equations  of  clas 

5  (P2,Z2)  are  of  the  form 

E.l 

y-  (i 

=  F(x,y) 

9 

E.2 

y- 

2 

*2 

^  +  F(x,y) 

9 

E.3 

y- J_ 

h 

^  +  b^  +  F 

! 

(x,y)  , 

E.4. 

y- J_ 

h. 

JL.  =  b  3 

V 

F(x,y) 

9 

1 

where 

F(x,y)  =  aj^  y  +  bg  y 

PI ;Z2.  Because  of  the  double  parametric  zeros  ,  one  has 
Sg  =  b^  =  0  ;  in  addition  ,  we  have  also  to  consider  ,  for 
n  /  0  ,  b^  =  0  and  for  n  =  4  ,  b^  ^  0 
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The  various  cases  to  be  considered  are  given  in  the 
following  table  : 


1 

PI  Pi 

1 

Z2  1 

P 

. 

n 

na'^4-b{n+2)  =  0  |0,-(n+2)  j 

*  1  1 

b^  =  0 

- 

2 

2  1 
a'^  =  2b  j  1  »  5  1 

b3  -  0 

- 

2 

a  =  0  ,  b?^0  2,2  ! 

(  j 

b3  =  ° 

- 

3 

3  =  2b  il  ,  3  i 

1  i 

o 

It 

CO 

- 

5 

p  1 

5  =  b  !  1  ,  2  1 

1  1 

b^  =  0 

- 

k 

+  9b  =  0  io,-6  1 

!  i 

b  ^  0 

0 

Thus  ,  the  stable  equations  of  class  (P1;Z2)  are  of  the 


form 


E.5. 

y*  ( 

f 

1  - 

I 

ayy  -  ^  y^  +  F(x,y) 

V 

n  • 

"y 

(n+2)^ 

E.6. 

y  - 

1 

? 

• 

,  =  ayy 

+  _  y^  +  F(x,y)  , 

2 

y 

2 

E.7. 

y  - 

1 

*2 

y 

1 

II 

o' 

+  F(x,y)  , 

2 

y 

E.e. 

y  - 

2 

• 

.  =  ayy 

+  ^  a^  y^  +  F(x,y)  , 

3 

y 

2 

E.9. 

y  - 

4 

• 

.  =  ayy 

+  5  a^  y^  +  F(x,y)  , 

5 

y 

E.IO 

.  y  - 

3 

? 

• 

,  =  ayy 

-  y^  +  b  +  F(x,y) 

4 

y 

9 

where 

F(x,y)  = 

= 

y  + 

bj  y  + 

bg  y  . 

9 
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P2;Z1.  We  assume  n  >  0  .  Because  of  the  double  parametric 
poles  ,  one  has  a  =  b  =?=  o  .  The  various  cases  to  be  considered 
are  given  in  the  following  table 


n 

P2 

1 

1  ^  ' 

1 

Z1  1 

! 

i  ^ 

n  ^  2 

b^  0  1 

- 

na|+(n-2)^bj^=0 

0,n-2 

n  =  2 

b,  *  0 

J. 

_  1 

aj,  =  0  ,  b^  ^  0 

0  ,  0 

4 

bj  ^  0 

2 

a|  +  -  0 

0  ,  2 

2 

b^  ^  0 

3 

1 

a2  =  0  ,  b^  0 

0,0 

Thus  ,  the  stable  equations  of  class  (P2;Z1)  are  of  the 

form 


1 

\  *2 

*  „  j- 

E.ll. 

y  ” 

- 

1  — 

=  a.  y  ^  2  _  +  F(x,y),n?^2, 

\ 

n 

/  y 

(n-2)=y 

E,12. 

y  - 

1 

=  \  +  F(x,y)  , 

“2”  ' 

y 

y 

E.13. 

y  - 

3 

=  a 

y  -  ^2  +  b^y^  +  F(x,y)  , 

y  y  1 

"T"  * 

y 

C 

E.14. 

y  - 

1 

? 

y^  +  ^4  +  F(x,y)  , 

”  ■ 

y 

i 

y 

where 

F(x 

»y)  = 

+  b 

_  y  4- 
Z  ^ 

^3  • 

;Z1 .  The  various  cases  to  be  considered  are  given  in  the 


following  table  ; 
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n 

PI 

1 

P 

1 

Z1 

j 

i 

p 

n.  Z 

na^+(n+27b=0 

0,-(n+2) 

na?+(n-2)^b,  =0 

0,n-2 

2 

a^  +  8  b  =  0 

0  ,  -k 

ag=0  ,  b^ii^O 

0  ,  0 

2 

1 

a^  =  2b 

1  ,  5 

a2=0  ,  b^^?^0 

0  ,  0 

2 

a=0  ,  b  ^  0 

2  ,  2 

a2=0  ,  bj^f^O 

0  ,  0 

3  ! 

3a^  =  2b 

1  ,  3 

0,1 

5 

5a^  =  b 

1,2 

54  +  9b^-0  1 

0  ,  3 

Thus  ,  the  stable  equations  of  class  (P1,Z1)  are  of  the 

form 

E.15.  y  -  (  JL.=  ( a-y+  ^2  \  y-  y^-  ^^2  +F(x.y) , 

^  n  /  y  ^  y  /  (n+2)^  (n-2)^y 

E.16.  y-_  7  =  ayy  -  ^  y^  +  +  F(x,y)  , 

Z  y  8  y 

*2  .  2 

E.I7.  y  -  >1.  ^  =  ayy  +  ^  y^  +  +  F(x,y)  , 

2  y  2  y 

*2 

E.ie.  y  -  -I _  =  b  y^  +  ^24-  +  F(x,y)  , 

2  y  y 

E.I9.  y  -  ^  ^  ^ ay  +  ^2  y+  ^  a^y^-  ^^2  4-F(x.v)  , 

3  y  V  y  y  2  y 

E.20.  y  -  JL  ^  =  f  ay  +  ^2  \  y+’Sa^y^.  ^^2  +  F(x,y)  , 

5  y  \  y  /  9y 

where 

F(x,y)  =  a^^  y  +  bj  y^  +  bg  y  +  b^  . 

Comparing  the  values  of  p  for  PI  in  the  tables  (P1,Z2)  and 
(PljZl)  I  one  is  led  to  study  together  the  equations  E.6  and 
E,17  i  E.7  and  E.I8  ,  E,8  and  E.19  >  E,9  and  E,20. 
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Theorem  IV  gives  additional  conditions  for  stability  , 

We  always  begin  by  considering  the  smallest  values  of  p  . 

36.  Equation  E, 1,  Using  a  transformation  (  ^>9  )  ,  we  may 
assume  =  b^  =  0  ,  Then  ,  equation  E,1  is 


its  general  solution  is  given  by  y  =  (Kx  +  H)^  and  is  stable  . 
On  setting  y  =  2^  ,  equation  (36.I)  reduces  to  z  ==  0  , 

394  Equation  E,2.  For  P2  ,  one  has  p*3,b^s==24.;hya 

transformation  (  ,  <p  )  ,  one  may  assume  a^  ==  0  ,  b^  »  /+  or  s«=l. 

“2  * 

Now  ,  set  y=z  ,z=l+uz  and  determine  u  by 

z  u  =  3  u  -  ^2  z  +  z  u^  . 

2 

To  obtain  the  condition  for  stability  ,  set  u  =  Pz  +  yz  +  z  v 
and  determine  p,y  by 

•  • 

4p=b2  ,  p=Y  ,  Y  =  0  ; 

therefore  , y  is  a  constant  and  b^  =  K  x  +  H  . 

The  stable  equation  E.2  is  thus  of  the  form 

*2 

(39.1)  y  -  ^  +  (Kx  +  H)  y  . 

2  y 

According  to  the  values  of  K  and  H  ,  one  finds  one  or  the 
other  of  the  three  following  canonical  equations  , 
a,  K=H='0  .  The  equation 
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*2 

(39.2)  y  =  _  ^  +  4  y 

2  y 

may  be  integrated  by  the  method  of  variation  of  parameters. 

One  finds 

y^  =  y  (4y^  • 

thus  ,  y{x)  is  an  elliptic  function  and  is  stable  , 

2  **  3 

One  may  also  set  y  =  and  determine  z  by  z  =  2  z"^  . 

b, K=0,H^0,A  linear  transformation  on  x  reduces 
(39.1)  to 

*2 

(39.3)  y  “  JL _ +  4  y^  +  2  y  . 

2  y 

This  equation  may  be  integrated  by  the  method  of  variation  of 
parameters  .  One  finds 

y^  =  4  y  (y^  +  y  +  K^)  ; 

thus  ,  y(x)  is  an  elliptic  function  and  is  stable  . 

One  may  also  set  y  ==  z  and  determine  z  by  z  *  2  +  z  . 

c,  K.H  ^  0  .  A  trivial  change  of  variable  brings  equation 
(39.1)  to 

*2 

(39.4)  y  =  J_  +  4  y^  +  2  X  y 

2  y 

2 

which  on  setting  y  =  z  becomes 
••  3 

z  =  2z-'  +  xz  , 
i,e.  equation  II  of  Table  I  . 


So. 


40,  Equation.  E.3«  Foi'  P2  ,  one  has  p  =  2  ,  b^s  =  4  •  By  a 

transformation  (  X,<p  )  ,  one  may  assume  =  0  ,  =  4  and 

-2  * 

s  =  1  .  Now  set  y=z  ,z=l+uz  and  determine  u  by 

z  u  =  2  u  -  ^  (bpZ  +  b.^z^)  +  ^  u^  , 

2  ^  2 

We  obtain  the  condition  for  stability  by  setting  u=  P  z+z  v, 

where 


2  p  =  b^  ,  P*  =  0  ; 

therefore  ,  b^  =  2  K  , 

For  Z2  ,  we  have  according  to  the  table  of  equivalence  , 

p  =  0  ,  b.-s  +  1  =  0.  Set  y  =  >  z  =  1  +  u  z  and  determine 

^  s 


u  by 

.u .  -  _L+f_L  -  5 

2s  \  2s  S  2 

s 

• 

The  condition  for  stability  is  s  =  0  or  b^  =  H  , 

The  stable  equation  E.3  is  thus  of  the  form 

*2 

(40.1)  y  =  _  +  3y^  +  2Ky  +  H  . 

4  y 


z+K 


uz  - 


z 

T 


u 


The  general  solution  of  (40,1)  may  be  determined  by  the  method 
of  variation  of  parameters  .  One  finds 

y^  =  4  [y3  +  K  y^  -  H  y  +  y3/2  j 

2 

or  on  setting  y  =  z  , 

*2_  4iv  2  ii.v 

z=z  +Kz-H  +  Kj^z; 
thus  ,  y(x)  is  an  elliptic  function  and  is  stable  . 


4l.  Equation  E.4.  By  a  transformation  )  ,  we  may 

assume  a^^  =  =  0  .  There  is  no  condition  for  stability  for  P2. 

For  Z2  ,  one  has  p  =  0  ,  b^s  +1=0;  the  condition  for  stabi- 
lity  is  s  =  0  or  b^  =  K  . 

The  stable  equation  £.4  is  of  the  form 

o  *2 

(41.1)  y  ^  +  K  . 

4  y 


If  K  =  0  ,  we  have  a  particular  case  of  (3B,1)  ;  y(x)  is  stable 
If  K  ^  0  ,  set  y  =  -  K  u  and  find  the  canonical  equation 

-  1  . 


3 

u  = 


*2 

u 


4  u 

The  general  solution  may  be  determined  by  the  method  of  varia- 

2  *'  *2 

tion  of  parameters  or  by  setting  u  =  z  .  One  finds  2zz  =  z  -  1 


and  upon  differentiation  z  —  0  ,  Therefore  ,  y(x)  is  stable  . 


42.  Equation  E.5.  By  a  transformation  (  4  >  )  ,  we  may 

assume  a^  =  b^  =  0  .  Then  ,  the  equation  E,5  is 


ayy  - 


na 


{n+2)' 


y^  +  1^1  y^  . 


The  only  condition  for  stability  is  produced  by  PI  ,  p  =  0  , 

nas  +  n  +  2  =  0  and  according  to  (37.4)  is 

b^  s  +  a  s  +  =-0 

ns 


or 


n 


n  +  2 


The  stable  equation  E.5  is  of  the  form 


&2, 


(42.1 


V  n  , 


ayy  - 


n  a 


y3  + 


n  I  2 

_  a  y  , 


n  +  2 


^  (n  +  2)*^ 

where  a  is  an  arbitrary  analytic  function  of  x  . 

To  integrate  (2+2,1)  ,  note  that 

(42.2)  ^  ^  y  =  n  V 

y  n  +  2 

is  regular  at  the  parametric  poles  .  Then  v  is  given  by 


v  +  v  =  Oorv  = 


X  +  H 


;  therefore  ,  y  is  determined  by  the 


Bernouilli  equation 

•  n 

y  = _ 

or  on  setting  y 
the  first  order 


y  + 


nay 


X  +  H 
=  1 
w 

w  + 


n  +  2 

,  by  the  linear  differential  equation  of 


nw 


n  a 


0  . 


X  +  H  n  +  2 

On  setting  (x  +  H)’^  w  =  W  ,  one  obtains 

W  +  (x  +  H)^  =  0  ; 

n  +  2 

thus  ,  y(x)  is  stable  . 

More  generally  ,  suppose  that  a^  and  bg  are  not  Identically 
zero  .  Then  ,  the  condition  of  stability  corresponding  to  PI, 
p  =  0  is 

’  n  +  2  , 

a  -  a^  a  -  _____  b^  =  0 

^  n 

and  gives  b^  . 

The  equation  E.5  becomes 


63, 


(42. 


1  1  \  *  •  n  3.  n  ^  2, 

3)  y  -~i  1  -  _ )  Ji _  =ayy+a.y-_  -  -n-  , 

n  ^  y 


y^+  ^  ( a-  a  )  y  ^-f  b  ^y 


n  +  2 


(n  +  2] 

To  integrate  this  equation  ,  determine  again  y  by  (42.2)  and  set 
w 

V  =  _____  ;  w  is  given  by 

w  •  b 

w  -  a^  w  -  ^2  w  =  0  , 

n 

and  y  is  a  solution  of  the  Bernouilli  equation 

*  w  .  n  a  2 

y  =  n  _  y  +  _  y  . 

w 


n  +  2 


n 


w 


On  setting  y  =  ”  ,  one  finds 

u 


u  + 


n  a 


n 


A*  M 

w  =  0 


n  +  2 

therefore  ,  u  and  y  are  stable  , 


43.  Equation  E.IO,  Both  PI  and  Z2  give  a  condition  for 

2 

stability  .  For  PI  ,  we  observe  that  the  condition  a  +  9b  =  0 
is  Identical  to  the  condition  given  for  equation  E,5  where 
n  =  4.  Then  ,  p  =  0  ,  2as  +  3  =  0  and  by  a  transformation  (  X  ,  cp  )  , 

3 

one  may  assume  a,  =  0  ,  a  =  -  and  s  =  1  .  The  condition  of 

^  2 

stability  is  b^  =  0  . 

For  Z2  ,  on  using  the  table  of  equivalence  ,  one  has  p  =  0, 

1  +  b^  s  =  0  ;  the  condition  for  stability  is  s  =  0  or  b^  =  K  /O. 
Thus  ,  the  stable  equation  E.IO  is  of  the  form 

(43.1)  y  -  ^  yy  -  ^  y^  +  b  y  +  K 

4  y  2  4 

where  is  an  arbitrary  analytic  function  of  x  . 


f 


To  integrate  (43.1) 


f 


64. 


observe  that  it  is  equivalent  to  the 


differential  system 


•  2 

y  +  y  =  2  V  y 


y  = 


K 


2  V  +  V*  «  b. 


On  eliminating  y  between  these  relations  (  consider  ^  )  ,  one 

y 

obtains 

v+3vv+v^-b^v-  ^  (b^  +  K)  =  0  , 

2  —  2 

This  is  an  eqviation  of  type  A(x,y)  =  0  and  class  PI  ,  p  =  0 

see  ^  23  t  20.4'!  ;  on  setting  v  =  ^  ,  it  reduces  to 

w 

the  linear  equation 

w*  -  b„  w  -  _ (b'_  +  K)  w  =  0  . 

2  —  2 

Therefore  ,  y  is  stable  , 

Remark  .  We  may  also  choose  (  R  ,  (p  )  such  that  =  b^  , 

a  =  -  ^  ,  I  Gambler  ,p.  28,eq.2^  ;  then  ,  one  finds  the 

2 


equation 


y  = 


3  *  y" 

^  y  y  -  -i. 


q 


+  (y  +  y*^)  +  ry  +  q  , 


where  q  and  r  are  arbitrary  analytic  functions  of  x  . 

This  equation  is  equivalent  to  the  differential  system 


y  +  y  =  2  V  y 


[  2  V 


+  V  - 


one  finds 


On  eliminating  y  between  these  relations  , 

V+3VV-  ^  ^  v+v^-  v^-  A  ^  ^  +r  ^  V  -  ^  (r-tq~  )=0. 

2  q  2q  ^  2q  ^2  y  2  q 

* 

On  setting  v  =  ^  ,  this  equation  of  type  A(x,y)  =  0  and  class 

w 


PI  ,  p  =  0  reduces  to  the  linear  equation 


44.  Equation  E.ll.  Suppose  n  ^  Z  and  consider  the  equation 
i  n  /  y 


(44.1) 


a  y  - 


na. 


(n  -  2)> 

For  P2  ,  there  is  no  condition  for  stability  . 
For  Z1  I  one  has 

n  a^  E  =>=  -  (n  -  2)  ,  p  =  0  ; 

na2S=(n-2)  (n-1)  ,  p* 


+  a^y  +  bg  y  +  b^  . 


n  -  2  . 


A  transformation  (  X ,  q)  )  may  be  chosen  such  that  a^  *  0  , 


n  -  2 


.  The  condition  for  stability  for  Z1  ,  p  =*  0  , 


(s  =  1)  is  then  b^  =  0  . 

Equation  E.ll  is  then 

(44.2)  y  -  (^1  - 


‘2 

y 


n  -  2  y 

y" 


+  by 


y  n  y  ny 

where  b  is  determined  by  the  condition  for  stability  for  Z1  > 
p=n-2,  s=l-n,To  obtain  this  condition  (  call  it 
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n*  2 

condition  A)  ,  set  u  *  P^_^{z)  +  2  *^11  and  write  n-2  linear 

equations  . 

For  instance  ,  suppose  n  =  4  ;  then 

*  z  2 

zu=2u  +  bz  -  _  u 

2  ' 

2  ' 

On  setting  u=pz+z  U,  one  finds  p+b=P,  p=0  and 
hence  the  condition  of  stability  b  =  K  . 

To  integrate  (44.2)  ,  observe  that  ^  ^  is  regular 

y  y 

for  the  parametric  zeros  Z1  ,  p  =  0  and  set 
y  =  1  +  n  V  y  ; 

then  ,  V  is  given  by  the  Riccati  equation 

;  t  v2  =  ^  b 
n 

w  ••  t) 

or  if  V  =  _____  by  w  =  _  w  .  Accordingly  ,  y  is  given  by 

w  n 

y  =  K  + 

If  w(x)  has  a  simple  parametric  zero  (  x  *  )  ,  then 

y(x)  =  (x  -  x^)  s(x)  ;  for  y(x)  to  be  stable  ,  b  must  satisfy 
the  condition  of  stability  A  . 

45.  Equation  E,12.  Consider  the  equation 

*2 

(45.1)  y  -  ^  ^  =  ^4  +  a-  y  +  b  y  +  b_  . 

2  y  y  ^  ^ 

For  P2  ,  there  is  no  condition  for  stability  . 

For  Z1  ,  one  has  according  to  the  table  of  equivalence  , 

2 

2b ^s  +  1=0  ,  p=0  ,  0, 


A  transformation  (  X,  <P  )  may  be  chosen  such  that  b„  =  0  , 


b,  =  -  _____  ;  then  the  condition  for  stability  for  Z1  is 
^  2 


a.  +  b  s  =  0  or  ,  because  s  =  +  1 


b^  =  0  . 


The  equation  E.12  is  then 

1 


(45.2] 


y  - 


_ 


2y 


and  may  be  integrated  by  the  method  of  variation  of  parameters, 
*2  1 

One  finds  y  =  uy  ,  u  =  +  2K  and  hence 


(45.3)  y^  =  1  +  2  K  y  . 

On  differentiating  (45«3)  >  one  obtains  y  *  K  and 

(45.4)  y(x)  =  ^  +  K^x  +  Kj,  ; 

2 


therefore  ,  y(x)  is  stable  . 

On  substituting  y(x)  given  by  (44.4)  into  (45*3)  »  one 
finds  =  1  +  2  K  . 


46.  Equation  E.13.  For  P2  ,  one  has  p=2  ,  b^s*3  ;  a 
transformation  (  X,qp  )  may  be  chosen  such  that  a^  =  0  ,  b^  *  3 
and  s  =  1  .  The  equation  giving  u  is  then 

zu  =  2u  -  ^  bpZ  +  a^z^  -  ^3  z^  +  ^2  z^  +  a-z^u  +  ^  u^ 
2  2  2  2 

2 

To  find  the  condition  for  stability  ,  set  u  =  p  z  +  z  v  and 
obtain 

2  p  =  bg  ,  P  =  a^  . 

•  • 

For  convenience  ,  set  b^  =  12  q  ,  a^  =  6  q  ;  then  ,  equation 


E,13  becomes 


(46.1)  y  - 


^ _  =  6  q  +  3  +  IE  q  y  +  b  - 

4  y  y  ^7 


For  Z1  ,  one  has  p=0,12qs+,l  =  0;p**2,  12qs  =  3, 
From  the  table  of  equivalence  ,  one  finds  that  n  =  -  4  and  that 
for  p  =  0  ,  the  condition  for  stability  is 

b„  s  -  6  q  s  +  ^  =  0  ; 

2s 


because  s  =  - 


,  one  has  b^  =  -  12  q  . 


12  q 

The  condition  for  stability  corresponding  to  Z1  ,  p  =  2  , 
determines  q  .  Indeed  ,  because  4  q  s  =  1  ,  one  deduces  from 
(37.4)  *  the  equation 


{46.2)  zu  =  2u  -  2  ^  +  Az  +  12  q  - 


z  u  - 


zu 


2q 


with 


A  =  - 


-L  _ 


+  12  q  . 


To  find  the  condition  for  stability  ,  set  u  =a  +  p  z  +  z'^v 
into  (46.2)  ;  one  obtains 


a  =  ^ 


P  =  2 


•2 

q 


-  12  q  ,  p  =  12q  -  2p  ^ 


so  that  ,  on  eliminating  p  , 

q^^^  =  12  q^  +  12  q  q 


Therefore  ,  the  stable  equation  of  this  class  is  of  the  form 


(46.4)  y  =  2.JL 
4  -  y 


+  6  q 


+  3  y^  +  12  q  y  -  12  q  - 


where  q  is  a  solution  of  (46,3)  • 

To  Integrate  (46,4)  >  we  use  a  method  which  gives  also  the 

•  • 

values  of  q  .  Observe  that  y  +  12  q  is  regular  for  the  parame¬ 
tric  zeros  corresponding  to  p  =  0  .  Therefore  ,  the  equation 

(46.4)  is  equivalent  to  the  differential  system 

•  « 

y+12q  =  -2vy  , 

3y  =  “2v-fv^-  12  q 

which  shows  that  v  is  stable  if  y  is  stable  and  conversely  . 

On  eliminating  y  between  these  relations  ,  it  is  readily  seen 
that  V  is  determined  by  the  equation 

v  =  -vv  +  v^-  12qv+12q 

which  is  of  type  A(x,y)  s  0  and  class  PI  ,  p  =  2  and  5  Lsee  (25«5)1 
For  this  equation  to  be  stable  ,  it  is  necessary  that  q 
satisfies  (46.3)  »  the  general  solution  of  the  corresponding 
equation  is  then  given  by 


V 


» 


where  V  is  a  solution  of  (46.3)  distinct  from  q  .  Therefore  , 
y  is  stable  , 

According  to  the  values  of  K  and  H  ,  one  finds  one  or  the 
other  of  the  three  following  canonical  equations  , 

a.  K=H=0  .  Then  (46.3)  is  q  =  6  q^  or  q^  =  4  q^  +  . 
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When  =  0  ,  one  has  q  =  (x  +  K^)  ;  a  trivial  change 

of  variable  (  x  +  -•  x  )  brings  equation  (if6,4)  to 

(46.5)  5  -  J_  il  +  3y^  -  ^  i_  -I-  ^  -  ''*'*  . 

4  y  X  y  X  ^4  “6“ 

When  ^  0  ,  one  has  q  =;|3(x;0,K^)  where  ^  is  the  well  known 
elliptic  function  of  Weierstrass  ;  equation  (46,4)  is  then 


(46.6)  y 

b.  K  = 
therefore 


=  —  ^  +3y^+6  ^  JL  +12^  y-12  £  -  . 

4  y  '■  y  y 

0  ,  H  ^  0  .  Then  q  =  6q^  +  H  ,  q^  =  4  q^  +  2Hq  +  ; 

,  q  *=  ^  {x;-2H,H^)  and  (46,4)  has  the  same  form  as 


(46.6)  . 


c .  K  ^  0  ,  On  using  a  trivial  change  of  variable  ,  one  may 

** 

assume  K  =  1  ,  H  =  0  so  as  to  reduce  (46.3)  to  q  =  6q  +  x  . 


47.  Equation  E,14.  For  P2  ,  one  has  p  =  3  ,  b^s  =  4  and  for 

p 

Z1  ,  p  =  0  ,  2b^s  +  1  =  0  .  By  a  transformation  (  X  ,  cp  )  ,  one 

may  assume  b.,  =  4  and  b,  =  -  ^  , 

^  2 

First  ,  we  consider  Zl  ,  The  condition  for  stability  is 

a^  +  b^  s  =  0  ,  s  =  +  1  so  that  a^  *  0  ,  b^  ==  0  . 

For  P2  ,  one  has  s  =  1  and 

zu  =  3  u  -  _______  b-  z  +  _____  z-^  +  z  u  . 

2  '4 

To  obtain  the  condition  for  stability  ,  set  u  =»  pz  +  y  +  z^v 
and  find 

•  • 

4p=b2,  P=Y  »  Y=0 
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so  that  . bg  =  K  X  +  H  . 

The  stable  equation  E.14  is  of  the  form 

*2 

(47.1)  y  =  _Z _  +  4  y^  +  (K  X  +  H)  y  -  _  . 

2  y  .  2y 

According  to  the  values  of  K  ,  we  have  the  following  canonical 
equations  . 

a.  K  —  0  .  One  has 

*2 

(47.2)  y  =  ^  -  _ L  +  4  y^  +  2  H  y  ; 

2  y  2y 

the  method  of  variation  of  parcimeters  gives 

y^  =  4  y^  +  4  H  y^  +  2  H^y  +  1  . 

Therefore  ,  y(x)  is  an  elliptic  function  and  is  stable  . 

b,  K  ^  0  ,  A  trivial  change  of  variables  (Kx+H  -►  x,K^y  y  ) 

reduces  (47.1)  to 

*2 

(47.3)  y  =  _  -  J. _  +  4Hy^-xy  ,  H^O. 

2y  2y 

To  integrate  (47.3)  ,  observe  that  ^  ^  =  2  v  is 

y  y 

regular  for  the  parametric  zeros  corresponding  to  s  =  -1  ,  The 
equation  (47.3)  is  equivalent  to  the  differential  system 

y  =  2.V  y  -  1  , 

.  2 

4Hy  =  2v  +  2v'^+x  . 

V 

On  eliminating  y  between  the  relations  ,  one  finds 

v  =  2v^+xv-2H-  ^ 

2 

so  that  V  is  a  solution  of  a  Painleve  equation  II  (  see  Table  I). 
Therefore  ,  y(x)  is  stable  , 


92 


46.  Equation  E.I5.  Consider  the  equation 

*2.2  * 

(46.1)  y  -  ^ _  =  ayy-  _  y^+  ^4..  +a-y+b^y^+b  y+b  . 

2  y  y  I  I  <  J 

For  PI  ,  one  has  p=0,as+2=0;  there  is  no  condition 

for  stability  for  PI  ,  p  =  -  4  . 

2 

For  Z1  f  one  has  2b^s  +l=0,p=0,0. 

A  transformation  (  X  >  (p  )  may  be  chosen  such  that  a  =  -  2  , 

b,  =  -  ^  ,  The  conditions  for  stability  are 

4 

for  PI  ,  a^  =  b^  ; 

for  Zl  ,  a^  +  b^  s  =  0  ,  s  =+  1  and  hence  a^  -  b^  =  0  , 
Therefore  ,  the  equation  (46.1)  oen  be  written  in  the  form 

*2  .3 

(46.2)  y  -  J, _  «  -  2  y  y  -  .  +  f  y  , 

2  y  2  2  y 

where  f  is  an  arbitrary  analytic  function  of  x  . 

To  integrate  equation  (46.2)  ,  observe  that  ^  +  y  is 

y 

regular  for  the  parametric  poles  PI  ,  p  =*  0  and  that  ^  ~  ^  is 

y 

regular  for  the  parametric  zeros  Zl  .  Accordingly  ,  set 

(46.3)  y=-y^+l+2vy  ; 

substitution  in  (46.2)  gives 

(46.4)  V  +  v^  =  (f  -  1) 

2 

and  y(x)  is  determined  by  the  two  Riccati  equations  (46.3-4). 

To  prove  that  y(x)  is  stable  ,  set  y  *  ^  ;  then  from 

w 

(46.3)  ,  one  obtains  2  v  w  =  w  -  w  and  (48,4)  becomes 

•  2  *? 

2  w  w  =  w*^  -  w*^  +  2  f  w*^ 
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whence  ,  upon  differentiation  , 

(it)  ••  *  * 

w'^^  -2fw-fw+w=0  . 

Thus  ,  w(x)  is  stable  and  also  y(x)  , 

More  generally  ,  if  a  ^  -  2  and  ^^4.  ^  “ 


1 

T 


one  finds 


the  conditions  for  sts.billty 

for  PI  : 


(46.5: 


for  Z1 


2  b^  =  a  -  a  a^  ; 

b-  =  0  and  b,  =  2  a.  b, 
3  414 


On  setting  b^  =  r  ,  the  stable  equation  E,15  becomes 


(46.6) 


y  “ 


=  ayy  + 


1 

T 

r 

FT 


y 

a^ 

y  - 


If  a  - 

2  \  ■  2  r  / 


+  b^y  + 


r 

T 


49. 


(49.1) 


Equation  E.14.  Suppose  n  ^  2  and  consider  the  equation 


V  - 


n  a 


y3  - 


(n+2)' 


^  ^2 _  +  a^  y  +  b^  y^  +  bg  y  +  b^  . 

(n  >  2)^  y 


Note  that  E,14  reduces  to  E.ll  when  a  *  b^  =  0  .  One  has 
for  PI  : 

n  a  s  =  -  (n  +  2)  ,  p  =  0  ; 

n  a  s  =  -  (n  +  2)  (n  -  1)  ,  p  =  -  (n  +  2)  ; 

for  Z1  : 

na2S=-(n-2)  ,  p=0  ; 

nagS=(n-2)  (n-1)  ,  p=n-2. 
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n  "f"  2 

A  transformation  (  \  ,  9  )  may  be  chosen  such  that  a  =  -  ** 

n  -  2 


n 


and  a^  -  - 


.  Then  ,  equation  (49* 1)  becomes 


’-(•■4-) 


_ 


-  +  2  y;  .  n  -  2 


n 


n 


7 

T" 


(49.2) 


y- 

n 


ny 


+  a^y  +  ^3 


The 

conditions 

for  stability  are 

for 

PI  , 

P  = 

0 

f  s  "  1  •  a  "  b 

for 

Z1  , 

p  = 

0 

,  s  =  1  :  a^  +  b^ 

for 

PI  , 

P  = 

- 

(n+2)  p  there  is  no 

set 

==  b^ 

«s 

-  b^  =  b  and  b^  = 

Now- 

can  be  re-written 


(49.3) 


'■(■■4-) 


. 


n  +  2  •  n  -  2 

_ yy  -  _ 


n 


n 


y 

T" 


y- 

n 


■tb(  y  +  y'^-  1  )  +  fy  . 


ny 


It  remains  to  consider  the  condition  for  stability  for  Z1  , 
p=n-2,s=l-n;  this  condition  (  call  it  condition  A  ) 
gives  a  relation  between  b  and  f  and  their  derivatives  .  To 
find  this  condition  A  ,  set  as  usual 
u  =  Pj^.3(z)  +  U 

and  apply  our  basic  theorem  IV, 

For  instance  ,  consider  n  =  2  and  n=  3  . 

For  n  =  2  ,  one  has 
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y 


2 _ ]L  =  -  2yy  -  _  +  b(y  +  -  1)  +  fy 

2  y  2  2y 


and  condition  A  is  b  =  0  ;  one  finds  equation  E.15  (see  48»2) , 
For  n  =  3  ,  one  has 


y 


^  yy  “  ^  ^  ^  ^  +b(y+y^-l)+fy 

3  3  y  3  3  y 


and 


zu  =  u  +  3b  +  /  __+f  j  z  - 


bz‘ 


^  +bzu  +  ^  z^u 

12  “S" 


the  condition  for  stability  is 


(49.4) 


3  b  +  f  + 


-T 


-6  b* 


and  determines  f  . 

Suppose  that  the  condition  A  for  stability  has  been  determined. 

To  integrate  (49.3)  »  set 

•  2 

y  =  -y  +  l+  nvy 

so  that 

V  +  =  b  V  +  ^  ^  ~ 

n  \  n 

•  • 

u.  w 

Now  ,  set  y  =  _  ,  v  =  _  ;  w  and  u  are  given  by 

u  w 

(49.5)  w  -  b  w  -  J_(f  - 

n  V,  n 

(49.6)  u  -  n  ^  u  -  u  =  0  . 

w 

The  singular  points  of  w  are  fixed  .  Let  x  =  be  a  simple 
parametric  zero  of  w(x)  ;  then  x  =  x^  is  a  regular  singular 
point  (  in  the  sense  of  Fuchs  )  of  (49.b)  .  Note  that  the  roots 
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of  the  corresponding  indlcial  equation  are  0  and  n+1  ;  therefore  , 
u(x)  may  have  a  logarithmic  term  except  if  a  certain  condition 
(  call  it  condition  B  )  is  satisfied  .  In  this  case  ,  equa¬ 
tion  (49.6)  has  two  regular  solutions 

u^(x)  =  (x  -  £  1  +  . . ,  J  , 

/  \  1  .  (x  -  X  ) ^  , 

u^(x)  =  1  +  _ o  +  ... 

^  2(1  -  n) 

to  which  correspond  respectively 

y  1  ( X )  _________  +  . . .  > 

X  -  X 

o 

y^vX)  ~  o  +  ...  . 

1  -  n 

Therefore  ,  conditions  A  and  B  coincide  and  y(x)  is  stable  . 


50.  We  now  recall  the  method  used  by  B.  Gambler  to  deter¬ 
mine  the  stable  equations  E.14.  Set 


(50. i: 


y  =  -  y  - 


+  n  V  y 


n  -  1 

and  substitute  into  equation  (49.3)  ;  then 


(50.2)  yl  V  +  v^-bv- 


n 


n^-2n+2 


n'^  ( n  -  1 ) 


■] 


n  -  1 


n  -  2' 
n  -  1 


V. 


To  abbreviate  ,  set 
n  b 


oc  =  - 


p  =  - 


n  -  2 

M 

b 


V  =  Z  +  g 
n 


+  ”  ^  b^  -  ^ 


"  "  ^  (n  -  2)^ 


n 


n"-  -  2n  +  2 


n*^  ( n  -  1 ! 
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P(z)  ~  +  a  z  +  P 

so  that  (50.1)  and  (50.2)  become  respectively 

(50.3)  y  =  -  y^  -  ^  +  a  y  +  n  z  y  , 

n  -  1 

(50.4)  yfz  +  P(z)l  =  -  ^  z  . 

n  -  1 

Hence  ,  if  y(x)  is  stable  ,  v(x)  is  also  stable  and  vice 


versa 


On  eliminating  y  between  (50,3)  and  (50.4)  >  one  finds 


n ( n- 1 ) 


zz  - 


n  -  2 


n  a 


n  -  2 


z  + 


n  -  4 
n  -  2 


(50.5) 


(n-1)^  ,3  n(n-l)a  2  ^  „  p  ^ 


Z“^  - 


n  -  2 


n  -  2 


n  -  2 


where 


n_  *.n-l  2,n  -2n  +  2„  n-2 

=  a  +  _  a  +  _ p  -  _  , 


n-2 


n-2 


n  -  1 


Bq  =  P  +  «P  . 

3  n  -  2 

Equation  (50.5)  is  an  equation  of  the  type  E,14  with  n  repla¬ 
ced  by  n-2.  Therefore  ,  we  have  to  consider  two  cases  according 
as  to  whether  n  is  even  or  odd  , 

a)  n  =  2k  .  Suppose  n  =  4  .  The  equation  (50.5)  becomes 

•  .  .  ^  2 

(50.6)  z  »  ^  ^  -6zz-2  az-  ^  z^-6 a  z^-BpZ-(  p  +2  ap  )-  P 

2  z  2  2z 

If  P  =  0  ,  this  equation  is  stable  and  of  the  type  E.5 


for  n  =  2  L  see  42,3^  . 


If  p  0  ,  equation  (50,6)  is  an  equation  of  the  type  E.15 
IJsee  i+6,6j  and  is  stable  if  the  condition  for  stability 
P  +  2  «P  =*  0  is  satisfied  . 

For  n  =  2k  ^  1+  ,  repeated  applications  of  the  method  give 
equations  of  types  E,5  or  E.15  and  hence  k  distinct  classes  of 
stable  equations  . 

b)  n  =  2k  +  1.  Suppose  n  =  3  .  One  has  p  =  0  j^see  (49.4)J  . 

For  n  =  2k  +  1  ^  3  and  P  =  0  ,  equation  (50.5)  is  an  equation 
of  the  type  E. 11  (see  42,3)  and  is  stable  ,  If  p  ^  0  ,  repeated 
applications  of  the  method  give  again  equations  of  the  type 
E.ll  and  hence  k  distinct  classes  of  stable  equations  , 

51.  Equations  E,6  and  E.17.  Consider  the  equation 

*2  .  2 

(51.1)  y  -  ^  ^  =  ayy+a  y+  ^  y^+b^y^+b^y+b  +  ^4  . 

2y  ^  2  ^  '"^y 

If  b^  0  ,  we  have  equation  E,17  and  if  b^  ==  b^  ==  0  ,  equation 

E.6. 

Assume  b^  7^  0  .  Then  for  Z1  ,  one  has  p  =  0  ,  2b^s^  +  1  =  0. 
The  condition  for  stability  is  +  b^s  +  ^  =  0  so  that 

3 

•  • 

b„  =  0  ,  2a.  =  4  .  For  convenience  ,  set  a.  =  q  =  ^  and 

2 

b^  =  -  72  H  r  f  where  H  is  a  constant  .  Note  that  H  =  0  gives 
equation  E.6, 

Now  ,  consider  PI  ;  one  has  p=  1  ,  as=-l  and  p  =  5  » 
as  =  3  .  By  a  transformation  (  4  ,  cp  )  ,  one  may  assume  a  =  -  1 
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or  s  *  1  when  p  =  1  ;  then 

•  3  z  2 

2U  =  u  +  q-bj^-  bgZ  -  b^z"^  +  q  z  u  +  _____  u 

2 

To  obtain  the  condition  for  stability  ,  set  u  =  a  +  z  v 
and  find 


a  =  b^  -  q 


a  =  -  bp  +  a  q  +  __2. 
^  2 


therefore  , 

(51.2)  b. 


1  /,  2  ^2^  C  « 

••  CJ  )  ■•  D  2^  "i*  Q  • 


For  Pl,p=5>as-3>  one  has 


The  condition  for  stability  determined  by  this  equation 
is  complicated  ;  we  use  another  method  . 

A  transformation  ( X  ,  t  )  may  always  be  chosen  such  that 
b^  +  2a^  =  0  ;  then  ,  b^  =  -  2q  and  b^  =  3^  q  + 

from  (51. 2)"^  .  Therefore  ,  (51.1)  may  be  rewritten  as 

(51.3)  y  >  il  -y;+q;+J_  y3-2qy^3/q  +  j1)  y-  . 

2  y  2  \  2  y  y 

Now  ,  set 

(51.4)  12w=y+y^-3qy 

and  note  that  w  is  regular  for  the  parametric  poles  correspon¬ 
ding  to  p  =  1  and  has  a  double  pole  at  the  parametric  poles 
corresponding  to  p  =  5  . 

Equation  (51.3)  gives 
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(51.5) 


w  -  q  w 

so  that  (51.3)  is  equivalent  to  (51.4-5)  5  moreover  ,  w  is 
stable  if  y  is  stable  and  conversely  . 

On  eliminating  y  between  (51.4-5)  >  one  obtains 

w=6w^+(q+q^  )w-6Hr^ 

.  2*22 
observe  that  r=qr,qw  -Hrr=  q(w  -  Hr  ) 

*  2  r 

Because  q  +  q  =  _  ,  one  finds 

r 


w  =  6  ^  w  -  6  H  r^  . 

r 

This  is  an  equation  of  type  A(x,y)  *  0  and  class  P2  ,  p  *  5. 

To  obtain  the  canonical  form  of  this  equation  ,  set  12  g  =  ^ 

r 


and  w  =  V  -  g  so  that 

(51.6)  v=6v^+g-6g^-6Hr^. 

We  have  to  consider  two  cases  according  as  to  whether 
H  »  0  or  H  51^  0  . 

i,  H  =  0  ;  we  obtain  equation  E.6  .  In  order  that  (51.6)  be 
stable  ,  g  and  v  must  be  solutions  of 

(51.7)  V  =  6  v^  +  K  X  +  Kj  . 

If  Vj^  and  v  are  respectively  a  particular  and  the  general 
solution  of  (51.7)  f  one  has  g  =  v^  ,  r  =  12  v^  r  and 


therefore  ,  y(x)  is  stable  , 
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ii,  H  ^  0  ;  we  obtain  equation  E.17  .  We  may  evidently  assume 
H  =  1  .In  order  that  (51.6)  be  stable  ,  one  must  have 


(51.61 


g=6e^+6r^+Kx+Kj 
V  =  6  v^  +  K  X  +  . 


Because  r  =  12  g  r  ,  one  may  write 

.2 


g  +  v  =  6(g  +  r)'^  +  Kx  +  K^  ; 


therefore  ,  if  v  ,  v^  ,  v^  are  distinct  integrals  of  (51.7)  > 
one  has  g  +  r  =  Vj^  ,  g-r  =  Vg  and 


2  g  =  Vi  +  Vg  . 


2r  =  Vi  -  V2 


"^1  ■  ^2 


Moreover  , 


w  =  V  -  , 


(v  -  Vj) (v  -  Vg) 

- ^ - ^ - 6— 


’  v,  + 
V  -  1 


la  -< 


V.  +  v^ 
V  -  1  2 


therefore  ,  y(x)  is  stable  . 


52.  Equations  E.7  and  E.16.  Consider  the  equation 


*2  . 

(52.1)  y  -  ^  ^  =  a,y  +  by^  +  b^y^  +  by  +  b  +  ; 

2  y  ^  ^  y 

if  b^  7^  0  ,  we  have  equation  E.I6  and  if  b^  =  b^  =  0  ,  equation 

E.7. 

By  a  transformation  (  X  ,  cp  )  ,  we  may  assume  a^  =  0  ,  b  = 

2  •  2 
For  Z1  ,  one  has  p  ■■=  0  ,  b^s  +1  =  0.  The  condition  for 
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stability  isb^  +  s==0  ,  i.e,  b^^O  ,  s  =  0orb^=K?^0, 

For  PI  ,  one  has  p=2  ,  s=+l  and 

•  2  3  z  2 

zu  =  2u  -  b.s  -  b^z  -b„sz  -b,z+  u  , 
123  4 

2 

To  obtain  the  condition  for  stability  ,  set  u=o  +Pz+zv 
and  find 


2a-  b^s  =0  , 


a  =  p  -  b  +  “ 


2 


p  =  0 


so  that  2b^s  =  -  /fb^+b^b^  ,i.e.  b^  =  0  ,  i+b^  =  b^b^  or 


b^  =  X  +  H  , 


d  bg  =  +  e  Kg  . 


We  have  to  consider  two  cases  according  as  to  whether 
=  0  or  7^  0  . 

a.  =  0  ;  then  '  ^4  constants  and  one  finds 

the  equation 

*2 


(52.2)  y  =  Jl_ 
2y 


y^  +  4Kj  y^  +  2K  y  - 
2  2y 


Using  the  method  of  variation  of  parameters  ,  one  obtains 

y^  =  y^  +  4  K^  y^  +  4  Kg  y  +  K  +  y  ; 
y(x)  is  an  elliptic  function  and  is  stable  . 

{[  One  may  also  multiply  (52.2)  by  2  ^  ,  integrate  and 

take  into  account  2  ^  ^  ^  =  d 


dx 


-2. 

y 


b,  K^  ^  0.  Then  ,  by  a  trivial  change  of  variable  ,  one  may 

2 

assume  b^  =  k.x  and  bg  =  2x  +  Kg  ;  the  equation  (50.1)  becomes 

(52.3)  y  =  JL  +  -L  +  ^xy^  +  (2x^  +  K  )y  +  ^ 

2  y  2  2y 
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which  is  the  equation  iv  of  Table  I  . 

For  K  “  0  ,  one  obtains  the  corresponding  stable  equations 

E.7. 


53.  Equations  E.6  and  E,19.  Consider  the  equation 

*2 


(53.1)  y  - 


^  ^  ay+a^+  ^2  ^  y+  ^ 


a^y^+b^y^+b^y+b^-  ^ 


if  0  ,  v/e  have  equation  E.19  and  if  *=  b^  =  0  ,  equation 

E.e. 

p 

By  a  transformation  (  X  ,  <[*  )  ,  we  may  assume  a  =  «  _____  . 

3 

Then  ,  one  has 

for  Pl;p=l,s=l  and  p=3,s=-2; 
for  Z1  ;  p  =  0  ,  3^2®  *  -  1  and  p  =  1  ,  3agS  =  2  . 

The  condition  for  stability  for  Z1  ,  p  =  0  is 

• 

•  2  s 

a.  +  b„s  -  a^,  s  +  ____  _____  =  0 
^  3  s 

or 

bs  .  3  aj  .  3  . 

For  convenience  ,  let  us  write  a^  =  r  ,  3  a^  =  2  q  so  that 
(53.2)  b^=2qr-3r  . 

For  PI  ,  p  =  1  ,  one  has 

2 

•  P  '3 

zu  =  \x+a^-h^+z{&^-h^)-h^z  \x+  ______  . 


The  condition  for  stability  is 


aj  -  bj  +  ag  .  bg  -  aj(aj-bj) 


0 
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To  simplify  the  notations  ,  we  may  assume  +  5  =  0  , 

on  using  a  transformation  (  X , f  )  .  Then  3b^  =  -  10  q  and 

(53.3)  b  =  4  q  +  r  +  . 

•  3 

Equation  (53«1)  becomes 

(53.4)  y  -  ^  ^  ^  yy+  ^  qy+r  ^  ^  y^  -  qy^ 

3  y  3  3  y  3  3 

3 

+  b  y  +  b  -  ^  , 

^  ^  y 

where  b^  is  given  by  (53.3)  and  b^  by  (53.2)  , 

The  conditions  for  stability  for  Z1  ,  p  =  1  and  PI  ,  p  =  3 
are  complicated  ;  we  use  another  method  • 

Observe  that  ^  +  y  -  4q  is  regular  for  the  parametric 

y 

poles  PI  ,  p  =  0  and  that  ^  is  regular  for  the  zeros 

y  y 

Z1  ,  p=0  .  On  setting 

(53.4)  ^  +  y  -  4  q  +  =  V 

y  y 

one  obtains 

,  ^  •  v^ 

(53.5)  V  =  -  _  -  2  V  q  +  V  y  , 

3 

Equation  (53.4)  and  the  differential  system  (53.4“5)  are 

thus  equivalent  ,  On  eliminating  y  between  (53.4-5)  >  one  finds 

'•212  2  * 

V  ®  _____  v^  +  2v  q  +  4vq  -  2vq  +  3rv  ; 

9 

this  is  an  equation  of  type  A(x,y)  =  0  which  ,  on  setting 
V  =  3(w  -  q)  ,  reduces  to  the  canonical  form 


For  this  equation  to  be  stable  ,  it  is  neces- 


Q  see  (24.3)]  . 
sary  that 


3  r  -  2  q  -  2  q^=  K  x  +  H  , 
q  +  2qq-3rq  =  K^  . 

On  eliminating  r  between  these  relations  ,  one  obtains 
(53.7)  q  =  2  q3  +  (Kx  +  H)  q  + 

equation  II  of  Table  I  ;  then  , 

(53. S)  3r-Kx+H+2q+2q^  . 

The  stable  equations  for  this  class  are  of  the  form 


L  yy  +  J-  qy  +  ^  JL  +  ^ 

3  3  y  3 


(53.9) 

-  qy^  +  (4q+r+  ^  q^ )y+( 2qr-3r )  -  ^  ^ 

3  3  y 

where  q  is  given  by  (53.7)  a-nd  r  by  (53.6)  . 

The  general  Integral  of  (53.9)  is  given  by 


*  *  2  2 

w  -  q  +  w  -  q 

y  =  ^ » 

w  -  q 

where  q  and  w  (  ^  q  )  are  respectively  a  particular  solution 
and  the  general  solution  of  (53.7)  . 

To  obtain  equation  E.8  ,  one  has  to  assume  r  =  0.  Equa¬ 
tion  (53.6)  is  then 

w=  2w^  -  2(q+  q^)  w+  q+  2qq  ; 


for  this  equation  to  be  stable  ,  it  is  necessary  that 
(53.10)  q+q^=Kx+H 
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so  that  q  +  2qq  =  K  . 

The  stable  equation  E.6  is  then 


*2 

2  y 

y  =  _ 

3  y 


2  *  ,  2  *,2  3  10  2  ^ 
_ yy+ _  qy+ _ y  -  _  qy  + 


(53.11) 


where  q  is  a  solution  of  (53.10)  .  Its  general  integral  is 

*  *  2  2 

^^w-q  +  w‘'-q‘=' 

w  -  q 

where  w  is  a  solution  of 

w  =  2  w^  -  2(Kx  +  H)  w  +  K 


54.  Equations  E,9  and  E,20.  Consider  the  equation 

*2 


(54.1)  y  -  ^  ^  =^ay+aj^+  ^2  ^  y+5a^y^+b^y^+bgy+b2-  ^  ^ 


5  y 


9y 


if  a  ^  0  ,  we  have  equation  E,20  and  if  a_  =  b_  =  0  ,  equation 


E.9. 


By  a  transformation  (  X  ,  tp  )  ,  we  may  assixme  a  =  - 


Then  ,  one  has 

for  PI  :  p  =  1  ,  s  =  1  and  p  =  2  ,  s  =  -  ^  ; 

2 

for  Z1  :  p  =  0  ,  ^a^s  =  -  3  and  p  =  3  ,  =  12, 

We  consider  the  simplest  conditions  of  stability  . 
For  Z1  ,  p  =  0  ,  one  finds 


^2  ^1^2  ■ 


^3  * 


ro  ro 


/ 
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For  PI  ,  p  =  1  ,  one  has 

•  2  3  z  Z 

zu  =  u+aj^-b^+(a2-b2)  z-b^z  -b^z^+Ca^+a^z )  zu+  _____  u 

On  setting  u  =  a  +  z  v  ,  one  obtains 

a  4-  a,  -  b.  =  0  , 

.  ^  ^  2 
«  =  a  +  «  . 

5 

For  convenience  ,  on  using  a  transformation  (  \  ,ip  )  ,  we 


assume  b^^  +  14  a^^  =  0  and  write 


5  a^  =  -  q 


,  5  b^  -  14  q  , 


^2  =  ^  • 


Hence  ,  the  conditions  for  stability  are 
for  Z1  ,  p  =  0  :  3  b^  =  -  (qr  +  5r  )  ; 
for  PI  ,  p  =  1  :  bg  =  r  -  3  q  +  . 


Then  ,  equation  (54«1)  is 

*? 

••  4  y  2  *  1  *  y 

y  =  _  ^  _  yy-  _  qy+r  _i 

5  y  5  5 


_ +  y 

y  5 


^  y3+  14q 


(54.2) 


The  conditions  for  stability  for  PI  ,  p  =  2  and  Z1  ,  p  =  3 
give  the  values  of  q  and  r  ;  because  these  conditions  are 
complicated  ,  we  use  another  method  . 

Set 

JL  +  y  +  3q  +  ^ _ L.  =  -  V  . 

y  3  y  2 

Note  that  v  is  regular  at  the  parametric  poles  PI  ,  p  =  1 
and  at  the  parametric  zeros  Z1  ,  p  =  0  and  has  simple  poles 


r-3q+ 


6 

T 


y  - 


1 

T 


(qr  +  5r)  -  ^  ^ 


loe 


for  PI  ,  p  =  3  and  Z1  ,  p  =  2  .  The  equation  (54.2)  is  equiva¬ 

lent  to  the  differential  system 


5  2  5 

vy-y  -  3qy-  r 


y  = 


V 


-  q 


V  2 

On  eliminating  y  between  these  relations  ,  one  finds 

,3 


V  =  -  (v+q)v  +  +  3  q 


q  +  2  q^  -  r  )  V  . 


This  equation  is  of  type  A(x,y)  =  0  and  on  setting  v  +  q  =  w 
reduces  to  the  canonical  form 


w 


-  ww  +  +  (2q  -  q^  -  ^  r)w  +  q  -  qq  +  ^  r  q  , 

3  3 


For  this  equation  to  be  stable  ,  it  is  necessary  that 


(54.3) 


2  q  -  q^  -  r  =  -  12  V. 


q  -  qq  +  r  q  *  12  v  ^  , 


where  is  a  solution  of 


(54.4) 


V  =  6v‘^  +  Kx  +  H. 


From  (54.3)  >  one  deduces 
(54.5)  q  =  -  qq  +  q^  -  12  v^^  q  +  12 
The  general  solution  of  (54.5)  is  ("see  (25.5)1 


V  -  V, 


where  V(  /  v^)  is  the  general  solution  of  (54.4)  ;  then  , 


(54.6) 


r  =  24  +  12  V  -  3 


/ 

i 


IJLij 

V  -  V, 


109 . 


The  general  solution  of  (62)  where  q  and  r  are  given  by  (54. 5-^6)  , 
is 

•  « 

^  _  w  -  q  w  q 

w  -  q  2  2 

where  w  q)  is  the  general  solution  of  (54.5)  • 

To  obtain  equation  E.9  >  one  has  to  assume  r  =  0  .  Then  , 
instead  of  (54.3)  >  one  finds 

2  q  -  q*^  =  -  12  v^  ,  q  -  qq  =  12  v^ 

so  that  v^  =  K  is  a  constant  ,  The  general  solution  of  the ■ 
corresponding  equation  E,9  is 


w  -  q  2  2 


where 

w  =  -  w  w  +  w3  -  12  K  w 
2  q  =  q^  +  K  . 
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VIII.  Equations  of  the  type  A(x,y) 


55 

(55.1) 


.  The  stable  equations  of  this  type  are  of  the  form 


/  1  +  1  ^  »  y)  V  +  N(x  ,  y) 

\  2y  y  -  1  /  y(y-i)  y(y  -  l) 


* 


where 


M  ^  ,  M 

_  =  a  y  +  +  _ o_ 

D  D 


^  =■  b  y^  +  b.  y^  +  bj,  y  +  b  + 

D  ^  ^  ^  B 


M^(x,y)  ,  N^(x,y)  are  polynomials  of  the  first  degree  in  y  . 

The  values  of  y  for  which  the  general  existence  theorem  of 
Cauchy  does  not  apply  are  y  =  oo  ,  0  and  1  . 

If  a  and  b  are  not  both  zero  ,  y(x)  has  simple  parametric 
poles  ;  if  a  =  b  =  0  ,  y(x)  has  double  parametric  poles  . 

Q  • 

If  y{x)  has  a  simple  parametric  pole  ,  set  y  =  »  z  =“  1+uz 

z 


and  determine  s  by  (note  that  n  =  -  2) 

2bs^-2as-l  =  0  ; 

because  p  =*  as  and  must  be  an  integer  ,  one  has  a  =  0  ,  b  0  . 

If  y(x)  has  a  double  parametric  pole  ,  one  has  a  =  b  =  0  ; 
because  p  =  -  1  ,  one  also  has  b^  =  0  , 

To  investigate  the  coefficients  of  M  and  N  in  a  neighborhood 
of  a  parametric  zero  ,  replace  y  by  w“^  in  (55.1)  so  as  to  obtain 
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This  equation  has  the  same  form  as  (55.1).  Therefore'  ,  if  y{x) 
has  simple  parametric  zeros  ,  then  M^(x,0)  =  0  ,  N^(x,0)  ^  0  ; 
if  y(x)  has  double  parametric  zeros  ,  then  M^(x,0)  =  0  ,  N^(x,0)=0 
and  the  coefficient  of  y  in  N(x,y)  is  also  zero  . 

To  sum  up  ,  and  changing  the  notations  ,  equation  (55.1) 
may  be  rewritten  as 

(55.2)  y  ^  ^  y^  +  ~  ^  y  +  y(y-l)  Q(x,y)  , 

\2y  y-l/  y-1 

where 

(55.3)  Q(x>y)  ■■  e  y  +  f  +  _f _  +  Jl  +  **  +  _  • 

y  (y.l)2  -  1 

Note  that  e  =  0  implies  f  =  0  and  that  g  =  0  implies  h  =  0  . 

The  transformation  y  =  w  -  1  brings  equation  (55.2)  to  an 
equation  of  the  same  form  according  to  the  equivalence  table 
(y)  ;ab  efghk£ 

(w):  b  a  g  h  e  f  kk-£ 

The  substitution  t  =  q)  (x)  transforms  (55.2)  into 


M  +  1  ') 

1  v*^+  r  (a-  $ )y-b+4  1 

1  1  +)'<y-l>  O(x.y) 

\  2y  y  -  1  ^ 

'y-1 

1  •p 

Cp  cp^ 

Observe  also  that  y  =  1  does  not  play  the  same  role  as 
y  *  00  or  y  =  0  ;  accordingly  ,  y  =  1  must  be  considered  separately. 


56,  Parametric  poles  .  Set  y  =  ^  ,  z  =  1  +  u  z  and  deter- 

z 

2 

mine  s  by  2es  =  1  ;  then  p  =  0  and 
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z  u  =  ^  ~  ^  +a+('e-f)s  +  0(z)  , 

s 

The  condition  for  stability  is 

s  -  1  +  as  +  (e  -  f)  s^  =  0  ; 

•  p 

therefore  ,  s  +  as  =  0,  (e-f)  s  *=l  and  finally 

(56.1)  e  =  2  e  a  , 

(56.2)  e  +  f  =  0  . 

Note  that  for  P2  ,  e  =  f  =  0  so  that  (56,2)  still  holds  , 

Parametric  zeros  ,  From  the  equivalence  table  ,  one  deduces 
at  once  the  conditions  for  stability 

(56.3)  g  =  2  b  g  , 

(56.4)  g  +  h  =  0  . 

Note  that  for  Z2  ,  g  =  h  =  0  so  that  (56.4)  still  holds  . 

57.  Parametric  unities  .  We  have  to  consider  separately  the 
parametric  unities  of  the  first  order  and  the  parametric  unities 
of  the  second  order  . 

For  the  parametric  unities  of  the  first  order  ,  set  y  =  1+sz, 
z  =  1  +  u  z  ;  one  has 

# 

*  s 

y  =  s  ( 1  +  zu  +  _  z )  , 

s 

y  -  s(zu+u+2  ^  ^  z  +  2  ^  zu  +  zu^)  , 

s  s  s 

^  ^ ^  ^  ^  +  O(z^)  , 

2y  y  -  1  sz  2  2  2 

ay  -  b  ,  a  -  b 

_ =  a  +  _ 


(57.1) 


y  -  1 


sz 


9 
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■(y-1)  Q(x»y)  =  ^  +  k  +  1  +  (e+f+g+h+-2  )sz 


sz 


(57.2) 


+  (2e  +  f  ~  g)s^z^  +  O(z^)  . 

H 

Because  of  (56.2)  and  (56.4)  *  one  finds 

« 

zu  * 


a  -  b  ^  k 


+  1  1  +  ^  \  u 


(57.3' 


+  +  A^z  +  uz(a  +  s)+0(z'')  , 


JJthe  coefficient  of  z  is  zero^  where 
A  = 


=  +  a  +  (a  -  b)  + 


A.  =  -  +  a  4- 

s  s 


£  - 


,,2  .  *2 

®  4.  o  a.  ® 

+  S  +  _ 


We  determine  s  by  setting 


1  +  ^  -  fa  +  ^  =  0 


SL  •  b 

p  =  1  +  _  must  be  an  integer  .  Therefore  ,  in  order  that 

s 

(55.1)  be  stable  ,  one  must  have  one  or  the  other  of  the  following 
distinct  possibilities 


i.  a  -  b  ^  0 


k  =  0 


p  =  0 
p  =  1 


s  =  b  -  a 


s*^  +  k 


il.  a-b=0  ,  k^O  , 

iii.  a“b=0  ,  k=0. 

[["see  (16.15)  and  (I6.I6)  ;  observe  that  n  =  00  and  change  s 
into  ^  ^  . 


Case  iii  corresponds  to  parametric  unities  of  the  second 


order 


When  a  -  b  =  0  ,  one  may  assume  a  =*  b  *  0  ["  choose  qs  such 
that  q)  =  a  qi  . 

To  determine  the  condition  for  stability  corresponding  to 
double  parametric  unities  ,  observe  that  because  a  =  b  -  0  , 


e  and  g  are  constants  see 

(56.1)  and 

(56.3)]  . 

The  equation 

(55.1)  is  then 

(57. y  ^  ^  ^ 

\  y^+y(y-i) 

re{y-l)+g 

( J_-i)  +  -L 

V  2y  y  -  1  . 

I 

W  1 

>-  ,2  ^  y  - 

Now  ,  set  y=l+sz,z=l+zu  and  determine  s  by  l+2s=0 
then 


y  z  u  =  u  + 


s 

s 


+  O(z^)  +  0(zu)  , 


The  condition  for  stability  is 


s 

s 


_  =0  or  s  =  He^ 


and  finally  , 
(57.5) 


p  ^  „  -Kx 


H  e 


56.  Because  y  =  00  and  y  *  0  play  the  same  role  ,  we  have 
to  consider  the  following  nine  distinct  possibilities  ; 

P2  ,  ZZ  and  U2  or  U1  ,  p  *  0  or  U1  ,  p  =  1  ; 

PI  ,  Z2  and  U2  or  U1  ,  p  =  0  or  U1  ,  p  =  1  ; 

PI  ,  Z1  and  U2  or  U1  ,  p  =  0  or  U1  ,  p  =  1 

We  consider  these  cases  separately  and  begin  with  the 


last  three 
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56.  PltZltUl.p  =  0.  Equation  (55.2)  is 


(56.1) 


y  =  (  — i-  +  — i — y 

\  2y  y  -  1  / 


y2  +  ay  -  b 


y  -  1 


y  + 


+  y(y-l)  r  e(y-l)  +  ^  (1-y)  +  _ 

y^ 


Vi- 


The  condition  for  stability  is  =  0  or  because  s  =  b  -  a  , 


(56.2: 


b  +  a-  b+  ^  (b‘^-a'^)=0. 


To  simplify  the  notations  ,  set 
a 

then  ,  from  (56.2)  ,  (56.1)  and  (56.3)  » 


2  2 

2c  ,  b  -  2d  ,  2e  =  q’  ,  2g  =  r 


(56.3) 


^  q  =  2cq  ,  r  +  2dr  =  0. 
Equation  (56.1)  may  be  rewritten  as 


. (56.4) 


(  1  ^  1  \ 

\  ^  2(cy  +  d)  ;  ^  (y-l)2  / 

V  2y  y  -  1  1 

^  y-1  2  ' 

2 

q  y  - 


-■  <tve  -  d  -  c  -  d)  y  . 


To  integrate  equation  (56.4)  »  observe  that 
V  =  ^  +  2(c  +  d) 

y-1  y  -  1 

is  regular  for  U1  ,  that  v  +  q  y  is  regular  for  PI  and  that 
v+qy+t+2(q+r )  is  regular  for  Z1  .  Accordingly  ,  set 


(56.5) 


2(c  +  d] 


+  q  y  =  -  2  w 


y-1  y-1 

on  substituting  in  (56.4)  »  one  finds 
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i 


(56.6) 


y  -  1 
4y 


w  +  2  w  d 


L  write  y  =  -  2w(y  -  1)  -  2(c+d)y  -  qy(y  -  1)^  . 

The  differential  system  (56.5-6)  is  thus  equivalent  to  the  equa¬ 
tion  (56.4)  5  if  y  is  stable  ,  w  is  also  stable  and  conversely. 

On  eliminating  y  between  (56; 5)  and  (56,6)  ,  we  obtain 
an  equation  for  w  .  For  convenience  ,  set 
A*w+2dw  ,  B=/4.w^-r^ 


so  that 


B  +  4  A 


y  -  1 


4A 

B  +  4  A 


Now  ,  rewrite  (56.5)  as 

+  2w  +  2(c+d)  +  q(y-l)  -  ^  «  0 

y  y 

and  observe  that 

B-6Aw  =  -4Bd  ; 

then 

(56.7)  B+4A+4Aq+  2(d  -  c)  (B  +  AA)  =  0 


or 

*  *  B 

w  =  -  2  w  w  +  a,  w  +  w  +  w  +  a. 

12  3  4 

i.e,  an  equation  of  type  A(x,y)  “  0  and  class  PI  ,  p  «  1  ;  thus  , 
w  and  y  are  stable  . 

To  integrate  (56.7)  »  observe  that  because  q  =  2  c  q  ,  one 
also  has 


Aq  =  qw+2dqw  = 


^  (qw)  +  2(d  -  c)qw 
dx 
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hence  ,  on  setting 


4Z  =  B  +  <4.A  +  4qw 

=  4^/'^  +  +  w(q  +  2d) 


f 


(5^*7)  becomes 

Z  =  2(c  -  d)  Z  . 

Therefore  ,  Z  is  determined  by  Z  =  K  q  r  and  w  by  a  Riccati 
equation  . 


59.  Pl;Zl;Ul.p  =  1.  Equation  (55.2)  is 

(59.1)  y  =f  2. _  +  — +y(y-l)  [  e(y-l)  +  ._i_(l>y) 

V  2y  y-1  /  „2 

^  ^  k  ^  X 1  "7~ 

where  e  and  g  are  constants  . 

p 

One  has  k7^0,p=l,s  +k  =  0  and 

•  2 

zu  =  u  +  A  +AtZ  +  suz  +  0(z  ) 
o  1 

the  coefficient  of  zu^  is  zero^  ;  the  condition  for 
stability  is 

^  r  ^  -  J_  +  si  =0 

dx  L  2s  s  J 

or 

2 _ 1,  =  0  ,  _l_f  ^  +  +  s 

ds  s  dx  2s 

Therefore  , 

(59.2)  s  =  H  e*^ 
and  k  +  2  2  +  2s^  =  K^^s  or 
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(59.3)  =  Kj  s  -  =  Kj  s  +  k  . 

Now  ,  we  consider  two  cases  according  as  K  is  or  is  not  zero  , 
i.  K  =  0  ;  then  s  ,  k  ,  2  {  and  e  ,  g)  are  constants  . 

Equat ion  ( 55 . 2 )  is 


(59.4)  y  J. _ +  ^  ]  y^+y(y-l)r  H(y-i)+H 

\  2y  y-1  j  ^ 


1  +  ^2  +  1 
/  r-i 


(y-1) 


24 


This  equation  may  be  integrated  by  using  the  method  of  varia¬ 
tion  of  parameters  ;  one  finds 


H, 


=  2y(y  -  1)^  r  Hy  -  ^^1  -  _  _ 

y  y-1  (y.i)2 

As  a  particular  case  ,  one  has 


+  H, 


1- 


y  = 


1  +  1 


2y 


y  -  1 


*2 

y 


whose  integi*al  is  y  =  Th  (H  x  +  H^)  . 

ii.  K  ^  0  .  By  a  transformation  t  =  <P  (x)  ,  one  may  assume 
e^*  =  t  and  k  =  2Ht^  »  B  ~  Ht^  +  H^t  ;  hence  ,  on  writing 
again  x  instead  of  t  ,  one  has 

V  J 


(59.5) 


+  H^xy  + 


Hx 


2  y(y  +  1) 


y  -  1 


This  is  the  irreducible  equation  IV  of  Table  I. 

A  transformation  t  =  op  (x)  brings  this  equation  to 


(59.6) 

^  H^y  Hy(y  +  1) 


X 


y  -  1 
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60.  P1:Z1;U2.  Equation  (55.2)  Is 

(60.1)  y  ^ 

where  e  and  g  are  constants  . 

The  condition  for  stability  Is  P  —  H  ^see  (57.5)1  • 
Equation  (60.1)  is  thus  a  particular  case  of  (59.4)  or  (59.5). 


f  1  '  U^y(y-l) 

\  e(y-l)+  ^  (1-y)  + 

\  2y  y  -  ij 

y2  y  - 

6l.  On  using  the  preceding  results  ,  the  other  cases  are 
readily  disposed  of  . 

P2;Z2 ;U2.  One  has  ,  for  P2  ,  e  -  f  =  0  ;  for  Z2  ,  g  =  h  =  0 
for  U2  ,  a=b=k=0  .  Equation  (55.2)  is 


(61.1) 

y  =  [  ^  ^  \  y^  +  £  y 

V  2y  y  -  1  / 

with  the  condition  for  stability  2  “  H  e^^  ;  (6l,l)  is  a 
particular  case  of  (55.4-5)  . 

P2;Z2;Ul.p  =  0  .  Equation  (55.2)  is 


and  is  a  particular  case  of  (5^.1)  . 

P2;Z2;Ul.p  =  1  .  Equation  (55.2)  is 

(61.3 

and  is  a  particular  case  of  (59.1)  . 

P1:Z2;U2.  Equation  (55.2)  is 

(61,4)  y  =(  ^  \  y^  -t  y(y-l)f  e(y-l)  +  ^  ~| 

V2yy-ly/  L  y-l-* 
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where  e  is  a  constant  .  This  equation  is  a  particular  case  of 

(60.1)  . 

Pl:Z2;Ul.p  =  0.  Equation  (55.2)  is 


(61. 


5)  y  =/  ^  \  y^+  ~  ^  y  +y(y-l)r  e(y-i)+  2 

v2y  y-i/  y-i  y-i 


where  e  =  2  e  a  .  This  equation  is  a  particular  case  of  (5^.1). 
P1;Z2!U1.p  =  1.  Equation  (55.2)  is 
1 


(61.6)  y  =, 


_  +  __ 

2y  y 


J; _ \  y^+y(y-l)  I  e(y-l)H-  ^  +  '2 

-ij  (y  .  1,2  y  -  1 


which  is  a  particular  case  of  (59.1) 
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IX.  Equations  of  the  type  A(x,y) 


2 

T 


y  -  1 


62,  The  stable  equations  of  this  type  are 


(62.1) 

where 

(62.2) 

(62.3) 


y  = 


y  -  1 


+  P(x,y)  y  +3y(y-l)  Q(x,y)  , 


P(x,y)  =  ay  +  _  -  ^ 

y  y  -  1 

Q(x,y)  =  ey  +  f _  -  _ 


+  d 


e 


+  h  +  4- 

,2  y  y  -  1 


y"  (y  -  i)‘ 

The  values  of  y  for  which  the  general  existence  theorem 
does  not  apply  are  0  ,  1  ,  oo  , 

The  transformation  y  =  u”^  [l^resp.  y  =  1  -  wj  changes  y  =  0  , 
oo  ,  1  Into  u  =  oo  ,  0  ,  1  [[resp.  w  =  1  ,  oo  ,  oj  and  brings 
equation  (62.1)  to  an  equation  of  the  same  form  according  to 
the  equivalence  table 

(y):abc  d  e  f  g  h  k  ^ 

(u)  :  b  a  -c  d+c  f  e  g  k  h-iP-g 
(w)  :  -a  c  b  a+d  e  g  f  -e-h  £  k 
Therefore  ,  the  values  y  =  0  ,  1  ,  oo  play  the  same  role  , 

Note  also  that  the  transformation  defined  by  (l-y)u  =  1 
brings  equation  (62,1)  to  an 'oquatlofi  of  the- s$me  fora 
according  to  the  table  of  equivalence 
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(y)  ;abc  d  efgh  k  £ 

(u)  ;  c  -a  -b  a+b4d.  g  e  8  -h-e  -f-h  , 

The  solution  y(x)  of  (62.1)  has  simple  parametric  poles  if  a 
and  e  are  not  both  zero  ;  y(x)  has  double  parametric  poles  if 
a  =  e  =  0  . 

Suppose  y(x)  has  a  simple  parametric  pole  .  Set  y  =  ^  ,  where 

z 

2  2  1 
s  is  given  by  3es  -  as  -  =  0  ;  then  n  =  -  3  and  p  ==  _____  +as. 

3  3 


P 

In  order  that  y(x)  be  stable  ,  one  must  have  e  =  3  a  and  p  =  0, 
3as  +  l==0orp=l  ,  3as  =  2  . 

If  y(x)  has  a  double  parametric  pole  and  is  stable  ,  then 
p  =  -  ^  so  that  h  =  0  .  Note  that  when  a»=e  =  h*0. 


y(x)  has  in  fact  parametric  poles  of  the  third  order!  . 


for  PI  : 


for  Z1 


To  sum  up  ,  according  to  the  equivalence  tables  ,  one  has 
_2 


f  *  b‘ 


for  U1  :  g  =  c 
Therefore  , 


for  P2  ; 
for  Z2  : 
for  U2  ; 


a  =  e  =  h  =  0  ; 

b  =  f  =  k  =  0  ; 

c  =  g  =  £  =*  0  . 


(62. k)  Q(x,y)  *  a^y  +  ^ 


£ 


+  h  +  + 

,2  y  y  -  1 


y"'  (y  -  1) 

where  h  =  0ifa  =  0  ;k  =  0ifb  =  0;£  =  0if  c*=0  . 

The  transformation  t  =  cp  (x)  brings  equation  (62.1)  to 

(62.5)  y”=  y’^+  (  P(x,y)-  t  \  Q(x»y) 

2  \  y  y  -  iJ  V  /  ^ 

on  choosing  cp  such  that  op  =  d  ,  one  may  assiune  d  =  0  . 
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63.  Because  y  =  oo  ,  0  ,  1  play  the  same  role  ,  we  have 
to  consider  only  four  distinct  possibilities  ,  namely 
(P2;Z2;U2)  ,  (P2;Z2;U1)  ,  (P2;Z1;U1)  ,  (P1;Z1;U1). 

Equation  (62.2)  corresponding  to  (P2;Z2;U2)  is 


for  one  may  assiAme  d  =  0  by  a  proper  choice  of  <p  ,  Then  ,  y 
is  given  by 

P  =  K  y^(y  -  1)^  ; 

y(x)  is  thus  an  elliptic  function  and  is  stable  . 

^  On  setting  2y  =  1  +  u  ,  one  has  u^  =  K(u^  -  1)^3  . 

For  the  other  cases  ,  at  least  one  of  a  ,  b  ,  c  is  not  zero. 

We  may  assume  c  ^  0  .  Indeed  ,  the  transformations  y  =  u  , 

y=u“^  ,  y=  1-u  and  their  products  y  =  ^  ~  ^  ,  y  =  ^ 

u  1-u 

y  =  ^  bring  equation  (62,1)  to  equations  of  the  same  form, 

u  -  1 

and  according  to  the  equivalence  tables  ,  commute  a  ,  b  ,  c. 

We  consider  first  the  equation  (P1;Z1;U1)  . 


64,  For  a  simple  pole  of  y(x)  ,  set  y  =  ,  z  ®  1  +  uz 

z 

and  note  that 

P(x,y)  =  a  ^  +  d  +  (b-c)  ^  +  0{z^)  , 

z  s 

Q(x,y)  =  a^  ^  +  h  +  (k  +  0  )  ^  +  O(z^)  ; 


P  p 

then  ,  substitute  in  (62,1)  ,  determine  s  by  3a  s  -as  - 


2 


0 
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and  set  p  *  +  as  so  that 

3 

•  12  2 
zu’=pu  +  A  +  B2  +  Cuz  -  ____  zu  +  0(z  )  , 


where 


B 


A  = 


d  X 


^  +d-as+3(a^-h)s  , 


s 

”  / 


3s 

*2 

1 


2  ^  4 


^  +  f3(h-k-h, 


3s‘ 


C  = 


+  d  , 


3s 


i.  Suppose  p  =  0  ,  3a-s  +1  =  0,  The  condition  for  stability 
is  A  =  0  or 

(64.1)  h+a(a+d)-a=0. 

From  the  equivalence  tables  ,  one  has  the  two  other  con¬ 
ditions  for  stability  corresponding  to  Z1  ,  p  =  0  and  U1  ,  p  =  0 
respectively 

(  k  +  b{b  +c+d)-b=0  , 

(64.2)  \ 

C  +  c(a  +  b  +  c  +  d)-c  =  0. 

ii.  On  using  the  relations  (64.1-2)  ,  we  shall  reduce  our 
problem  to  another  of  the  same  type  involving  an  equation 
with  only  two  singular  values  0  and  co  .  To  do  this  ,  observe 

that  the  parametric  zeros  and  the  parametric  poles  of  y(x) 

* 

are  poles  of  the  first  order  for  ^  ;  hence  ,  the  singular 

.  y 

y 

values  of  are  oo  and  the  value  corresponding  to  y  =  1  . 

y 
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Now  ,  ^  -  3ay  is  regular  for  the  parametric  poles  PI  , 

•  y 

V  Id 

p  =  0  ;  +  3  is  regular  for  the  parametric  zeros  Z1  , 

y  y 

« 

p  =  0  and  ^  -  3c  is  regular  for  the  parametric  unities  U1  , 

y 

p  =  0  ,  Therefore  ,  ^  -  3a-y  +  -  3c  is  regular  for  PI  , 

y  y 

p=0;Zl,p  =  0;Ul,p  =  0. 

Accordingly  ,  set 

(64,3)  JL  =  3a{y  -  1)  +  3b(  1  -  +  3c  +  3v  ; 

y  ^  y  ' 


equation  (62,1)  gives 

(64,4)  y  -  1 


Q(v) 

'  *  . . 

V  -  P(v) 


where 

(  P(v)  =  v^  +  (2a  +  2b  +  2c  +  d)  V  , 

(64.5)  \  P 

J  Q(v)  =  2v‘"  +  3cv  . 

Equation  (62,1)  and  the  differential  system  (64,3-4)  are 
thus  equivalent  ;  if  y  is  stable  ,  v  is  also  stable  and  con¬ 
versely  . 

To  determine  v  ,  one  has  to  eliminate  y  between  (64*3) 
and  (64,4)  ,  To  do  this  ,  rewrite  (64,3)  as 

^  =*  3a(y-l)  +  3 (''•''  ’^)  +  3(a  +  ^3  +  ^  +  y)  ^ 

y  -  1  y  -  1 

take  the  logarithmic  derivative  of  (64*4)  and  observe  that 

3(v  +  c)  ^  1  ^  1  , 

q( v)  V  2v  +  3c 
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then 


V 

-  p  +  3^2  .  ;  +  3  c  ^  P  .  _ 

(a+b+c-d) 

V  -  P  2v  +  3c  2v  +  3c 

i  .  6  .  f 

V  is  a  solution  of  an  equation 

(64.6) 

V  =  +  B(x;v)  V  +  C(x;v) 

2v  +  3c 

of  type 

A(x;v)  =  ^ _  . 

2v 

To 

reduce  (64.6)  to  a  canonical  form  ,  set 

^  3c 

u  =  V  +  , 

2 

(64.7) 

a=2(a  +  b)-c  +  d  , 

(64.6) 

p  =  3c{a  +  b)  +  (c  +  2d)  - 

3c 

4 

“2 

so  that 

*  *  2 

v-P  =  u-  u  -  au  +  p  , 
Q(v)  =  2u|^u  -  ^  . 


A  simple  calculation  shows  that 

*2  -  3  -  . 

u  =  ^  +  (a+b+c+2d)  u  +  3(b-a-c) 

2u  2 

(64.9) 

2 

+u  r  o+9ac-  p  -  ^  ”  (c+d)]|  +  p  (a+b+c+2d)-  p  -  P  . 

2  2u 

Now  ,  the  transformation  t  =  q>  (x)  brings  equation  (62,1) 

to  (62.5)  .  If  u  is  replaced  by  w  (x)  ,  (64.3)  and  (64,4) 

remain  invariant  if  the  coefficients  of  (62.1)  are  replaced  by 

the  corresponding  coefficients  of  (62.5)  .  I^*  <P  (x)  is  chosen 
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so  that  27=  (a+b+c+2d)  7  ,  and  if  w  is  to  be  stable  (  and 

thus  also  y  )  ,  then  equation  (64.9)  must  reduce  to  the  canoni¬ 
cal  form  see  ^  52  ,  eq.  E.7  and  E.ieJ  • 

(64.10)  w  =  ^  +4Dw‘^+2Ew-  ^  , 

2w  2  2w 

where  D  ,  E  ,  F  are  given  by  (52.2)  and  (52.3)  . 

Therefore  ,  to  determine  a,b,c,d,  one  has 


(64.11) 

a  =  2(a  +  b)  -  c  +  d  , 

• 

(64.12) 

p  =  3c(a+b)  +  (c+2d)  - 

3c 

> 

4 

2 

(64.13) 

a  +  b  +  c  +  2d  =  0  , 

(64.14) 

3(b-.c-a)  =  4D  , 

(64.15) 

a+9ac-p  -  (c  +  d) 

2 

• 

=  2  E  , 

(64.16) 

P  =  0  , 

(64.17) 

P  =  F  . 

Set  3c 

=  2  W  ,  Then  ,  from  (64.12-^13) 

and  (64.17) 

(64.16) 

3(a  +  b)  =  2  ; 

W  W 

from  (64.14) 

and  (64.13)  respectively  , 

(64.19) 

3(b  -  a)  =  4D  +  2  W  , 

3d  =  -/^  +  ^  +w)  . 

V  w  w  / 

Therefore  , 

« 

a 

=  3(a  +  b  +  d)  =  ^  ^  ..  W 

W  W 

so  that  ,  from  (64.15)  » 
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*2 

(64.20)  W  =  T ^  W^+4DW^+2EW-  ^ 


2W 


2 

2~ 


65.  From  the  preceding  ,  it  follows  that  the  stable  equa 
tions  of  the  class  (P1,21,U1)  are 


y  = 


.(  -L  +  ay  t  _i.  -  _1_  t  d]  ; 

\  7  y-i;  y  y-l  J 


r.2 


+  3y(y-l)  a  y  + 


(y  -  1) 


^  +  h  + 

2  y  y 


JL  +  _L_1 

y  y  -  1  J 


where 


}a  =  ^  ^ 

-  2D  -  W 

9 

3b 

3; 

W 

W  W 

“IT 

3c  =  2  W  , 

3d  =  -  f 

w 

+  ^ 

F 

_  +  w 

V 

w 

"w" 

h  =  a  +  ^  (b 

+  c  -  a) 

2 

k  =  b  -  ^  (b 

+  c  -  a) 

9 

2 

£  =  c  -  (a 

+  b  +  c) 

f 

+  2  D  +  W  , 


)  ■ 


W  is  a  solution  of  (64.20)  and 


y  -  1  = 


2  w(  w  -  W) 

t  '4' 


w-w*'-  OW  +  p  w-W  -(w-W)(w+W+  a  ) 


because 


(3  =  F  =  (  a  +  W)W  -  W  . 


Therefore  , 
(65.1) 


•  • 


_  w-W  +(w  -  W)  (w  -  W  -  a  ) 

y  .  "j  . . 

w-W  -  (w-W)  (w+W+  a  ) 
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and  is  stable  . 


66,  Remark  .  For  PI  ,  p  =  1  ,  3as  =  2  ,  one  has 

zu=u+A+Bz+  Czu  -  ^  zu^  +  O(z^)  , 

3 

where 

Vi  * 

D  >=  +  a  +  d  so  that  a  =  a  D  ,  (  see  64.1)  ; 

a 

A  =  D  -  ,  3C  =  3d-6d-2D  , 

a 

B  =  -  D  -  _J_  -  (2a  -  d)D  +  a(b-c-a)  +  3(h  -  k  -  8  )  , 

3  2 

The  condition  for  stability  is 

,  »2 

A  +  B  A  C  -  _ _  =  0 

3 

or 

_ -  _J_  a{b-c-a)  +  k  +  8  -  hd  +  D  «  0 

dx  a  2  a  a 


2(ah  -  2ah)  +  (a^  +  2h)^  -  a^(b  -  c)  +  2a^(k  +  t  )=  0  , 

There  are  two  similar  conditions  for  stability  for  Z1  ,  U1  , 

namely 

2(bk  -  2kb  )  +  (b^  +  2k)^  -  b^(a  +  c)  +  2b^(h  -  0  -  c^)  =  0  , 

2(c  0*  -  2  0  c)  +  {c^  +  2  £  )^-c^(b-a)-2c^(h+k+a^+b^)  =  0  . 

These  relations  are  not  easy  to  handle  and  thus  the  reason 
why  another  method  was  used  . 
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67.  P2;Z2;U1.  One  has  a  =  h  =  0  for  P2  ,  b  =  k  =  0  for  Z2 
the  stable  equation  for  this  class  is 


where 

t=c-c(c+d)  . 

We  use  the  above  method  and  assixme  c  =  2d  =  0  ;  then  w  is 
determined  by  (64.10)  and  c  and  d  by  (64.11-17)  .  One  easily 
finds 

3c  =  2W  ,  W  +  F  =  0  ,  2D  +  W=  0.,  3d  =  -.W, 

a=-W  ,  P«F  ,  E  =  D^  ; 


therefore  , 

W  =  -  F  X  -  H  . 
Equation  (67.1)  is  then 


w 

y  + 

• 

^  -  2y  1  F  + 

y  + 

•  — 

1 

3  1 

y  y  -  1/ 

3 

y  - 

T  ^ 

3 

y  - 

T 

its  general  solution  is 


y 


1  + 


2w(w  -  W) 


w  —  w  +  Ww  +  F 


*  2 

w  +  w  -  Ww+F 

*  2 

w-w  +  Ww  +  F 


_  w  -  W  +  w(w  -  W) 
w  -  W  -  w(w  -  W) 


where  w  is  determined  by 

-2^  ~ 


w 


-  2  W  w^  +  2  w 


F*- 

2w' 
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6B,  P2 ;Z1 ;U1 .  One  has  a  =  h  —  0  for  F2  .  The  stable  equation 
for  this  class  is  equation  (62.1)  ,  where  a  *  h  =  0  .We  use 
again  the  above  method  and  ass\ime 

b  +  c  +  2d  =  0  f  3c  =  2W, 

Then  ,  one  has 

(68.1)  W=W^  +  2DV/-F 

so  that  W  is  determined  by  a  Riccati  equation  . 

Moreover  , 

3b  =  4D  +  2W  ,  3d  =  -2W-2D  , 

a  =  2  D  ,  P  =  F 

and 

2D  +  2D^-F  =  2E. 


The  stable  equation  is  thus 


y  = 


2 

T 


1 


+ 


*2 

y  + 


/  b 


+  3y(y-l) 


[ 


y  -  1 


I 


(y  -  l)‘ 


9 


where  k  =  b-fbd  ,  L=c  +  cd, 

The  general  solution  is  given  by  (65.I)  ,  where  W  is  a  solution 
of  the  Riccati  equation  (68,1)  and  where  w  is  the  general  solu¬ 
tion  of  (6/+., 10)  . 

Note  that  if  W  is  replaced  by  -  ^  in  (68.1)  ,  one 

u 


obtains 
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u-2xu-Fu  =  0  ; 

if  F  =  -  2n  ,  n  an  integer  ,  this  equation  is  satisfied  by 
H^(x)  ,  the  Hermit e  polynomial  of  order  n  . 
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X,  Equations  of  the  type  A(x,y)  *  ^  ^  ^ 

4  V  y  y  -  1 

69.  The  stable  equations  of  this  type  are 

(69.1)  y  =  ^  ^  +  __i_  1  y^  +  P(x,y)  y  +  y(y-l)  Q(x,y)  , 

4  V  y  y  -  1; 

where  P(x,y)  ,  Q(x,y)  are  given  by  (62.2-3). 

The  values  of  y  for  which  the  general  existence  theorem 
does  not  apply  are  0  ,  1  ,  00 . 

The  solution  y(x)  of  (69. 1)  may  have  simple  or  double 
parametric  poles  . 

If  y(x)  has  a  simple  parametric  pole  ,  set  y  =  ^  , 

z 

•  2 

z  =  1  +  uz  ,  where  s  is  given  by  2es  =  1  ;  then  n  =  -  2  (case  vi), 
a  =•  0  ,  e  ^  0  and  p  =  0  ,  Then  ,  substitution  in  (69»1)  yields 

zu  «  ^  +  d  +  (e  -  h)s  +  0(z)  . 

s  4s 

The  condition  for  stability  is  thus 

s  -  ^  +  d  s  +  (e  -  h)  s^  »  0 

4 
2 

or  ,  because  2es  =  1  , 

s  +  d  s  =  0  ,  4(e  -  h)  s^  =  3  ; 

finally 

(69.2)  2h+e=0  ,  e=2ed. 

If  y(x)  has  a  double  parametric  pole  ,  then  a  =  e  *  0 
and  ,  as  a  consequence  ,  h  =  0  . 


134 


Changing  the  notations  and  taking  into  account  the  relations 

(69.2)  ,  we  shall  rewrite  equation  (69.I)  as 

(69.3)  y  =  \  +  P*(x,y)  y  +  y(y'l)  Q*(x,y)  , 

4  V.  y  y  -  1/ 

where 


P*(x,y) 

Q*(x,y) 


4d^(2y-l) 


> 


Note  that 

(69.4)  d  =  a  d  X 

according  to  (69.2)  ;  for  PI  ,  one  has  d  ^  0  and  for  P2  ,  d  «  0. 

The  exceptional  values  y  =  0  and  y  »  1  play  the  same 
role  ;  Indeed  ,  on  setting  y  =  1  -  w  in  (69.3)  »  one  obtains 
an  equation  of  the  same  form  according  to  the  equivalence  table 
(y)  :  abcdfgkp 

(w)  :  a  c  b  d  g  f  £k. 

It  is  readily  seen  that  y(x)  has  simple  parametric  zeros 
if  b  and  f  are  not  both  zero  and  double  parametric  zeros  if 
b  *  f  =  0  and  k  ^  0  , 

In  agreement  with  the  equivalence  table  ,  one  sees  at  once 
that  y(x)  has  simple  parametric  unities  if  c  and  g  are  not 
both  zero  and  double  parametric  unities  ,  if  c  =  g  =  0  arid 
£  0  ,  Therefore  ,  we  have  to  consider  only  six  distinct 


possibilities  namely 
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(P2;Z2;U2)  ,  (P2;Z2;U1)  ,  (P2;Z1;U1) 
(P1;Z2;U2)  ,  (Pl;Z2;Ul)  ,  (P1;Z1;U1) 

We  begin  by  considering  the  case  (P1;Z1;U1)  . 

Furthermore  ,  note  that  a  transformation  t  = 
equation  (69.3)  to 
7"  - 

4  V  y  y  -  1 


(p  (x)  brings 


y(y-l) 


* 

Q 


(x,y) 


70.  Suppose  that  y(x)  has  simple  parametric  zeros  .  Set 
y=sz  ,  z=l+uzso  that 

P*(x,y)  =  ^  +a+c+csz+  0(z)  , 

sz 


Q*(x,y)  =  ^  ^  -  (ifd^  +  e  +  £  )  +  sz(Sd^-2g-  £  ) 

2  2  sz 

s  z 


-  s^  z^(,3g  +  £  )  +  0(z3)  . 

On  substituting  in  (69.3)  »  one  obtains 

(70.1)  z  u  =  ^  u  +  A  +  0(z) 

z  \  k  s  j  \  Z  s  J 


where 


(70.2)  A  -  a  +  c  +  b  ^  .  -  JL  +  _£ 

g2 

Now  ,  determine  s  by 


then  p  =  +  _____  must  be  an  integer 

2  s 

to  Appendix  I  ,  vl  Fsubstitute  s  with 


s 


Therefore  ,  according 
,  one  finds 
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(70.3) 

and 

p  =  0 

For  Z1  ,  p 

1 «  a .  f 


f  = 

,  2b  +  s  =  0  ;  P  =  2  ,  2b=3s 

0  f  the  condition  for  stability  is  A  =  0  ,  s 


(70.4)  k==2b-^Sb(a  +  b4-c)  . 

In  agreement  with  the  equivalence  table  ,  one  obtains  for  U1  , 
p  =  0  ,  the  condition  for  stability 


(70.5) 


( 


g  =  c  , 

^  =  2  c  “  2c (a  +  b  +  c)  . 


71.  Before  considering  the  case  Z1  ,  p  =*  2  ,  we  find  the 
condition  of  stability  for  Z2  . 

For  a  double  parametric  zero  ,  one  has 

(71.1)  b  =  f  =  0  ,  kT^O. 

2  * 

On  setting  y=sz  ,z=l+uz,  one  finds 


a  -  ^  =a+c+csz^+  O(z^)  , 

y  -  1 

Q*(x,y)  =  ^  -  (  ^  +  g  +  4d^)  +  O(z^) 

2 


and 


z  :  =  _i_/l  - 
2z  V 

Note  that  p  *  0  1  . 


^  +  a  +  c  -  ^  ^  +  0(z) 

s  y  2  s 


Determine  s  by  s  =  k  ;  because  p  *  0  ,  the  condition  for 


2b  , 


stability  is 


; 
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(71.2)  k  =  2k{a  +  c)  . 

For  a  double  parametric  unity  ,  the  condition  for  stability  is  , 
according  to  the  equivalence  table  , 

(71.3)  c=g=0  ,  ,  £=2£(a+b). 


72.  PI ;  Z1  ;U1 .  We  have  d  0  for  PI  .  On  using  the  conditions 
(70.3-5)  t  we  reduce  the  problem  to  another  one  involving  an 
equation  with  only  one  singular  value  ,  namely  oo  . 

V  2b 

To  do  this  ,  note  that  +  _  is  regular  for  Z1  ,  p  =  0 


and  that  ^  ~  is  regular  for  U1  ,  p  =  0  [  set  1  -  y  *  sz 

1  -  y 

and  note  that  s  =  -  2c  ,  V/e  don’t  pay  particular  attention  to 

PI  because  one  ought  to  consider  ^  ^  with'^s  =  +  ^  ; 

y  s  "■  4d 

•  • 

^  +  4dy  is  regular  for  a  set  of  parametric  poles  and  ^  -  4dy 

y  y 

is  regular  for  the  other  set  of  parametric  poles ^  . 


Now  ,  set 


y  +  +  y  ~  =  4  V 


y  1  -  y 

on  substituting  in  (69.3)  >  one  finds 

(2y  -  1)  (v^  -  d^)  =  V  -  a  V  . 

Hence  ,  the  differential  system 

(72.1)  y  =  2b(y  -  1)  +  2  c  y  -  4  v  y(y  -  1)  , 

(72.2)  2y  -  1  *  V  -  a  V 

v^  -  d^ 
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is  equivalent  to  equation  (69.3)  ;  if  y(x)  is  stable  ,  v  is 
also  stable  and  conversely  . 

To  determine  v  ,  eliminate  y  between  (71.1-2)  .  To  do 
this  ,  set 

v-av’=A  ,  V  -  a  —  B 

and  note  that  d  =  a  d  [|see  (69.4)1 
B-2vA  =  2aB  ; 


then  , 


A  =  2(a+b+c)A  +  2(v  -  b  +  c)B 


or 


(72.3) 


V  =  (3a  +  2b  +  2c)  V  +  2v^  +  2(c  -  b)v^ 

+  a-2a(a  +  b  +  c)  2d^v  +  2d^(b  -  c) 


i.e,  ,  an  equation  of  type  I  see  equation  20,5  and  §  24]  . 

1 


2/1 

Now  ,  by  the  method  used  in  type  _  _ 

3  y 


+ 


y  -  1 


[]i.e.  ,  a  proper  choice  of  q>  "]  »  one  may  assume 
(72.4)  3a+2b+2c=0 

so  that  ,  on  setting  c  -  b  =  3  ^ay  rewrite  (71.3)  as 

V  =  2v^  +  6  D  v^  +  (a  +  a^  -  2d^)v  -  6  D  d^  , 

The  transformation  v  =  w  -  D  brings  this  equation  to  the 
canonical  fora 

w=2w^+(a  +  a^-'  2d^  -  6  D^)  w 


(72.5) 


+  D  +  4D^  -  (a  +  a^  -  2d^)D  -  6  D  d^ 


In  order  that  w  be  stable  ,  equation  (72.5)  must  ba  identical  to 
(72.6)  V=2v3+SV+T 
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so  that 

(72.7)  a  +  -  2d^  -  6D^  *=  S  , 

(72.6)  D  +  -  (a  +  a^  -  2d^)  D  -  6  D  d^  *  T  . 

On  using  d  =  a  d  ,  one  obtains  from  (72.7)  and  (72,6)  after 
simple  calculations 

(72.9)  d-2d^-Sd-6D^d=0  , 

(72.10)  D-2D3-DS-T-6Dd^=0. 

From  (72,9-10)  ,  it  is  clear  that  D  -  d  *  and  D  4-  d  = 
are  solutions  of  (72.6)  ;  therefore  ,  on  using  (69.4)  and  (72.4), 
one  finds 


2  d  *  Vg  -  Vj  ^  0  , 


(72.11) 


D  =  c  -  b 


1  ) 


(Vi  +  Vg)  , 


f 


•  • 


The  general  solution  of  equation  (69.3)  ,  where 


(72.12) 

is  thus  given  by 


f  -  b^  ,  g  =  e^  , 

k  =  2b  +  ab  ,  {^=2c  +  ac 


2  y  -  1  = 


(72.13) 


w  -  D  -  a(w  D) 
(w-D-d) (w-D+d) 


-  .  a  V  -  Vj  -  -  a(2v  -  .  v^)  _ 

2(V  -  Vj){V  -  Vg) 

where  V  is  the  general  solution  of  (72.6)  ;  therefore  y  is 


stable 
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73. 

(73.1) 


P2 ;Z1 ;U1 .  For  P2  ,  one  has  d  =  0 


2y  -  1 


V  -  a  V 


then  (72,2)  becomes 


On  eliminating  y  between  (72.1)  and  (73*1)  >  one  finds 
again  (72.3)  ,  where  d  =  0  , 

The  method  used  above  shows  that  w  defined  by  v  =  w  -  D  , 
3D  =  c  -  b  ,  33-  +  2b  4-  2c  =  0  satisfies 


w  =  2  w^  +  (a  +  a^  -  6D^)w  +  D  +  415^  -  (a  +  a^)D  . 

In  order  that  w  be  stable  ,  this  equation  must  be  identi¬ 
cal  to  (72,6)  so  that 

(73.2)  a  +  a^  -  6  D^  =  S  , 

(73.3)  fi  +  4  -  (a  +  a^)  D  *  T  . 

Hence  ,  D  is  a  solution  of  (72,6)  and  a  is  given  by  the 

Riccati  equation 

a  +  a^  =  S  +  6  V? 

X  • 

or  by  the  linear  equation  f  set  a  =  |  , 

u  / 

(73.4)  u  =  (S  +  6  V^)  u  . 

The  general  solution  of  the  corresponding  equation  (69.1) 
is  given  by  (72,13)  and  is  stable  . 


74.  P2  ;  Z2  ;  U2  ,  Because  d>=k  =  B=0  ,  the  stable 
equation  for  this  class  is 


(74.1) 


y 


3  /  1  4.  1  \  '2  ^ 

_  _  +  _  y  +  a  y  . 

4  V  y  y  -  I  J 


1 


Ul. 

By  a  transformation  f  ,  one  may  assume  a  =  0  .  The  solution 
of  (74.1)  is  given  by 

y^  =  K  y^(y  -  1)^  ; 

therefore  ,  y(x)  is  an  elliptic  function  and  is  stable  . 


75.  PI  ;  Z2  ;  US  .  The  stable  equation  for  this  class  is 

(75.1)  y  =  ^  ^  ^  \  y^+ay+y(y-l)  j  4d^(2y-l)+  ^  ^  j 

4\yy-lj  yy-l"^ 


k  \  y  y  -  IJ  y; 

By  a  transformation  9  ,  one  may  assume  a  =  0  ;  consequently  , 
d  ,  k  ,  C  are  constants  , 

To  find  the  solution  of  (75.1)  ,  set  y  *  and  obtain 
(75.2)  w  =1  _} _  +  ^  1  w^+(w-l)f  4d^f  ^  +kw  -  f 


'  ^  4- 

3 

i  +kw 

w  2w 

4(w~l) 

J  -  \  w 

/ 

The  solution  of 


w  =  r  +  3  1 

L  2w  4(w  -  1)  J 


2w  4(w  -  1)  J 
is 

(75.3)  =  K  w(w  -  1)3/2  ^ 

On  using  the  method  of  variation  of  parameters  ,  one 
determines  K  as  a  function  of  w  by 

JL  -  ^  r  (a  -  w)  +  k  -  _L_1  . 

dw  TTTTTP^  *-  •  W  -  1  J 


(2  -  w)  +  k  - 


Now  ,  set  w  =  1  +  andt>(u)  “  K(l  +  u^)  ;  then  , 

-  4  r  4d2.  ^  ^  4  k  +  J.] 

Hn  I-  ?.P  ?  -i 


(1  -  U^)^ 


12+2 


and 

%{m)  *  2+  r  k  u  +  ^  ^  ^  +  H  1  . 

u  n  ^  2 

1  +  u 

Thus  ,  u  is  an  elliptic  function  and  y(x)  is  stable  , 


76.  PI  ;  Z2  ;  U1  .  For  convenience  ,  we  slightly  change 
the  notations  and  replace  a  by  2a  in  P  (x,y)  ,  Then  ,  the 

conditions  for  stability  are 

# 

i.  for  PI  ,  d  =  2  a  d  , 

ii.  for  Z2,  b=f=0  ,  k  ^  0  , 

k  =  2k(2a  +  c)  , 

iii.  for  U1  ,  £  =  2c  -  2c(2a  +  c) 

so  that  it  remains  only  to  consider  the  condition  for  stability 
for  U1  ,  p  -  2  , 

Equation  (69.3)  inay  be  rewritten  as 


(76.1) 


+ 


y(y-l)  [  2+d^(2y-l) 


2a  - 


Now  ,  set 

(76.2)  y=2cy-4dy  (y-  1)  -  2v(y  -  1) 

and  substitute  in  (76,1)  ;  on  taking  into  account  the  values 

of  d  and  £  ,  one  obtains 

2  . 

(76.3)  =  2v  +  2v^-2+(a  +  c  +  d)v  +  k. 


y 
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On  eliminating  y  between  (76.2-3)  »  one  finds 


V  - 


*2 

■yr 


=  V  +  2(d  -  a)  - 


4E  ^  k 


■1- 


(76.4)  +  _f_  -  f^Zc-Zd  -  j 

» 

+  ^  E  k  -  (c  +  2d)  k  -  - 

~  2 


^  -4E(c+2d)-2d 


k^ 

” 


where  E=a+c+d. 

This  equation  (76.4)  is  of  type  III  and  clasj?^  E.19  S  to 
obtain  a  canonical  form  ,  we  use  the  usual  method  (  see  5  53  ) 
and  assume  b^^  +  5  “  0  (  notations  of  y  53  )  »  l.s*  » 

(76.5)  c  +  3  a  =  0  . 

Then  E==d-2a,k  =  -,2akso  that  (76.4)  may  be  rewritten 
as 


V  - 


*2 

2 


(76.6) 


+  2  ^3  +  10 


0*9  *  V  * 

^  V  V  -  ^  (a  +  d)  V  -  _  V 

3  3  3v 

( a+d ) v^- r  ^  +4(a  +  d)  -  ^  ( a-td ) ^ 


V-  3 


+  __f__(a  +  d)k  +  k-  ^ 


3v 


According  to  (53 *4)  »  set 


q=-(a+d)  ,  r-- 

then  f  equation  (76.6)  becomes 
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V  - 


*2 

2 


2  V  ;  +  2  +  r  ^  ^ 


qv 


(76.7) 


where 


+  b^  v^+  b3  - 


3r 


b,,  *  ^  q^  +  r  +  4  q  , 


2  q  r  -  3  r 


For  V  (  and  y  )  to  be  stable  ,  one  must  have  (  see  ^  53) 
q  =  2  q^  +  (  Kx  +  H)  q  +  , 

3r  =  Kx  +  H  +  2q  +  2q^  ; 

the  general  solution  of  (76.7)  is 


(76.8) 


•  *  2  2 

w-q  +  w  -  q 

V  —  .  > 


w  -  q 

where  w  is  a  solution  of 

w  =  2  q^  +  (Kx  +  H)  q  +  . 


77 •  PS  ;  Z2  ;  U1  .  For  P2  ,  one  has  d  =  0  .  The  stable 
equation  for  this  class  is 


(77. 


1)  y  =  _L.(  _i_  +  ^ \  y^+(  2a  -  _ \  y+y(y- 

4  V  y  y-1  y  V  7-1  J 


1) 


JL+-L.-. 

L  y  y-i 


(y 


where 

k  =  2k(2a  +  c)  , 

B  =  2  c  -  2c(2a  +  c)  . 

The  solution  of  (77.1)  is  obtained  on  setting 
y  =  2  c  y  -  2  v(y  -  1) 

and  is  given  by  (76.7)  and  (76.6)  ,  where  q  ** 


-  a  ,  r  =  - 
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^  1 

XI,  Equations  of  the  type  A(x,y)  =  +  . 

3y  2(y  -  1) 

78.  The  stable  equations  of  this  type  are 

(76.1)  y  =r  ^  ^  ,1  P  +  P(x,y)  y  +  y(y-l)  Q(x,y)  , 

^  3y  2{y  - 

where  P(x,y)  ,  Q(x,y)  are  given  by  (62.2-3)  • 

The  values  of  y  for  which  the  general  existence  theorem  does 
not  apply  are  0  ,  1  ,  co  . 

The  solution  y(x)  of  (78.1)  may  have  simple  or  double 
parametric  poles  , 

If  y(x)  has  a  simple  parametric  pole  ,  set  y  =  ^  , 

z 

•  25 

z  =  1  +  u  z  ,  where  s  is  given  by  es  -  as  -  =  0  and 

’T” 

2 

p  =  _____  +  as  ;  then  n  =  -  6  (  see  14.4  >  case  vii)  and 
3 

2 

3e  =  3a  ,  Accordingly  ,  y(x)  may  have  two  families  of  simple 
pareimetric  poles  corresponding  respectively  to  p  =  0  ,  3as  =*  -  2 
and  p  =  4  ,  3as  =  10  .  Then  ,  substitution  in  (78.1)  yields 

zu  =  ^  ^  ^  +d-as+(e-h)s+0(z)  . 

3  s  2s 

The  condition,  for  stability  for  p  =  0  is  thus 

(78.2)  2h  =  3a  -  a^  -  3ad  , 

2 

If  y(x)  has  a  double  parametric  pole  ,  then  a  *  e  =  0 
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and  as  a  consequence  h  =  0  Qsee  §  17j  . 

Consider  now  y  =  0  ;  the  parametric  zeros  of  y(x)  may  be 
simple  or  double  . 

If  y(x)  has  a  simple  parametric  zero  ,  set  y  =  sz  , 
z=  1+uz  in  (76.1)  and  determine  s  by 


2 


3 


+ 


0 


Then 

V  •  • 

1  .  2  s  s,,,  .bs,f-k,-/'  v 

+  _ I  u-  _ _ -  _  +c+d  +  ____  +  _  +  0  ( z  )  . 

3  s  /  3  s  Z  g2  s 

Because  p  =  ^  ^  must  be  an  integer  ,  one  sees  at  once 

3  s 

see  ^  l6  ;  note  that  n  =  -  3I  that  f  =  3b^  and  that  y(x) 
may  have  two  families  of  simple  parametric  zeros  corresponding 
respectively  to.p=0  ,  s+3b=0  and  p  =  1  ,  2s  =  3b  . 

The  condition  for  stability  for  p  =  0  is 


(76.3)  k  =  3b-_^b^-  3b(c  +  d)  , 

2 


If  y(x)  has  a  double  parametric  zero  ,  then  b  =  f  =  0  and 
as  a  consequence  k  ==  0  . 

Finally  ,  consider  y  =  1  ;  the  parametric  unities  of 
y(x)  may  be  simple  or  double  . 

If  y(x)  has  a  simple  parametric  unity  ,  set  y  =  1  +  sz  , 
z  =  1  +  z  u  so  that  (76.1)  yields 

f _ )  -  _1_  u  +  A  +  0  ( z  )  , 

1  ^ 

where 
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A  =  -  ^  ^  s  +  a  +  b  +  d+  ^  ~  ^  . 

°  s  3  s' 

Now  ,  determine  s  by  ^  ^  ^  =  0  ;  because  p  =  -  ^ 

2  s  2  s 

s 

■p 

must  be  an  integer  ,  one  has  c-0,g7^0,s'  =  2g  and  p  =  0,0. 
The  condition  of  stability  is  then  A  =0  and  according  to 
=  2g  , 

J__  s^  +  ^  -  g  =  0  ,  =  a  +  b  +  d 

3  s 

or 

(76.4)  3  E  +  g  *  0 

(76.5)  g  =  2g(a  +  b  +  d)  . 

If  y(x)  has  a  double  parametric  unity  ,  then  c  »  g  =  0  and 
consequently  E  *=  0  . 

Note  that  c  is  always  zero  .  Because  of  the  asymmetric 
role  of  y  =  0  ,  1  ,  oo ,  one  has  to  consider  eight  distinct 
possibilities  ,  namely 

(P2;Z2;U2  or  Ul)  ;  (P2;Z1;U2  or  Ul)  ; 

(P1;Z2;U2  or  Ul)  ;  (P1;Z1;U2  or  Ul) 

The  conditions  for  stability  are  given  above  except  for  PI  , 
p  =  4  and  Z1  ,  p  =  1  . 

We  begin  by  considering  the  case  (PI  ;  Z1  ;  Ul  )  . 

79.  PI  ;  Z1  ;  Ul  .  On  using  the  conditions  for  stability 
already  obtained  in  paragraph  76  and  ,  for  convenience  ,  changing 
■E  into  3  E  »  rewrite  equation  76.1  as 


ud. 


y 

,  “  ,1 

+  f  ay  +  ^  +d  )  y  + 

^  3y 

2(y  -  1) J 

V  y 

J 

+  3y(y-i) 

,  3?  1 

+ 

1 

. 

y^  (y  -  1)^ 

+  y(y-l)r 

h  +  + 

3l  T  , 

L 

y 

y  -  1  J 

,  k  and  g  are  given  by 

76.2-3-5  respectively  ; 

+  g 

=  0  . 

(79.2) 

Now  ,  we  reduce  the  problem  to  another  one  involving  an 
equation  with  only  one  singular  value  ,  namely  y  =  co  . 

To  do  this  ,  note  that  ^  y  is  regular  for  Pl,p  =  0 

y  2 

and  that  ^  is  regular  for  Zi  ,  p  =  0  ;  accordingly  ,  set 

y  y 


(79.3) 


y  =  y(y-l)  +  3b(y-l)  +  3vy 

2 


y  =  3b(y  -  1)  +  3w^y 


and  substitute  in  (79.1)  .  For  convenience  ,  set  first 
(79.4) 
and  find 

•  r- 

L'  * 


w. 


4^1 


w 


1  +  (b+d+ay)  w^^  +  ab(y  -  1) 


(79.5) 

+ 

then  ,  write 

(79.6) 


(y-i)[ 


.  3? 


+ 


-  +  £ 


(y  -  1) 


2  2 


2  y 


hi 


w.  «  _ (y  -  1)  +  V 

2 
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so  that 


.2  ^  ~ 

(79.7)  V  -  _  -(a+b+d)v  -  6  =  ^ 

2  y  -  1  L 


e 


Now  set 


a  +  b  +  d  =  -r 


(79.8)  ^  A-gv-v^  +  Srv-2^  , 

B  -  +  2  e  ; 

then  ,  according  to  (76,4-5-7)  »  one  obtains 

(79.9)  E  =  -  2  g  r 

(79.10)  y  -  1  =  ^  ^  . 

A 

To.  eliminate  y  between  (79.3)  and  (79.10)  ,  rewrite  (79.3) 


as 

^  =  3  +  3b  +  3v  +  _^(y  -  1)  +  ^  ^ 

y-1  .2  ,  2  y-1 

and  note  that 

(79.11)  B  -  V  A  -  (v  -  2r)B  ; 

then  V  satisfies 

V  =  -  vv  +  v-^  -  (a+2b+2r)v 


(79.12) 


V  -  ^  ^  (a  -  b)+rj 

-  ^  r  +  ^  ( a+2b ) r-2  E  +2r^ ]  ^  “  ^  ^  ^  ^  * 


To  integrate  this  equation  of  type  I  and  class  20.6  ,  set 

V  =  w  -  q  and  detei’mine  q  so  that  the  coefficient  of  w  be  three 

2 

times  the  coefficient  of  w  ;  one  finds 


(79.13) 
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(79  *14)  w  =  -  ww  +  w^+  ^  ^  +dj  (3w+w^)'  +  C^w  +  Cg  » 

where  ,  Cg  are  Independent  of  w  , 

By  the  method  already  used  for  the  equations  of  type 

^  M  -t  ^ 

3  \  y  y  -  1 


.  ~i 


j^i.e.  a  proper  choice  of  (p  ;  see  ^  64  J 


one  may  assume 

(79.15) 
so  that 


9a  *  -  10 (q  +  r) 


(  C, 


q  -  3q  -  (r  +  qr  - 


(79.16) 


Cg  *  q  +  2q^  4-  (r  +  rq  - 


If  equation  (79.14)  is  to  be 
may  have 


9b  *  4(2q  -  r)  , 

r^  .  2  i>  )  , 

r^  -  2  8  )q+6  f  q+2r  0 
stable  [*  see  ^  25.0"^ 


(79.17)  =  -  12  ,  Cg  -  12 

where  is  a  solution  of 

(79.18)  V=6V^+Kx+H  ; 


one 


equation  (79.14)  is  then 

(79.19)  w  =  -  ww  +  w^  -  12  V^w  +  12  . 

Now  set 

(79.20)  2  0  =  -  s^ 

and  observe  that  ,  because  of  (79.9)  » 

(79.21)  s  +  sr  =  0  ,  s  +  s(r.-r)*0 


one  has 
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s{r  +  qr  -  -  2  £  )  =  -  s  -  q  s  +  ; 

then  ,  from  (79.16-17)  ,  one  deduces 

s  =  -  s  q  -  s  q  -  12  s  +  +  3  s  q  , 

q  =  _  q  q  _  s  s  -  12  q  +  q-^  +  3  s  q  +  12 

so  that  q  +  s  and  q  -  s  are  solutions  of  (79,19)  .  According 

to  §  25. e  ,  one  has 


where  distinct  solutions  of  (79.  Id)  . 

The  solution  of  (79.1)  is  given  by 


2v-v^-t2rv+s^ 

where  r  =  -  ____  ,  v  =  w  -  q  and  where 
s 


is  a  solution  of  (79.19)  . 

do.  We  now  consider  the  other  cases  . 
i.  P2  ;  Z2  ;  U2  .  The  equation  is 

(dO.l)  y  =  r  ^  1  y^  +  d  y  . 

L  3y  2(y  -  1)  i 

A  transformation  y  -*  cp  (x)  may  be  chosen  so  that  d  *  0 
£  set  q)  =  d  <p  ]  ;  then  ,  the  integral  of  (dO.l)  is  given  by 
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(80.2)  =  K  y^(y  -  1)^ 

where  K  is  as  usual  ,  an  arbitrary  constant  ,  Therefore  ,  y(x)  is 
an  elliptic  function  of  x  and  is  stable  . 
ii.  P2  ;  Z2  ;  U1  .  The  equation  is 


(80.3) 


-[■ 


3y  2(y  -  l) 


y^  +  d  y  -  +  P  y 

y  -  1 


by  a  proper  choice  of  the  transformation  x  -»  cp  (x)  ,  one  may 
suppose  d  =  0  ;  hence  ,  see  78.4“5  J  g  =  0,g  =  -3B  » 

£  a  constant  . 

To  integrate  (80.3 )  »  where  d  =  0  ,  set  y  =  w”'*’  ;  then 


w=  [ 


w  2(w 


■^1 


,  3  g 


w’~  + 


w 


w  -  1 


-  £  w 


3 

or  on  setting  w  =  , 


*z 

L  2z 


^  (z3 


— ]  ' 
-  1)  -• 


2+  .  e 

.  1  3 


The  integral  of  this  equation  may  be  obtained  by  the  method 

of  variation  of  parameters  .  To  this  end  ,  set 

=  z(z^  -  1)  K(z) 

and  determine  K(z)  by 

dK  ^  _  2£  d  (  z  ^ 

dz  3  dz 


z3  -  1 


so  that 


*2  /  3  ,  >  2  P  2 

z  =  z  (z^  -  Ij  -  u  z 


where  is  an  arbitrary  constant  .-  Therefore  ,  z(x)  is  an 
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elliptic  function  and  w(x)  is  stable  . 

81.  P2  ;  Z1  ;  U2  ,  The  equation  is  (78,1)  ,  where  a  *  h 
=  £  =  0  ,  i.e.  , 


1' '  "  ,  '  >1 

y^  +  (  ^  y+(y-l)|  +k 

L  3y  2(y  -  l)j 

1  y  /  L  y  J 

where  k  is  given  by  (78.3)  . 

We  use  the  notations  of  paragraph  78,  According  to  (79.4)  > 

set 

y  =  3b(y  -  1)  +  3vy 
so  that  see  (79.10)"} 


(81.2) 

y  -  1- 

.  • 

2v-v‘^+2(b+d)v 

By 

a  proper  choice  of  x  —  cp  (x) 

,  one  may  suppose  b  +  4d  —  0 

then  , 

(81.3) 

4  q  *  3  b  ,  r 

+ 

1 

o 

(81.4) 

<! 

II 

< 

1 

9 

(81.5) 

w  =  -  ww  +  w^  +  w 

^  ^2 

(81.6) 

=  2  q  -  q^  =  -  12 

^1 

(81.7) 

•  « 

Cg  =  q  -  q  q  =  12  Vj 

• 

Differentiation  of  (81.6)  and  comparison  with  (61.7)  show 
that  is  a  constant  j  therefore  ,  equation  (79.18)  reduces  to 

(81.8)  V  =  6  +  H  . 
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•  Q 

On  setting  q  =  -  2  ^  ,  one  sees  from  (61,6)  that  Q  is 

Q 

determined  by 

(61.9)  Q  =  3  Q  . 

The  integral  of  (81.1)  is  given  by  (81.2)  and  (8l.4)  ,  where 

(81.10)  w  =*  ^  ; 

V  - 

in  (8l,10)  ,  V  ,  the  solution  of  (81.8)  ,  is  an  elliptic  function 


82,  P2  ;  Z1  ;  U1  .  The  equation  is  (79.1)  with  a=h=0,i.e,, 


(62.1; 

1  y-r  " . H-  11 

y+3y{y-i)r  ^ 

L  3y  2(y  -  1)J 

^  ^  ‘  ^  y"  (y-1 

+k(y  -  1)  +  3  £  y  . 

Using  again  the  method  of  paragraph  79  >  one  finds 

q  +  r  =  0  , 

(82.2)  =  2  q  -  q^  +  2  2  =  -  12  , 

(82.3)  Cg  =  q  -  q  q  +  2  q  »  12 

so  that 

3q-3qq  +  2  £+2£q  =  0. 

Because  =  2  ^  q  see  (79*  9)1  ,  one  also  has 

(82.4)  q-qq+2^q-0 

so  that  Cg  =  0  and  is  a  constant  ;  therefore  ,  equation 
(79.16)  reduces  to 

(62.5)  V  =  6  +  H 
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Now  it  follows  from  (62,2)  and  (62.4)  that 

(62.6)  q  ”  3  q  q  +  -  12  q  »=  0  , 

On  setting  q  *  -  ^  ,  equation  (62,6)  becomes  the  llnea.r 

Q 

equation  of  the  third  order 

(62.7)  Q  -  12  Q 

and  ^  is  given  by 

(62.6)  £ 

where  is  an  arbitrary  constant  , 

The  integral  of  (62,1)  is  given  by 


where  V  ,  the  solution  of  (62.5)  »  is  an  elliptic  function  , 


®3.  PI  ;  Z2  !  U2  .  The  equation  is  (76.1)  with  b=f=k»g=  2  =0, 

i-  •  6  •  $ 

y  =r _ L.  +  ^  1  y^+(ay+d)y  +  3  _  y^(y  -  1) 

L  3y  2(y  -  l)J 

(63.1) 

+  (y  -1)1  a  -  _  -  adl  . 

2  ‘-2 


We  again  use  transformation  (79.3)  and  find 
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(63.2) 

1  3  -V 

y  -  1  * 

*  2 

V  +  2rv  -  V 

with  r  =  - 

(a  +  d)  .  Now  set 

(63.3) 

v  =  w-  q  ,  2q  = 

so  that 

(63.4) 

w=-ww  +  w^  +  CjW  +  C, 

where 

=  -  q  -  =  -  12  , 

•  • 

Cg  =  q  +  2  q  q  =  12  ; 


is  a  solution  of  (79.16)  . 

Therefore  ,  the  integral  of  (63,1)  is  given  by  (63.2-3-4)  > 
where  q  *  ^  is  given  by 

Q 

Q  =  12  e  . 


64.  PI  ;  Z2  ;  U1  .  The  stable  equation  for  this  class  is 

(79.1)  with  b=f~k=0  t  i.e.  , 

2 


(d/f.l)  y  =r  ^  ^  t  y^+(ay+d)y+3y(y-l)  1 

L  3y  2(y-l)  J  L 


y.  3E 


We  use  again  (79.3)  and  find 

(64.2) 


y  .  1  «  3(v^  +  2  £  ) 

2v-v^+2vr-2  £ 


with 

(64.3)  r*-(a  +  d)  ,  0=-2^r. 
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Set 

(64.4)  v»w-q  ,  2q=r 

80  that 

(64.5)  w»-ww4-vP  +  Cj^w+Cg  > 

where 

(64.6)  *  -  q  -  q^  +  2^-  “  -  12 

(84.7)  C2-q+2qq+80q*12V^  ; 

is  a  solution  of  (79.1^)  • 

On  eliminating  £  between  (84.6)  and  (84.7)  »  obtains 


•  o  * 

q  +  6  q  q  +  4  h  “  q  ”  1^  ^  * 

Or  ,  on  setting  2q  «  Q  , 

Q  «  .  3  Q  Q  -  Q  +  24  , 

This  equation  is  of  type  I  and  class  20.4  >  by  the  transfer- 
« 

mation  Q  *  ^  ,  it  reduces  to 

U 

•  9 

U  -  48  U  -  24  V^u  =  0  . 

From  (84,3)  ,  it  follows  that  £  ,  where  is  an 

arbitrary  constant  , 

The  integral  of  (84.1)  is  given  by  (84,2-4)  with  q  *  ^  . 


g5,  PI  :  Z1  ;  U2  .  The  stable  equation  for  this  class  is 

(79.1)  with  £  =  0  ,  i.e.  , 

(85.1)  y  »1  ^  y^+f  ay+  +ci)y+y (y-1 )  | 

2(y-l)4  V  y  '  8  y2  y 
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We  again  use  the  same  method  .  One  has 


(65.2) 


y  -  1  = 


3  v" 


'  T  ' 

2v  - 


V  +  2rv 


On  setting 


(65.3) 


V  =  w  -  q 


t 


one  finds 

(65.4)  w=-ww+w^+C^w+Cg  , 

where 

(65.5)  =  q  -  3  -  (r  +  q  r  -  r^)  =  -  12  , 

(85.6)  Cg  =  q  +  2  q^  +  (  r  +  r  q  -  r^)  q  =  12  ; 

Vj  is  a  solution  of  (79.16)  . 

•  p 

On  eliminating  r  +  q  r  -  r  from  (65.5-6)  ,  one  sees  that  q 
is  a  solution  of  (65.4)  .  Further  ,  on  setting 


U 


in  (65.5)  «  one  obtains 

U  +  q  U  -  (3  -  12  V^)  U  =  0 

which  determines  r  when  q  is  given  . 

The  integral  of  (65.I)  is  determined  by  (65.2-3)  . 
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XII,  Equations  of  the  type  A(x,y)  =  ^ 

Z  {  7 


1  +  1 


y-1  y-H  y 


66.  The  stable  equations  of  this  type  are 


(86.1)  ?  -  il  (  . 

2  \  y  y  -  1  y-H  J  y(y-l) (y-H)  y(y-l) (y-H) 

where  M(x,y)  ,  N(x,y)  are  polynomials  in  y  of  degrees  4  and  6 
respectively  ;  H  is  a  constant  distinct  from  0  and  1  or  H  =  x. 
For  convenience  ,  we  rewrite  equation  (66,1)  as 


.*2 


(86.2)  y  =  JL 
2 


/  ^  ^  ^  )  +  I  ay+a . 

\  y  y-1  y  -  Hj  L 


+y(y-i) (y-H) 


+h+  ^  +  t 


°  ^ 
y-1  y-H  j 

“I 


m 


y  y-1  y  -  Hj 


y  (y-1)  (y-H) 

where  a  ,  a^  ,  ,  m  are  analytic  functions  of  x  . 

The  values  of  y  for  which  the  general  existence  theorem 
does  not  apply  are  y  =  oo,0,l,K, 

If  a  and  h  are  not  both  zero  ,  y(x)  has  simple  parametric 
poles  ;  if  a  =  h  =  0  ,  y(x)  has  double  parametric  poles  . 

g 

If  y(x)  has  a  simple  parametric  pole  ,  set  y  =  _ _  , 

z 

z  =  1  +  u  z  .  According  to  (14.4)  >■  one  has  n  =  -  2  and  hence 
[  see  (16.22)] 

(86.3)  a  *0  ,  h^^O  ,  p  =  0  ,  2hs^  =  1  . 


The  corresponding  condition  for  stability  follows  from 


l6o 


z  u 


+  a. 


(k  +  ^  +  m)  s-fO(z) 


and  is 


+  a. 


-(k  +  £  +m)s 


p 

or  ,  because  2  h  .s  =  1  » 

(66.4)  k  +  ^  +  m  =  0 


0 


and 

(66.5)  ^  +  a^  =  0  or  h  =  2  h  , 


If  y(x)  has  a  double  parametric  pole  ,  one  has  a  »=  h  *  0 
and  thus  k  +  ^!  +  m  ==  0  (  see  ^  17  )  • 

The  values  y  =  0  ,  y  =  1  play  the  same  role  ;  the  trans¬ 
formation  y  =  1  -  w  brings  equation  (66,2)  to  an  equation  of  the 
same  form  according  to  the  equivalence  table 

e  f  g  h  k  p  m 

-f  -e  g  h  -  P  -k  »'iTi 

zero  ,  set  y  =  sz  ,  z=l+uz; 

+  ^  u  +  +  0(z)  , 

s 


(y)  :  H  a  a^  b  c  d 

(w)  :  1-H  -a  a^^+a  c  b  -d 

If  y(x)  has  a  simple  parametric 
then  (66.2)  gives 


where 

A 


-  1  + 


H 


s  .  ,  d  .  ,  s 

+a-+C“  +  b 

2  H  2 

s 


according  to  (66.3)  »  one  has  a  *  0  "j  . 


Now  determine  s  by 

1  +  ^  -  0 


s  s 
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because  p  =  _____  nrust  be  an  integer  ,  one  concludes  that 
s 

b  =  0  ,  e^O  p  p*0  ,  s^  +  EeH=0  ; 

p 

the  condition  for  stability  is  =  0  or  because  s  =  -  2  e  H  , 

(66.6)  -  ^  +  a.  4-  c  -  +  k  ^  =  0  . 

s  H  s 

2 

Because  s  *  -  2  e  H  ,  one  has 

(66.7)  k  =  0  , 

•  • 

(86.6)  2(a  +  c  -  )=  _f__  +  . 

H  e  H 

If  y(x)  has  a  double  parametric  zero  ,  then  b  *  e  *  k  *  0  , 
[  see  ^  17  1  . 

According  to  the  table  of  equivalence  ,  one  obtains  tne 
conditions  for  stability 

for  Ul:  c=0  ,  f  ^  0  , 

(66.9)  0=0  , 

(66.10)  2  f  a  +  b  +  ^  ^  ; 

V-^  1-H/  f  1-H 

for  U2  :  c  =  f  =  0  =  0  . 

67.  We  have  to  consider  two  cases  according  as  to  whether 
H  is  a  constant  or  H  =  x  , 

Suppose  H  is  a  constant  distinct  from  0  and  1  .  The 
transformation  y  »  H(1  -  v)  brings  equation  (66.1)  to  an 
equation  of  the  same  form  according  to  the  table-,  of  equivalence 
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cdefghk  £ 


(v)  ;  1  -aH  aH+a  _ _ ^ 

H  H  H  H 


g  -e  -f  hH^'  -mH  -kH  -Em 


1  Note  that  y  =  0  ,  1  ,  H  correspond  to  v  =*=  1  ,  1  -  ^  ,  0 

^  H 

respectively  ^  . 

Therefore  ,  y  »  H  plays  the  same  role  as  y  =  0  or  y  =  1  . 

The  conditions  for  stability  for  H  follow  readily  from  the 
preceding  paragraph  ;  one  obtains 
for  HI  :  d  =  0  ,  g  0  , 


(87.1) 

(87.2) 


m  =  0 


^  -  -^) 


for  H2  ;  d  =  g  =  m  =  0  , 


Moreover  ,  a  transformation  x  T  (x)  may  be  chosen  so  that 
=  0  £set  tp  =  a^  <p  3  ,  It  then  follows  that  c  ,  f  ,  g  , 
h  are  constant  or  eventually  zero  [see  (66.5-8-10)  ;  (87.2)1  . 

The  stable  equation  of  the  type  considered  is  then 


(87.3)  y  = 


^  ^  ^  \  +y(y- 

y  y-1  y-H  / 


\  +y  (y-1)  (y-H)  f 
>  ^  y^  (y-1 


^  ®  +h 


)^  (y-H)^ 


where  e  ,  f  ,  g  ,  h  are  constant 


To  integrate  this  equation  ,  we  use  the  method  of  variation 
of  parameters  .  One  obtains 

_ 2  _  /-i _ / _  1  \  / - tJT  \  I  V  _L  U -  ®  -L  '  S  1 


y^  *  2y(y-l)  (y-H)  I  K.  +  h  y  -  ^  - 

y  y-1 


y-H 


where  is  an  arbitrary  constant  ;  this  equation  shows  that 


y{x)  is  an  elliptic  function  and  is  stable  . 

The  result  holds  if  one  or  more  of  the  constants  e  ,  f 
g  ,  h  is  zero  , 


Suppose  H  *  X  .  Taking  into  account  the  conditions 
for  stability  already  obtained  for  y  =  oo,  0  ,  1  j^i.e. 
a  =  b  =  c  =  0  ,  k=^=m  =  ,  one  may  rewrite  equation 

(S6.1)  as 


(es.i) 

+y(y-l)  (y-x)!" ^  ^  +  hi  . 

y^  (y  -  1)^  (y  -  x)^ 

Set  y=x+sz  ,z=l+zu;  equation  (S6.1)  yields 

z  u  =  ^  j  (l+s)^+2d(l+s)+2x{x-l)g  I  u  +  A  +  0( 

2  s  z 

where 

Ajj.=  -  2  ^  +  f  +d  |+(2x-l)_f_  + 

(88.2) 

+  ^  '*'-**  r  +  (1  +  s)  ( ^ 

2s  L  s  V  X  X  -  1 

Determine  s  by 

(88.3)  (1+s)^  +  2d(l  +  s)  +  2x(x  -  l)  g  »  o 


and  note  that 
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(ee.4) 


p  = 


1  +  d 


must  be  an  integer  . 

Elimination  of  s  between  (66.3-4)  gives 


(66.5) 


1  + 


^  \  f  1  +  ^ 


I 


+  2x(x  -  1  )g 


Taking  into  account  the  value  of  g  given  by  (66.5) 
rewrites  (66.3)  s-s 

f  s  -  ^  \  ^  s  +  2  +  ^  +  2d 


I, 


SO  that 


1  +  d 


5  q 


2p  +  1 

The  integers  p  »  q  satisfy  the  Diophantlne  equation 
p  +  q  +  2pq  =  0 
so  that  p  =  q  =  0  . 

Accordingly  , 


(66.6) 

(66.7) 


d  -  -  1  , 

s^  =  1  -  2x(x  -  l)g  . 


The  corresponding  condition  for  stability  is  =  0 


(66.6)  2a^.+2(2x-l)g+(l+s‘^) 


V  X  X  -  1 


(66.9) 


^  ^  4-  ^ 

s  X  X  -  1 


=  0  . 


0  . 


one 


or 


From  (66.7-8)  ,  and  (66.9)  >  one  obtains 
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(es.io) 


X  -  1 


»=  0 


3  =  0  . 

Therefore  ,  1  -  2x{x-l)g  or  x(x-l)g  is  a  constant  and 

(ee.ii) 


2x ( X  -  iT 

where  K  is  an  arbitrary  constant  . 

The  values  of  h  ,  e  ^  f  are  now  given  by  (66.5-6-10)  and 

are 


(66.12)  h 


e  = 


x^(x-l)^ 


x(x-l)' 


x^(x-l) 


where  are  arbitrary  constants 

The.  equation  is  then 


*2 

y 


1  + 


y  -  1  y  - 


y  -  X 


X  -•  1 


(66.13 


y(y-i) (y-x) 


2x^'(x  -  1)^ 


r  k,+  lisL  +  . 

^  f-u-  -  1  f  V  ,  xV^ 


k,  ,  k^  ,  k„  ,  k,  are  again  arbitrary  constants  . 

1'2  3  4 

That  the  equation  (66.13)  is  stable  has  been  shown  by 
Palnlev4  .  Equation  (66.13)  is  the  irreducible  equation  VI 
of  Table  I  . 
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